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Results are reported from a joint analysis of Phase I and Phase II data from the Sudbury Neutrino 
Observatory. The effective electron kinetic energy threshold used is Tofi = 3.5 MeV, the lowest 
analysis threshold yet achieved with water Cherenkov detector data. In units of 10^ cm~'^ s~^, the 
total flux of active-flavor neutrinos from *B decay in the Sun measured using the neutral current 
(NC) reaction of neutrinos on deuterons, with no constraint on the *B neutrino energy spectrum, 
is found to be "^nc = 5.140 lQ ;[5g(stat) ^[^ ^^^(syst). These uncertainties are more than a factor of 
two smaller than previously published results. Also presented are the spectra of recoil electrons 
from the charged current reaction of neutrinos on deuterons and the elastic scattering of electrons. 
A fit to the SNO data in which the free parameters directly describe the total *B neutrino flux 
and the energy-dependent Ve survival probability provides a measure of the total ®B neutrino flux 
$SB = 5.046 

-o.i52(st^t) -o'i23(sys*)- Combining these new results with results of all other solar 
experiments and the KamLAND reactor experiment yields best-fit values of the mixing parameters 
of = 34.06 tl'Xl degrees and Amii = 7.59 tJ^ J? x lO"'^ eV^. The global value of $8b is extracted 
to a precision of I2 gs^o- ^ three-fiavor analysis the best fit value of sin^ ^13 is 
This implies an upper bound of sin^ ^13 < 0.057 (95% C.L.). 
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PACS numbers: 26.65.+t, 13.15.+g, 14.60.Pq, 95.85. Ry 



1. INTRODUCTION 

It is by now well-established that neutrinos are massive 
and mixed, and that these properties lead to the oscil- 
lations observed in measurements of neutrinos produced 
in the Sun [l|-p^. in the atmosphere by acceler- 
ators [H, [Tsj . and by reactors [iJI • The mixing model 
predicts not only neutrino oscillations in vacuum, but 
also the effects of matter on the oscillation probabilities 
(the 'MSW effect) To date, the effects of mat- 

ter have only been studied in the solar sector, where the 
neutrinos' passage through the core of both the Sun and 
the Earth can produce detectable effects. The model pre- 
dicts three observable consequences for solar neutrinos: a 
suppression of the i/f, survival probability below the aver- 
age vacuum value of 1 — ^ sin^ 26*12 for high-energy (®B) 
neutrinos, a transition region between matter-dominated 
and vacuum-dominated oscillations, and a regeneration 
of I'eS as the neutrinos pass through the core of the Earth 
(the day/night effect). In addition to improved precision 
in the extraction of the total flux of ®B neutrinos from 
the Sun, an advantage of the low energy threshold anal- 
ysis (LETA) presented here is the enhanced ability to 
explore the MSW-predicted transition region and, in ad- 
dition, more stringent testing of theories of non-standard 
interactions that affect the shape and position of the pre- 
dicted rise in survival probability [T7j-[23|. 

We present in this article a joint analysis of the data 
from the first two data acquisition phases of the Sudbury 
Neutrino Observatory (SNO), down to an effective elec- 
tron kinetic energy of Tes = 3.5 MeV, the lowest analysis 
energy threshold yet achieved for the extraction of neu- 
trino signals with the water Cherenkov technique. The 
previous (higher threshold) analyses of the two data sets 
have been documented extensively elsewhere 0, II] , and 
so we focus here on the improvements made to calibra- 
tions and analysis techniques to reduce the threshold and 
increase the precision of the results. 

We begin in Section [2] with an overview of the SNO 
detector and physics processes, and provide an overview 
of the data analysis in Section |31 In Section 2] we briefly 
describe the SNO Phase I and Phase II data sets used 
here. Section [S] describes changes to the Monte Carlo 
detector model that provides the distributions used to fit 
our data, and SectionlHldescribes the improvements made 
to the hit-level calibrations of PMT times and charges 
that allow us to eliminate some important backgrounds. 

Sections [7]-|9] describe our methods for determining ob- 
servables like position and energy, and estimating their 
systematic uncertainties. Section [10] describes the cuts 
we apply to our data set, while Section [TT] discusses the 
trigger efficiency and Section[T^presents the neutron cap- 
ture efficiency and its systematic uncertainties. We pro- 
vide a detailed discussion of all background constraints 
and distributions in Section [T^ 



Section 113] describes our 'signal extraction' fits to the 
data sets to determine the neutrino fluxes, and Section fTSl 
gives our results for the fluxes and mixing parameters. 



2. THE SNO DETECTOR 

SNO was an imaging Cherenkov detector using heavy 
water (^Il20, hereafter D2O) as both the interaction and 
detection medium [2^. SNO was located in Vale Inco's 
Creighton Mine, at 46°28'30" N latitude, 81°12'04" W 
longitude. The detector was 1783 m below sea level 
with an overburden of 5890 meters water equivalent, 
deep enough that the rate of cosmic-ray muons passing 
through the entire active volume was just 3 per hour. 

One thousand metric tons (tonnes) of D2O was con- 
tained in a 12 m diameter transparent acrylic vessel (AV). 
Cherenkov light produced by neutrino interactions and 
radioactive backgrounds was detected by an array of 9456 
Hamamatsu model R1408 20 cm photomultiplier tubes 
(PMTs), supported by a stainless steel geodesic sphere 
(the PMT support structure or PSUP). Each PMT was 
surrounded by a light concentrator (a 'refiector'), which 
increased the effective photocathode coverage to nearly 
55%. The channel discriminator thresholds were set to 
1/4 of a photoelectron of charge. Over seven kilotonnes 
(7x10® kg) of H2O shielded the D2O from external ra- 
dioactive backgrounds: 1.7 kT between the AV and the 
PSUP, and 5.7 kT between the PSUP and the surround- 
ing rock. Extensive purification systems were used to 
purify both the D2O and the H2O. The H2O outside the 
PSUP was viewed by 91 outward-facing 20 cm PMTs 
that were used to identify cosmic-ray muons. An addi- 
tional 23 PMTs were arranged in a rectangular array and 
suspended in the outer H2O region to view the neck of 
the AV. They were used primarily to reject events not as- 
sociated with Cherenkov light production, such as static 
discharges in the neck. 

The detector was equipped with a versatile calibration- 
source deployment system that could place radioactive 
and optical sources over a large range of the x-z and y-z 
planes (where z is the central axis of the detector) within 
the D2O volume. Deployed sources included a diffuse 
multi-wavelength laser that was used to measure PMT 
timing and optical parameters (the 'laserball') [2^, a 
^®N source that provided a triggered sample of 6.13 MeV 
7s 0], and a ®Li source that delivered tagged /3s with 
an endpoint near 14 MeV In addition, 19.8 MeV 

7s were provided by a '^H(p,7)'*He ('pT') source [2^ and 
neutrons by a ^^^Cf source. Some of the sources were 
also deployed on vertical lines in the H2O between the 
AV and PSUP. 'Spikes' of radioactivity (24Na and ^"^^Rn) 
were added at times to the light water and D2O volumes 
to obtain additional calibration data. Table [T] lists the 
primary calibration sources used in this analysis. 
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Calibration source 


Details 


Calibration 


Deployment Phase Ref. 


Pulsed nitrogen laser 


ooT 'Dcn oon 

oo7, oby, ooo, 


Optical & 


I & II 




('laserball') 


420, 505, 619 nm 


timing calibration 






16 AT 

N 


6.13 McV 7 rays 


Energy & reconstruction 


I & II 


[27] 


»Li 


spectrum 


Energy & reconstruction 


I & II 


[28] 




neutrons 


Neutron response 


I & II 


[25] 


Am-Be 


neutrons 


Neutron response 


II only 




3H(p,7)''He ('pT') 


19.8 MeV 7 rays 


Energy linearity 


I only 


[29] 


Encapsulated U, Th 


0-1 


Backgrounds 


I & II 


[25] 


Dissolved Rn spike 


0-1 


Backgrounds 


II only 




In-situ ^^Na activation 


0-1 


Backgrounds 


II only 





TABLE I: Primary calibration sources. 



SNO detected neutrinos through three processes f30|: 

+ ^ Vx + er (ES) 
I'e + d ^_p + p + e- (CC) 
v,^+d p + + (NC) 

For both the elastic scattering (ES) and charged cur- 
rent (CC) reactions, the recoil electrons were detected di- 
rectly through their production of Cherenkov light. For 
the neutral current (NC) reaction, the neutrons were de- 
tected via de-excitation 7s following their capture on 
another nucleus. In SNO Phase I (the 'D2O phase'), 
the detected neutrons captured predominantly on the 
deuterons in the D2O. Capture on deuterium releases a 
single 6.25 MeV 7 ray, and it was the Cherenkov light 
of secondary Compton electrons or e+e^ pairs that was 
detected. In Phase II (the 'Salt phase'), 2 tonnes of NaCl 
were added to the D2O, and the neutrons captured pre- 
dominantly on ^^Cl nuclei, which have a much larger neu- 
tron capture cross section than deuterium nuclei, result- 
ing in a higher neutron detection efficiency. Capture on 
chlorine also releases more energy (8.6 MeV) and yields 
multiple 7s, which aids in identifying neutron events. 

The primary measurements of SNO are the rates of 
the three neutrino signals, the energy spectra of the elec- 
trons from the CC and ES reactions, and any asymmetry 
in the day and night interaction rates for each reaction. 
Within the Phase I and II data sets, we cannot separate 
the neutrino signals on an event- by- event basis from each 
other or from backgrounds arising from radioactivity in 
the detector materials. Instead, we 'extracted' the sig- 
nals and backgrounds statistically by using the fact that 
they are distributed differently in four observables: effec- 
tive kinetic energy (Tes), which is the estimated energy 
assuming the event consisted of a single electron, cube 
of the reconstructed radial position of the event (R^), 
reconstructed direction of the event relative to the direc- 
tion of a neutrino arriving from the Sun (cos 9q ) , and 
a measure of event 'isotropy' (Pu), which quantifies the 
spatial distribution of PMT hits in a given event (SecO. 
Low values of /3i4 indicate a highly isotropic distribution. 

Figure [1] shows the one-dimensional projections of the 
distributions of these observables for the three neutrino 



signals, showing CC and ES in Phase II and NC for both 
data sets. The Phase II distributions are normalized to 
integrate to 1 except in Fig. [Ijc) , in which the CC and 
NC distributions are scaled by a factor of 10 relative to 
ES for the sake of clarity. The Phase I NC distributions 
are scaled by the ratio of events in the two phases, to illus- 
trate the increase in Phase II. In the figure, and through- 
out the rest of this article, we measure radial positions 
in units of AV radii, so that = {Rf^t/ RavY ■ Figure[2] 
shows the same distributions for some of the detector 
backgrounds, namely 'internal' ^^''Bi and ^os^pj (within 
the D2O volume) and 'AV ^oS'p^ (generated within the 
bulk acrylic of the vessel walls). While some of the ^^''Bi 
nuclei came from decays of intrinsic ^'^®U, the most likely 
source of ^-'^'^Bi was from decays of ^^^Rn entering the 
detector from mine air. The ^os^p^ nuclei came largely 
from decays of intrinsic ^^^Th. Near the Tes = 3.5 MeV 
threshold the dominant signal was from events originat- 
ing from radioactive decays in the PMTs. These events 
could not be generated with sufficient precision using the 
simulation, and so were treated separately from other 
event types, as described in Sec. 113.31 There were many 
other backgrounds; these are described in Sec. fT3l 

The energy spectra provide a powerful method for sep- 
arating different event types. The CC and ES spectra de- 
pend on the shape of the incident neutrino spectrum. We 
treated the CC and ES spectra in two different ways: in 
one fit we made no model assumptions about the under- 
lying spectral shape, allowing the CC and ES spectra to 
vary in the fit, and in a second fit we assumed that the un- 
derlying incident neutrino spectrum could be modeled as 
a smoothly distorted ^B spectrum. The shapes of NC and 
background spectra do not depend on neutrino energy 
and so were fixed in the fit, to within the systematic un- 
certainties derived later. Decays of ^^^Bi and ^osp^j jj^ ^j^g 
detector both led to 7 rays above the deuteron binding 
energy of 2.2 MeV, which created higher energy events 
when the photodisintcgration neutron was subsequently 
captured on either deuterium (Phase I) or predominantly 
■^^Cl (Phase II). A significant fraction of ^^^Bi decays pro- 
duce a 3.27 MeV-endpoint /3. These background events 
are therefore characterized by steeply falling energy spec- 
tra with a photodisintcgration tail, as shown in Fig. Ufa) . 
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FIG. 1: (Color online) The Monte Carlo-generated distribu- 
tions of (a) energy (Tcff), (b) radius cubed (J?^), (c) direction 
(cos^q), and (d) isotropy (/3i4) for signal events. The same 
simulation was used to build multi-dimensional PDFs to fit 
the data. In calculating R^, the radius 7? is first normalized 
to the 600 cm radius of the AV. The CC and NC cos Oq dis- 
tributions are scaled by a factor of 10 for clarity against the 
ES peak. 



FIG. 2: (Color online) The Monte Carlo-generated distribu- 
tions of (a) energy (Toff) on a log scale, (b) radius cubed 
(i?"^), (c) direction (cosSq), and (d) isotropy (/3i4) for back- 
ground events. The same simulation was used to build multi- 
dimensional PDFs to fit the background events. The back- 
grounds shown are internal ^^"Bi, internal ^"^^Tl, and AV 

208rpj 
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CC and ES events produced single electrons and, 
hence, the observed light from these events was fairly 
anisotropic, yielding a correspondingly high value for the 
isotropy parameter, /3i4. The /3i4 distributions show 
small differences due to the different energy spectra of 
the two event types, which affects /3i4 through the known 
correlation between energy and isotropy of an event. The 
isotropy of Phase I NC events looks similar to that of CC 
and ES events, because the 7 ray tended to produce light 
dominated by that from one Compton electron. By con- 
trast, the isotropy distribution of Phase II NC events 
is peaked noticeably lower because neutron capture on 
•^^Cl atoms nearly always resulted in multiple 7s, which 
could each scatter an electron and, hence, produce a more 
isotropic PMT hit pattern. Therefore, (3i4 provides a sen- 
sitive method for separation of electron-like events from 
neutron capture events in this phase, without requiring 
a constraint on the shapes of the CC and ES energy 
spectra, thus providing an oscillation-model-independent 
measurement of the flux of solar neutrinos. The isotropy 
distributions for ^^^Bi events and ^'^^Tl events inside the 
heavy water are noticeably different because, above the 
Teff = 3.5 MeV threshold, Cherenkov light from ^^^Bi 
events was dominated by that from the ground state /3 
branch while that from ^osrpj events was from a /3 and 
at least one additional Compton electron. The differ- 
ence allowed these events to be separated in our fit, as 
was done in previous SNO in-situ estimates of detector 
radioactivity 0, II]- 

The cos 9q distribution is a powerful tool for distin- 
guishing ES events since the scattering of i^e from the 
Sun resulted in electron events whose direction is strongly 
peaked away from the Sun's location. The direction of 
CC events displays a weaker correlation of ~ (1 — |cos 6q) 
relative to the direction of the Sun. The NC distribu- 
tion is flat since the 7s generated by neutron capture 
carried no information about the incident neutrino di- 
rection. Background events had no correlations with the 
Sun's location and, thus, also exhibit a flat distribution, 
as shown in Fig. [2Jc) . 

The radial position of events within the detector yields 
a weak separation between the three neutrino interaction 
types, but a much more powerful level of discrimination 
from external background events. CC and ES events oc- 
curred uniformly within the detector and hence have rel- 
atively flat distributions. NC events occurred uniformly, 
but neutrons produced near the edge of the volume were 
more likely to escape into the AV and H2O regions, where 
the cross section for neutron capture was very high due to 
the hydrogen content. Neutron capture on hydrogen pro- 
duced 2.2 MeV 7s, below the analysis threshold and thus 
less likely to be detected. Therefore, the radial profile of 
NC events falls off at the edge of the volume. This effect 
is more noticeable in Phase I, since the neutron capture 
efficiency on deuterium is lower than on "^^Cl and, hence, 
the neutron mean-free path was longer in Phase I than 
in Phase II. 



3. ANALYSIS OVERVIEW 

The 'LETA' analysis differs from previous SNO analy- 
ses in the joint fit of two phases of data, the much lower 
energy threshold, (which both result in increased statis- 
tics) and significantly improved systematic uncertainties. 

The neutrino signal rates were determined by creating 
probability density functions (PDFs) from distributions 
like those in Figs. [1] and [5] and performing an extended 
maximum likelihood fit to the data. The CC and ES 
spectra were determined by either allowing the fiux to 
vary in discrete energy intervals (an 'unconstrained fit') 
or by directly parameterizing the survival probability 
with a model and fitting for the parameters of the model. 

There were three major challenges in this analysis: re- 
duction of backgrounds, creation of accurate PDFs (in- 
cluding determination of systematic uncertainties on the 
PDF shapes), and extracting the neutrino signals, energy 
spectra, and survival probabilities from the low-threshold 
fits. 

Three new techniques were applied to reduce back- 
grounds compared to previous SNO analyses 0,I1]- First, 
we made substantial improvements to energy reconstruc- 
tion by developing a new algorithm that included scat- 
tered and reflected light in energy estimation. The inclu- 
sion of 'late light' narrowed the detector's effective en- 
ergy resolution by roughly 6%, substantially reducing the 
leakage of low-energy background events into the analysis 
data set by ^60%. Second, we developed a suite of event- 
quality cuts using PMT charge and time information to 
reject external background events whose reconstructed 
positions were within the fiducial volume. Third, we re- 
moved known periods of high radon infiltration that oc- 
curred during early SNO runs and when pumps failed in 
the water purification system. 

Creation of the PDFs was done primarily with a Monte 
Carlo (MC) simulation that included a complete model 
of physics processes and a detailed description of the de- 
tector. We made substantial improvements to the Monte 
Carlo model since our previous publications, and we de- 
scribe these improvements in detail in Sec. [5l 

Our general approach to estimating systematic un- 
certainties on the Monte Carlo-simulated PDF shapes 
was based on a comparison of calibration source data 
to Monte Carlo simulation, as in previous SNO analy- 
ses. In cases where the difference between calibration 
data and simulation was inconsistent with zero, and we 
had evidence that the difference was not caused by a mis- 
modeling of the calibration source, we corrected the PDF 
shapes to better match the data. For example, we ap- 
plied corrections to both the energy (Sec.|S]) and isotropy 
(Sec. HI) of simulated events. Any residual difference was 
used as an estimate of the uncertainty on the Monte Carlo 
predictions. Corrections were verified with multiple cal- 
ibration sources, such as the distributed 'spike' sources 
as well as encapsulated sources, and additional uncertain- 
ties were included to account for any differences observed 
between the various measurements. Uncertainties were 
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also included to take into account possible correlations 
of systematic effects with the observable parameters. So, 
for example, we allowed for an energy dependence in the 
fiducial volume uncertainty, and the uncertainty on the 
energy scale was evaluated in a volume weighted fashion 
to take into account possible variations across the detec- 
tor. 

The final extraction of signal events from the data 
was a multi-dimensional, many-parameter fit. Although 
marginal distributions like those shown in Figs. [T] and [2] 
could be used as PDFs, in practice there are non- 
trivial correlations between the observables that can lead 
to biases in the fit results. We therefore used three- 
dimensional PDFs for most of the backgrounds and for 
the NC signal, factoring out the dimension in cos 6*0, 
which is flat for these events. The CC and ES events 
had PDFs whose dimensionality depended on the type of 
fit. For the unconstrained fit, we used three-dimensional 
PDFs in (i?^/3l4 , cos^q), factoring out the Toff dimen- 
sion because the fit was done in discrete intervals, within 
which the T^s spectrum was treated as flat. For the di- 
rect fit for the Vf, survival probability, we used fully four- 
dimensional PDFs for the CC and ES signals. 

The parameters of the 'signal extraction' fits were the 
amplitudes of the signals and backgrounds, as well as sev- 
eral parameters that characterized the dominant system- 
atic uncertainties. A priori information on backgrounds 
and systematic uncertainties was included. To verify the 
results, we pursued two independent approaches, one us- 
ing binned and the other unbinned PDFs. We describe 
both approaches in Sec. [T31 

We developed and tuned all cuts using simulated 
events and calibration source data. Signal extraction al- 
gorithms were developed on Monte Carlo 'fake' data sets, 
and tested on a 1 /3-livetime sample of data. Once devel- 
oped, no changes were made to the analysis for the final 
fit at our analysis threshold on the full data set. 

In treating systematic uncertainties on the PDF 
shapes, we grouped the backgrounds and signals into 
three classes: 'electron-like' events, which include true 
single-electron events as well as those initiated via Comp- 
ton scattering from a single 7; neutron capture events 
on chlorine that produced a cascade of many 7s with a 
complex branching table; and PMT /3-7 decays, which 
occurred in the glass or envelope of the PMT assembly 
and support structure. The PMT /3-7 events were treated 
separately from other /3-7 events because they were heav- 
ily influenced by local optical effects near the PMT con- 
centrators and support structure, and are therefore hard 
to model or simulate. 

The analysis results presented here have substantially 
reduced uncertainties on the neutrino interaction rates, 
particularly for SNO's signature neutral current measure- 
ment. Although there are many sources of improvement, 
the major causes are: 

• The lower energy threshold increased the statistics 
of the CC and ES events by roughly 30%, and of 
the NC events by - 70%; 



• In a joint fit, the difference in neutron detection 
sensitivity in the two phases provided improved 
neutron/electron separation, beyond that due to 
differences in the isotropy distributions; 

• Significant background reduction due to improved 
energy resolution, removal of high radioactivity pe- 
riods, and new event quality cuts; 

• Use of calibration data to correct the PDF shapes. 

4. DATA SETS 

The Phase I and Phase II data sets used here have 
been described in detail elsewhere 0, II] • We note only a 
few critical details. 

SNO Phase I ran from November 2, 1999 to May 31, 
2001. Periods of high radon in Phase I were removed for 
this analysis based on the event rate. To minimize bias, 
we used Chauvenet's criterion to eliminate runs in which 
the probability of a rate fluctuation as high or higher 
than observed was smaller than 1/{2N), where N is the 
total number of runs in our data set (~ 500). With this 
cut, we reduced the previously published 306.4 live days 
to 277.4. Most of the runs removed were in the first two 
months of the phase, or during a period in which a radon 
degassing pump was known to have failed. This ~9% 
reduction in livetime removed roughly 50% of all ^^^Bi 
events from the Phase I data set. SNO Phase II ran from 
July 2001 to August 2003, for a total of 391.4 live days. 

SNO had several trigger streams, but the primary trig- 
ger for physics data required a coincidence of iVcoinc or 
more PMT hits within a 93 ns window. From the start 
of Phase I until December 20, 2000, A'coinc was set to 
18; it was subsequently lowered to 16 PMT hits. This 
hardware threshold is substantially below the analysis 
threshold, and no efficiency correction was required, even 
at 3.5 MeV (see Secfll])- 

5. MONTE CARLO SIMULATION 

SNO's Monte Carlo simulation played a greater role 
here than in previous publications, as we used it to 
provide PDFs of not only the neutrino signals, but for 
nearly all backgrounds as well. The simulation included 
a detailed model of the physics of neutrino interactions 
and of decays of radioactive nuclei within the detec- 
tor. Propagation of secondary particles was done us- 
ing the EGS4 shower code 31], with the exception of 
neutrons, for which the MCNP 32] neutron transport 
code developed at Los Alamos National Laboratory was 
used. Propagation of optical photons in the detector me- 
dia used wavelength-dependent attenuations of D2O and 
II2O that were measured in situ with laserball calibra- 
tions, and acrylic attenuations measured ex situ. The 
simulation included a detailed model of the detector ge- 
ometry, including the position and orientation of the 
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PSUP and the PMTs, the position and thickness of the 
AV (including support plates and ropes), the size and 
position of the AV 'neck', and a full three-dimensional 
model of the PMTs and their associated light concentra- 
tors. SNO's data acquisition system was also simulated, 
including the time and charge response of the PMTs and 
electronics. Details of the simulation have been presented 
in 0, II]; we describe here the extensive upgrades and 
changes that were made for this analysis. 

Ultimately, SNO's ability to produce accurate PDFs 
depends on the ability of the Monte Carlo simulation 
to reproduce the low-level characteristics of the data, 
such as the distributions of PMT hit times and charges. 
We therefore improved our timing model to more cor- 
rectly simulate the 'late pulsing' phenomenon seen in the 
Hamamatsu R1408s used by SNO. We also added a com- 
plete model of the PMT single photoelectron charge dis- 
tribution that includes PMT-to-PMT variations in gain. 
Gain measurements were made monthly with the laser- 
ball source at the center of the detector, and the simu- 
lation uses different charge distributions for each PMT 
according to these gain measurements. 

Addition of the more complete charge spectrum also 
allowed us to add a detailed model of each electronics 
channel's discriminator. On average, the threshold volt- 
age was near 1/4 of that for a single photoelectron, but 
there were large variations among channels because of 
variations in noise level. Over time, the channel thresh- 
olds were adjusted as PMTs became quieter or noisier; 
these settings were used in the simulation for each run. 
The discriminator model also provided for channel-by- 
channel efficiencies to be included, thus improving simu- 
lation of the detector's energy resolution. 

We made several important changes to the optical 
model as well. The first was a calibration of PMT efhcien- 
cies, which accounted for tube-to-tube variations in the 
response of the photomultipliers and light concentrators. 
These efficiencies are distinct from the electronics dis- 
criminator efhciency described above, as they depended 
on the PMT quantum efficiency, local magnetic field, and 
individual concentrator reflectivity, while the discrimi- 
nator efficiency depended upon PMT channel gain and 
threshold setting. The PMT efficiencies were measured 
using the laserball, as part of the detector's full optical 
calibrations, which were performed once in Phase I and 
three times in Phase II. The efficiencies in the simulation 
were varied over time accordingly. 

The light concentrators themselves are known to have 
degraded over time and the three-dimensional model of 
the collection efficiency of the PMT-concentrator assem- 
bly used in previous analyses had to be modified. We de- 
veloped for this analysis a phenomenological model of the 
effects of the degradation to the concentrator efficiency. 
Rather than modifying the concentrator model itself, we 
altered the PMT response as a function of the position at 
which the photon struck the photocathode. In effect, this 
produced a variation in the response of the concentrator 
and PMT assembly as a function of photon incidence 




Angle of incidence (degrees) 



FIG. 3: (Color online) Comparison of new model of photomul- 
tiplier angular response to data and the old model for Phase I 
at 365 nm. 



angle. A simultaneous fit was performed to laserball cal- 
ibration data at six wavelengths, with each wavelength 
data set weighted by the probability that a photon of that 
wavelength caused a successful PMT hit. The extraction 
of optical calibration data was extended to a larger radius 
than in previous analyses, in order to extract the PMT 
response at wider angles. Ex-situ data were also included 
in the fit to model the response at > 40° for events in the 
light water region. Time dependence was accommodated 
by performing separate fits in time intervals defined by 
the available calibration data: one interval in Phase I and 
three in Phase II. This change improved the modeling of 
any position-dependence of the energy response but did 
not affect the overall energy scale, which was calibrated 
using the ^^N source. We also made a global change to 
the light concentrator reflectivity based on measurements 
with the source. Figure [3] compares the new model 
of the PMT-concentrator response as a function of inci- 
dence angle to that used in earlier publications. 

The laserball calibration data were used as a direct 
input to the energy reconstruction algorithms, provid- 
ing media attenuations, PMT angular response measure- 
ments, and PMT efficiencies. For wavelengths outside 
the range in which data were taken, the Monte Carlo 
simulation was used to predict the response. 

6. HIT-LEVEL CALIBRATIONS 

The accuracy with which we know the charge and time 
of each PMT hit directly affects event position and en- 
ergy uncertainties. To calibrate the digitized charges and 
time, we performed pulser measurements twice weekly, 
measuring pedestals for the charges and the mapping of 
ADC counts to nanoseconds for the times. The global 
channel-to-channel time offsets and the calibration of the 
pulse risetime corrections were done with the laserball 
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source deployed near the center of the detector. These 
calibrations have been described elsewhere Q- 

Four significant changes were made to the calibration 
of PMT charges and times. The first was the removal 
of hits associated with channel-to-channel crosstalk. 
Crosstalk hits in the SNO electronics were characterized 
by having low charges, slightly late times, and being ad- 
jacent to a channel with very high charge. 

The second change was a correction to the deployed po- 
sitions of the laserball source to ensure that the time cal- 
ibrations were consistent between calibration runs. Prior 
to this correction, the global PMT offsets had been sen- 
sitive to the difference between the nominal and true po- 
sition of the source, which varied from calibration run to 
calibration run. The new correction reduced the time- 
variations of the PMT calibrations noticeably, but there 
was a residual 5 cm offset in the reconstructed z-position 
of events, for which a correction was applied to all data. 

There were a variety of ways in which PMTs could fail, 
and we therefore applied stringent criteria for a PMT to 
be included in position and energy reconstruction. The 
criteria were applied to both calibration and 'neutrino' 
data sets as well as to run simulations. 

The last improvement was a calibration to correct for 
a rate-dependence in the electronics charge pedestals. 
Crosstalk hits were used to monitor the pedestal drift 
and a time-varying correction was applied. With this 
correction we could use the PMT charge measurements 
to remove certain types of background events, and to 
substantially reduce systematic uncertainties on the en- 
ergy scale associated with variations in PMT gain, which 
affected the photon detection probability. 

Figure S] shows the distributions of PMT time-of-fiight 
residuals and measured photoelectron charges for a ^^N 
calibration run at the center of the detector compared 
to a simulation of that run. The simulation includes the 
upgrades discussed in Sec. [5j The time residuals show 
excellent agreement in the dominant prompt peak cen- 
tered near At = ns, as well as good agreement for the 
much smaller pre-pulsing (At ~ —20 ns) and late-pulsing 
{At ^ 15 ns and At ~ 35 ns) features. For the charge dis- 
tribution, the agreement is also excellent above 10 ADC 
counts or so, which corresponds to the majority of the 
charges used in the analysis. Thus, we are confident that 
the simulation models the behavior of reconstruction and 
cuts with sufficient accuracy. 



7. POSITION AND DIRECTION 
RECONSTRUCTION 

The primary reconstruction algorithm used in this 
analysis was the same as in previous Phase I publica- 
tions. We used reconstructed event position and direc- 
tion to produce the PDFs shown in Figs. [T] and [H and 
to reject background events originating outside the AV. 
Knowledge of event position and direction was also used 
in the estimation of event energy (see Sec. Below 
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FIG. 4: Comparison of simulation to data for (a) PMT hit 
time-of-flight residuals and (b) photoelectron charge spectra. 



we outline the reconstruction method, and then discuss 
the uncertainties in our knowledge of event positions and 
directions. 



7.1. Reconstruction Algorithm 

The vertex and direction reconstruction algorithm fit- 
ted event position, time, and direction simultaneously us- 
ing the hit times and locations of the hit PMTs. These 
values were found by maximizing the log-likelihood func- 
tion, 

Whit 

log Cire,Ve,te) ^ J^^OgV (t^ , Tf, fe, Ve, t,) , (1) 

with respect to the reconstructed position (re), direction 
(ve), and time (tg) of the event. V{tY^ ,ri;rf.,Ve,te) is the 
probability of observing a hit in PMT i (located at f^) 
with PMT time-of-flight residual if^ (Eq. given a 

single Cherenkov electron track occurring at time te and 
position Te, with direction Ve- The sum is over all good 
PMTs for which a hit was recorded. The PMT time- 
of-flight residuals relative to the hypothesized fit vertex 
position are given by: 



I 5 

c 



(2) 



where ti is the hit time of the ith PMT. The photons are 
assumed to travel at a group velocity with an 
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effective index of refraction averaged over the detector 
media. 

The probabihty V contains two terms to aUow for the 
possibihties that the detected photon arrived either di- 
rectly from the event vertex (Pdircct) or resuhed from 
reflections, scattering, or random PMT noise (T'other)- 
These two probabihties were weighted based on data col- 
lected in the laserball calibration runs. 

The azimuthal symmetry of Cherenkov light about the 
event direction dilutes the precision of reconstruction 
along the event direction. Thus, photons that scattered 
out of the Cherenkov cone tended to systematically drive 
the reconstructed event vertex along the fitted event di- 
rection. After initial estimates of position and direction 
were obtained, a correction was applied to shift the vertex 
back along the direction of the event so as to compensate 
for this systematic drive. The correction varied with the 
distance of the event from the PSUP as measured along 
its fitted direction. 

The reconstruction algorithm returned a quality-of-fit 
statistic relative to the hypothesis that the event was 
a correctly reconstructed single electron. This statistic 
was used later in the analysis to remove backgrounds and 
reduce tails on the reconstruction resolution. Details of 
the reconstruction algorithm can be found in Q. 



7.2. Uncertainties on Position and Direction 

Many effects that could produce systematic shifts in 
reconstructed positions were modeled in the simulation. 
Data from calibration sources deployed within the detec- 
tor were compared to Monte Carlo predictions, and the 
differences were used to quantify the uncertainty on the 
simulation. The observed differences were not deemed 
significant enough to warrant applying a correction to 
the Monte Carlo-generated positions, and so the full size 
of the difference was taken as the magnitude of the un- 
certainty. The differences between data and Monte Carlo 
events were parameterized as four types: 

• vertex offset: a constant offset between an event's 
true and reconstructed positions; 

• vertex scale: a position-dependent shift of events 
either inward or outward; 

• vertex resolution: the width of the distribution of 
reconstructed event positions; 

• angular resolution: the width of the distribution of 
reconstructed event directions relative to the initial 
electron direction. 

These uncertainties can have an impact upon the fiux 
and spectral measurements in two ways: by altering the 
prediction for the number of events reconstructing in- 
side the fiducial volume and by affecting the shape of the 
PDFs used in the signal extraction. 



Reconstruction uncertainties were determined primar- 
ily from ^^N source data. In previous analyses [3], the 
volume density of Compton-scattered electrons relative 
to the source location was modeled with the analytic 
function S{r) ~ exp(^)/(r^). Model improvements for 
this analysis allowed us to extract this distribution for 
each ^^N source run from the Monte Carlo simulation of 
that run, and take into account the exact source geom- 
etry, effect of data selection criteria on the distribution, 
and any time-dependent detector effects. 

The distribution of electron positions was convolved 
with a Gaussian, representing the detector response, and 
the resulting function was fit to the one-dimensional re- 
constructed position distribution along each axis, allow- 
ing both the mean and standard deviation of the Gaus- 
sian to vary for each orthogonal axis independently. An 
example of such a fit is shown in Figure [5] This fit was 
done separately for the ^^N data and the Monte Carlo 
simulation of each ^^N run. The difference in the Gaus- 
sian means gives the vertex offset for that run and the 
square root of the difference in the variances represents 
the difference in vertex resolution. 
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FIG. 5; Fit of the ^''N Compton-electron position distribution 
convolved with a Gaussian to the reconstructed z position of 
'^^N data events for a typical central run in Phase II. 



7.2.1. Vertex Offset 

Analysis of the differences between the reconstructed 
and true event vertex positions at the center of the detec- 
tor, or 'central vertex offset', was done using runs with 
the source within 25 cm of the center, where the source 
position is known most accurately. This avoids confu- 
sion with any position-dependent effects, which are taken 
into account in the scale measurement (Sec. 17.2.21) . A 
data— MC offset was determined for each run, along each 
detector axis. The offsets from the runs were combined in 
weighted averages along each axis, with the uncertainty 
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for each run offset increased to include the uncertainty 
in source position. Ahhough the results showed a small 
mean offset along each axis, the magnitude was com- 
parable to the source position uncertainty and therefore 
we did not correct the PDFs based on this difference. 
Instead, asymmetric double-sided uncertainties were for- 
mulated by using the uncertainty in the weighted aver- 
age, and increasing it by the magnitude of the weighted 
average itself on the side on which the offset was mea- 
sured. The effects of these uncertainties were determined 
during signal extraction by shifting the position of each 
event by the positive and negative values of the uncer- 
tainty along each axis independently, and recomputing 
the PDFs. The values of the uncertainties are given in 
Table in in Sec. 
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7.2.2. Vertex Scale 



FIG. 6: (Color online) Vertex offset along the three detector 
axes as a function of position within the detector. 



A potential position-dependent bias in the recon- 
structed position that can be represented as being pro- 
portional to the distance of the event from the center of 
the detector is defined as a vertex scale systematic. 

In previous SNO analyses, uncertainty in the position 
of the calibration source was a major contribution to re- 
construction uncertainties, especially away from the z- 
axis of the detector, where sources were deployed in a 
less accurate mode. A new method was derived for this 
analysis to reduce sensitivity to this effect. Although the 
absolute source position was known only to ~ 2 cm on 
the z-axis and ~ 5 cm away from this axis, changes in 
position once the source was deployed were known with 
much greater precision. By comparing the result from 
each ^^N run to a run at the center of the detector from 
the same deployment scan, possible offsets between the 
recorded and true source position were removed, thus re- 
ducing source position uncertainties. In addition, any 
constant offset in vertex position, such as that measured 
in Sec. 17.2.11 was inherently removed by this method, 
thus deconvolving the measurement of scale from offset. 
This method allowed data from different scans to be com- 
bined, providing a more representative sampling across 
the time span of the data set and improving the statis- 
tics of the measurement. 

Vertex scale was investigated by using the data— MC 
reconstructed position offset along each detector axis, as 
shown in Figure |6l using only runs within 50 cm of that 
axis to minimize correlations among the three. The runs 
were grouped into 50 cm bins along each axis by source 
position, and the weighted average of the offsets for the 
runs within each bin was found. A linear function was 
fit to the bins as a function of position along that axis. 
Since the method was designed to remove any central 
vertex offset, the function was defined to be zero at the 
center of the detector. 

The slope from the fit provides the scaling required to 
bring the simulation into agreement with data. We did 
not apply a correction, but instead treated it as an asym- 



metric uncertainty on the reconstructed positions of all 
events. The effects observed along the x and y axes were 
of a very similar magnitude and, therefore, were assumed 
to be due to a radial effect, possibly caused either by 
small errors in the modeling of the wavelength-dependent 
refractive index or residual PMT timing calibration er- 
rors. Conservatively, the larger of the x and y values was 
used to bound this effect. The resulting uncertainty was 
applied in our signal extraction fits by multiplying the x, 
y and z position of each event in our PDFs by the value of 
the scale uncertainty, thus shifting events either inwards 
or outwards in the detector, and taking the difference 
from the nominal fit. Since the effect observed along the 
z-axis was larger, the difference of this from the radial 
effect was treated as an additional uncertainty, applied 
only to the z position of events. The values used for each 
uncertainty are listed in Table HIl in Sec. 17.2.61 

Since only runs within 50 cm of each Cartesian axis 
were used to determine vertex scale, diagonal axis runs 
could be used for verification. The method described 
measured the scale for each Cartesian axis independently. 
The values obtained for the y and z axes, for example, 
could therefore be combined to predict the scaling for 
runs on the y-z diagonal. The prediction was shown to 
agree very well with the data, as illustrated in Figure [71 
demonstrating the robustness of the analysis and its ap- 
plicability to events everywhere in the fiducial volume. 

A similar analysis was performed using ■^''^^Cf source 
data in Phase II. The results were consistent with those 
shown here, verifying that the same uncertainties could 
be applied to both electron-like and neutron capture 
events. 

We investigated several other potential causes of vari- 
ation in reconstruction accuracy. The ^^N-source event 
rate during most calibration runs was high in comparison 
to our expected neutrino event rate, so the results were 
checked using low-rate ^^N data. The stability over time 
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FIG. 7: (Color online) Vertex offset along the y-z diagonal 
as a function of position along that diagonal. The dashed 
line shows the prediction from the y- and 2-axis values and 
the solid line shows the best fit scaling value for these data 
points. Observed variations at negative positions are likely 
associated with systematics in source position. 



was determined by comparing runs across the span of 
the two phases. As in previous analyses [7], calibration- 
source dependence was investigated by verifying ^^N re- 
sults using the ^Li source. This also provides a check on 
the energy dependence because the *Li data extended to 
higher energies than the ^^N data. The results were all 
consistent within the uncertainties presented here. 



7.2.3. Vertex Resolution 

The position resolution achieved in this analysis was 
~20 cm for data events. The difference in resolutions 
between data and Monte Carlo events was modeled as 
a Gaussian of standard deviation (or 'width') Coxtra, by 
which the Monte Carlo distribution should be smeared 



to reproduce the data. Coxtra '^^^ given by (cToata ~ "'mc) 
for each ^^N run. This procedure is only valid for 
CMC < CData, which was the likely scenario since any 
minor detector non-uniformities tend to cause a broader 
resolution in the data. In some cases, the simulation and 
data were close enough to one another that statistical 
variation caused CData to appear to be less than (Tmc- In 
these cases, Ko-Q^tg^ — ctmc)! was taken to represent the 
uncertainty in the comparison. The results from the runs 
were combined in a weighted average, independently for 
each detector axis. The resulting values for CTextra are 
listed in Table |ll] in Sec. 17.2.61 These were applied dur- 
ing the signal extraction by smearing the positions of 
all Monte Carlo events by a Gaussian of the appropriate 
width. This was achieved for the binned signal extrac- 
tion (Sec. 114. 2[) by generating a random number for each 
event from a Gaussian of the correct width and adding 
the result to the event's position and, for the unbinned 



method, by a direct analytic convolution (Sec. 114.3)) . 



7.2.4. Angular Resolution 

The source was used for this measurement by re- 
lying on the high degree of colinearity of Compton scat- 
tered electrons with the initial 7 direction. The mean 
of the distribution of reconstructed event positions was 
used to estimate the source position. The reconstructed 
event position was used as an estimate for the scatter- 
ing vertex. To reduce the effect of reconstruction errors, 
only events reconstructing more than 120 cm from the 
source were used. The angle between the initial 7 direc- 
tion (taken to be the vector from the source position to 
the fitted scattering vertex) and the reconstructed event 
direction was found and the distributions of these angles 
were compared for data and Monte Carlo events. 

The same functional form used in previous analyses Q 
was fit to the distributions for data and Monte Carlo 
events within each run. The weighted average of the 
differences in the fitted parameters was computed across 
the runs and the resulting value used as an estimate of 
the uncertainty in angular resolution (given in Table HIl 
SeclLMl)- 



7.2.5. Energy Dependent Fiducial Volume 

The energy dependence of the vertex scaling is of par- 
ticular importance since it could affect the number of 
events that reconstruct within the fiducial volume as a 
function of energy and, hence, distort the extracted neu- 
trino spectrum. Because the ^^N source provided mo- 
noenergetic 7s, giving rise to electrons around 5 MeV, 
whereas the *Li source sampled the full range of the neu- 
trino energy spectrum, the ^Li source was used for this 
measurement. The fraction of events reconstructing in- 
side the source's radial position, closer to the detector 
center, was used as a measure of the number of events 
reconstructing inside the fiducial volume to take into ac- 
count both vertex shift and resolution effects. Absolute 
offsets between data and Monte Carlo events have al- 
ready been characterized in Sections [7.2. lH7.2!3l so a dif- 
ferential comparison of this parameter between data and 
Monte Carlo events was used to evaluate any energy de- 
pendence. A fit from Phase II is shown in Figure [51 The 
energy dependence is given by the slope of a straight line 
fit to the ratio of the data and Monte Carlo parameters, 
averaged across calibration runs. The final uncertainty 
is quoted as an asymmetric, double-sided uncertainty to 
account for the non-zero value of the slope and its uncer- 
tainty. The values for each phase are given in Table [iTj 
The absolute shift, indicated in Fig. [S] by an intercept 
different from one, is a measure of the global vertex scal- 
ing. This effect has already been evaluated in Sec. 17.2.21 
It does not impact the energy dependence and therefore 
is not relevant to this present measurement. 
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FIG. 8: (Color online) Ratio of the fraction of events recon- 
structing inside the source position for data and Monte Carlo 
events, as a function of effective electron energy, for *Li source 
runs. 



An additional check was performed using neutrino data 
from outside the fiducial volume. All standard analysis 
cuts were applied, as described in Sec. [TUl as well as a 
5.5 MeV threshold to select a clean sample of neutrino 
events. A Hill function was fit to the radial distribution 
of the events, with the half-point of the function rep- 
resenting the position of the AV. Statistics in the data 
were limited, so the fit was performed in just three en- 
ergy bins. Monte Carlo simulation of the three types 
of neutrino interactions was combined in the signal ra- 
tios found in a previous SNO analysis [8] and the same 
fit was performed. The ratio of the resulting fitted AV 
position in the data and simulation is a measure of the 
radial scaling and, therefore, the energy dependence of 
this ratio is a check on the analysis described above. 
The results were in good agreement. In Phase II the 
energy dependence was 0.8 ± 2.1%/MeV, in comparison 
to —0.07 ± 0.41%/MeV measured using the ^Li source. 





Uncertainty, Si 


Transformation 


Parameter 


Phase I Phase II 


of observables 


X Offset (cm) 


+ 1.15 
— 0.13 


+0.62 
— 0.07 


X + Si 


y Offset (cm) 


+2.87 
-0.17 


+ 2.29 
-0.09 


y Si 


z Offset (cm) 


+2.58 
-0.15 


+3.11 
-0.16 


Z + Si 


R Scale (%) 


+0.10 
-0.57 


+0.04 
-0.34 


(1 ± Too)*^* 


z Scale (%) 


+0.40 
-0.0 


+0.03 
-0.25 


(1 + Too)^ 


X resn (cm) 


+3.3 


+3.1 


x + Af{0,S,) 


y resn (cm) 


+2.2 


+3.4 


y+M{0,S,) 


z resn (cm) 


+1.5 


+5.3 


z+J^{0,S,) 


Angular resn 


±0.11 


±0.11 


l + (cose0-l)(l + (5O 


EFV (%/MeV) 


+0.85 
-0.49 


+0.41 
-0.48 


W=l + ^{T,s~5.05) 



TABLE II: Systematic uncertainties in the reconstructed po- 
sition and direction of events. EFV is the energy dependent 
fiducial volume uncertainty. The column labeled "Trans- 
formation of observables" refers to the formulae used to 
propagate these uncertainties into the signal extraction fits. 
J\f{0,Si) refers to a convolution with a Gaussian distribution 
of mean 0.0 and standard deviation Si. Events that are pushed 
past cos 00 = ±1.0 are randomly assigned a cosOq value in 
the interval [-1.0, 1.0]. W is an energy-dependent fiducial 
volume factor applied around the midpoint of the ^^N energy, 
where Tcff is the reconstructed effective electron kinetic en- 
ergy and 5.05 MeV is the central Toff value for the ^''N data. 
This was applied as a weight for each event when creating the 
PDFs. ("Resolution" is abbreviated as "resn"). 



8. ENERGY RECONSTRUCTION 

We estimated the kinetic energy of an event after its 
position and direction were reconstructed. The energy 
estimate was used both to reject background events and 
to produce the PDFs shown in Figs. [T] and [2] Improv- 
ing the resolution of the energy estimation algorithm was 
critical because of the low energy threshold of the anal- 
ysis - a 6% improvement in energy resolution reduces 
the number of background events reconstructing above 
threshold by -60%. 



7.2.6. Summary of Reconstructed Position Uncertainties 

Table HIl summarizes the uncertainties in reconstructed 
position and direction. 

It is worth noting that in previous analyses Q the ra- 
dial scaling uncertainty was evaluated at ± 1%, which 
translates to a 3% uncertainty in fiducial volume. The 
improved analysis presented here has reduced the scale 
uncertainty to a little over 0.5% at its maximum and 
significantly less in most dimensions. The resolution dif- 
ferences observed previously were on the order of 9 cm 
[?[, whereas the differences measured here are roughly 
one third of that in most dimensions. The angular res- 
olution uncertainty of 11% is an improvement over the 
16% measured in previous work Q. 



8.1. Total-Light Energy Estimator 

A new algorithm, called "FTK", was designed to use 
all the detected PMT hits in the e nerg y estimate, includ- 
ing scattered and refiected light [33]. The look-up ta- 
ble approach of the prompt-light fitter used in previous 
publications was abandoned in favor of a maximum like- 
lihood method, in which photon detection probabilities 
were generated based on the reconstructed event position 
and direction. The best value of the effective kinetic en- 
ergy, Teff, was found by maximizing the likelihood given 
the observed number of hit PMTs, Vhit, and taking into 
account optical effects due to the reconstructed position 
and direction of the event. In principle, one could con- 
sider a more sophisticated approach in which both the 
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number and distribution of all hit PMTs are used along 
with the recorded time of each hit, but such an approach 
is much more time intensive and was judged to be im- 
practical for the present analysis. 

We considered five sources of PMT hits in an event, 
defined by the following quantities: 

the expected number of detected photons 



• n 



dir 



that traveled directly to a PMT, undergoing only 
refraction at the media boundaries; 



cxp 



the expected number of detected photons 
that were Rayleigh scattered once in the D2O or 
H2O before detection (scattering in the acrylic is 
neglected); 



cxp 



the expected number of detected photons 
that reflected off the inner or outer surface of the 
acrylic vessel; 

• "-cxp " the expected number of detected photons 
that reflected off the PMTs or light concentrators; 

• "-cxp'' " the expected number of PMT noise hits, 
based on run-by-run measurements. 

FTK computed the probabilities of a single photon be- 
ing detected by any PMT via the four event-related pro- 
cesses: pdir, Pscat, Pav, Ppmt • The dircct light probability 
was found by tracing rays from the event vertex to each 
PMT, and weighting each ray by the attenuation proba- 
bility in each medium, transmittance at each boundary, 
solid angle of each PMT, and detection probability given 
the angle of entry into the light concentrator. Scattering 
and reflection probabilities were found using a combi- 
nation of ray tracing and tables computed from Monte 
Carlo simulation of photons propagating through the de- 
tector. 

If Nj is the number of potentially detectable 
Cherenkov photons produced in the event given the in- 
herent PMT detection efficiency, then the expected num- 
ber of detected photons given these probabilities is: 

^cxp dir Pscat Pav ~^ Ppmt )■ (3) 

To be able to compare rioxp to the observed iVhit , we need 
to account for noise hits and convert from detected pho- 
tons to PMT hits, since multiple photons in the same 
PMT produced only one hit. Given the rarity of mul- 
tiple photons in a single PMT at solar neutrino ener- 
gies, FTK made a correction only to the dominant source 



term, n\ 



dir 

cxp 



iV^Pdir- Letting 7Vmpc('^oxp) be the multi 



photon corrected number of direct PMT hits, the total 
expected number of hits is: 

iVcxp(A^^) « ^MPcKtp) 

+iV^ X (pscat + Pav + Ppmt) + ""xp''- (4) 

The probability of observing iVhit hits when iVoxp are ex- 
pected is given by the Poisson distribution: 



P(iVhit|A^7) = 



cxp 



A^hit! 



(5) 



To obtain a likelihood function for Tcff, rather than A^^, 
we integrate over the distribution of given an energy 



xP{N^\T,s)dN.„ 
(6) 

where P{N^ \ T^s) is the probability of Cherenkov 
photons being emitted in an event with energy Teff. The 
negative log-likelihood was then minimized in one dimen- 
sion to give the estimated energy of the event. 



8.2. Energy Scale Corrections and Uncertainties 

We measured the energy scale of the detector by de- 
ploying the tagged 7 source at various locations in 
the x-z and y-z planes within the D2O volume. Al- 
though ^^N was a nearly monoenergetic 7 source, it pro- 
duced electrons with a range of energies through multi- 
ple Compton scattering and e^e~ pair production. As 
a result, the single 6.13 MeV 7 produced an 'effective 
electron kinetic energy' (Toff) distribution that peaked 
at approximately 5 MeV. 

Using the ^^N 7-ray source to determine the detector's 
energy scale is complicated by its broad spectrum of elec- 
tron energies. To separate the detector's response from 
this intrinsic electron energy distribution, we modeled 
the reconstructed energy distribution with the integral 



P(Teff) = N / 

-^source 



1 



~dE, 



(7) 
Pi + 
is the ap- 



where A'^ is a normalization constant, (j{E^- 
P2\/E^- is the detector resolution, and Psource 
parent electron energy distribution from the ^^N 7 rays 
without including the detector optical response, sets 
the displacement of the ^^N peak, and therefore the offset 
in energy scale at that source location. The Psource dis- 
tribution was computed from a Monte Carlo simulation 
of 7 propagation through the source container and pro- 
duction of Cherenkov photons from Compton-scattered 
and pair-produced e^e' . We translated the number 
of Cherenkov photons in each simulated event to a most 
probable electron (MPE) kinetic energy with the same 
tables that were used in the FTK energy estimation al- 
gorithrn, and generated the distribution, Psource, of event 
values [33| . Given this fixed distribution for the ^^N cali- 
bration source, we fit for A'', pi, P2, and p^ in each source 
run, for both data and for Monte Carlo simulation of the 
same source position and detector state. The parameter 
differences between data and Monte Carlo, run-by-run, 
determined the energy corrections and uncertainties. Pa- 
rameters pi and p2 measure the detector energy resolu- 
tion, and are discussed further in Sec. 18.31 Parameter 
was used here to define the spatial energy scale correction 
and uncertainties. 

The Monte Carlo was initially tuned by adjusting a 
global collection efficiency parameter in the simulation to 
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minimize the difference between data and Monte Carlo 
energy scales for ^^N runs at the center of the detector. A 
series of additional corrections were then applied to the 
estimated energy of all the data and Monte Carlo events, 
to remedy known biases. 

Approximations in FTK's handling of multiple hits on 
a single tube lead to a small energy non-linearity, and 
we derived a correction for this by comparing the recon- 
structed energy for Monte Carlo events to their true en- 
ergies. Similarly, the simple PMT optical model used by 
FTK produced a small radial bias in event energies and, 
again, comparison of reconstructed energies of Monte 
Carlo events to their true values were used to provide 
a correction. 

Two additional corrections were based on evaluations 
of data. The first was to compensate for the degradation 
of the PMT light concentrators, which changed the de- 
tector's energy response over time during Phase I. The 
degradation affected the fraction of light that was re- 
flected off the PMT array. We tracked the variation using 
^^N runs taken at the center of the detector, and created 
a time-dependent correction to event energies that shifted 
their values by up to 0.4% [33]. 

The final correction was applied to remove a variation 
in energy with the detector z-coordinate. Figure 
shows the difference between the average reconstructed 
energies of events from the ^^N source for each calibra- 
tion run, and the Monte Carlo simulation of the run, as 
a function of the radial position of the source. As can be 
seen, for events in the top (positive z) hemisphere of the 
detector, the Monte Carlo underestimated the event en- 
ergies by as much as 3% and, in the bottom hemisphere, 
it overestimated the energies by almost the same amount. 
The cause of the former was the simulation's poor optical 
model of the acrylic in the neck of the AV. The latter was 
likely caused by accumulation of residue at the bottom 
of the acrylic vessel and variations in the degradation of 
the PMT light concentrators. 

To correct for the z-dependence of the energy scale, we 
first split the ^^N calibration runs into two groups. One 
group contained runs on the x-z plane along with half 
of the runs on the z-axis, and was used to construct the 
correction function. The second group contained runs on 
the y-z plane along with the other half of the z-axis runs, 
and was used later to independently evaluate the spatial 
component of the energy scale uncertainty. 

We found that the variation in the energy scale best 
correlated with the vertical position of the event (z) and 
the direction cosine of the event relative to the z-axis 
{uz). AH of the ^^N events in the first group were binned 
in the (z, Uz) dimensions and the peak of the energy 
distribution was found for data and Monte Carlo events 
separately. We fit a second-order polynomial in z and Uz 
to the ratio of the data and Monte Carlo peak energies. 
This smooth function provided the spatial energy correc- 
tion for data events. Fig.[3Jb) shows the spatial variation 
after this energy correction. 

To evaluate the spatial component of the energy scale 
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FIG. 9: (Color online) DifFerence between data and 
Monte Carlo energy scales as a function of radius for Phase II 
source runs in the upper hemisphere, on the equatorial 
plane, and in the lower hemisphere. Panel (a) shows the sig- 
nificant variation in these three regions before the spatial en- 
ergy correction. Panel (b) shows the same runs after the spa- 
tial energy scale correction is applied. (The fiducial volume 
cut is at 550 cm). 



uncertainty, we assumed azimuthal symmetry in the de- 
tector, and divided the second group of ^^N calibration 
runs into regions based on radius and polar angle. Within 
each region, the RMS of the individual run differences be- 
tween the corrected data and Monte Carlo energy scales 
defined the uncertainty on the energy scale in that vol- 
ume. All regions were then combined into a volume- 
weighted measure of the uncertainty on the overall en- 
ergy scale in the detector due to spatial variation and 
non-uniform sampling of the detector volume. As a ver- 
ification of the procedure, we reversed the roles of the 
two calibration groups (using the y-z plane to construct 
the calibration function and the x-z plane to evaluate 
the uncertainties) and found very similar corrections and 
uncertainties. 

The energy scale uncertainty of the detector also in- 
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eludes uneertainty in modeling of energy loss in the ^^N 
source itself, uncertainties in the online status of PMTs, 
variation in the channel response between high-rate cali- 
bration data and low-rate neutrino data, and uncertain- 
ties in the data acquisition channel gains and thresholds, 
which affect the photon detection probability. Many of 
these uncertainties have been substantially reduced com- 
pared to previous publications by the improvements to 
the Monte Carlo model described in Sec. [5] and the rate- 
dependent correction to the channel pedestals described 
in Sec. ID 

The components of the energy scale uncertainties are 
summarized in Table Hill We take the source uncertainty 
as 100% correlated between phases, and the other uncer- 
tainties as uncorrelated. To verify the validity of the ^^N- 
derived energy corrections and uncertainties over a wider 
range of energies, we compared the data and Monte Carlo 
energy distributions for ■^^■^Cf neutron source runs and 
the D2 0-volume radon spike, for both of which events 
are more widely distributed in the detector than for the 
^^N source. In both cases, the agreement between the 
data and Monte Carlo was well within the uncertainties 
stated in Table Hill 



Uncertainty Phase I Phase II 

PMT Status ±0.01% ±0.01% 

Threshold/Gain ±0.18 -0.31% ±0.13 -0.07% 
Rate ±0.3% ±0.05% 

Source ±0.4% ±0.4% 

Spatial Variation ±0.18% ±0.31% 

Total ±0.56 - 0.62% ±0.52 - 0.51% 



TABLE III: Summary of energy scale uncertainties. 



8.3. Energy Resolution 

Energy resolution was a significant systematic uncer- 
tainty because of its impact on background acceptance 
above the 3.5 MeV energy threshold. Due to differing 
event topologies in the two phases, the resolution uncer- 
tainties were treated as three independent, uncorrelated 
systematic parameters: Phase I events (both electron- 
like and neutron capture events). Phase II electron-like 
events, and Phase II neutron capture events. In all cases, 
the resolution was found to be slightly broader in the 
data than for Monte Carlo events. The difference was 
parameterized as a Gaussian of width (loxtra, with which 
the Monte Carlo distribution was convolved to reproduce 
the data. The width of the Gaussian was given by the 
quadrature difference of the data and Monte Carlo resolu- 
tions: CTcxtra = \J ("'Data ~ '''mc)' ^ resolution corrcctiou 
was formulated using calibration source data and applied 
to the Monte Carlo events used in PDF generation. The 
uncertainties on this correction were then taken from the 



spread of the calibration data. 

8.3.1. Energy Resolution Uncertainties for Phase II 
Electron-like Events 

The ^^N source was the primary source for this mea- 
surement. We evaluated the uncertainties in two ways 
by measuring the resolution for the spectrum of Comp- 
ton electrons differentially and integrally. 

The MPE fit described in Sec. 18.21 unfolds source ef- 
fects from the event distribution, allowing the extraction 
of the intrinsic monoenergetic electron resolution as a 
function of energy. The fit was performed for both data 
and Monte Carlo simulation of ^^N runs and the resulting 
resolutions were compared differentially in energy. The 
energy resolution at threshold is the dominant concern 
for electron-like events, due to the exponential rise of the 
backgrounds, and the value at 3.5 MeV was therefore 
used as representative of the detector resolution. (Tcxtra 
at threshold was found to be 0.152 ± 0.053 MeV. In terms 
of the fractional difference: 



CTMC 

this translates to atrac =2.4 ± 1.6% at threshold. 

To measure the integrated Compton electron resolu- 
tion using the monoenergetic 7 rays produced by the ^^N 
source, the reconstructed energy distribution for Monte 
Carlo-simulated 7s was convolved with a smearing Gaus- 
sian and the result was fit directly to the data, allow- 
ing the mean and width of the smearing Gaussian to 
vary. The resulting aextra of the smearing Gaussian was 
0.0 ± 0.046 MeV. This measurement represents a higher 
average energy than the 'unfolded' MPE value since the 
-'^^N provides 7s at 6.13 MeV. The value of fTftac from this 
7-ray measurement is 0.00 ± 0.08%. 

Two ^^^Rn spikes were deployed during Phase II, one 
in the D2O and one in the H2O volume. These provided 
a low energy source of /3s and 7s, below the analysis 
threshold and, therefore, all observed decays appeared 
due to the detector energy resolution, making the spikes 
particularly sensitive to this effect. The unbinned sig- 
nal extraction code (Sec. 114.31) was used in a simplified 
configuration to fit the data from each spike. 

The internal spike was fit with 3 PDFs in two di- 
mensions: energy and isotropy. The PDFs were ^^"^Bi 
electron-like events (primarily /3s) in the D2O volume, 
^^^Bi photodisintegration neutrons, and a 'quiet' data 
set drawn from neutrino runs near the date of the spike. 
The latter provides the energy distribution of all 'back- 
ground' events to the spike measurement, including other 
radioactive decays such as PMT /3-js as well as neutrino 
interactions. An analytic convolution parameter was also 
floated, defining the width of the convolving Gaussian 
applied to the Monte Carlo electron-like events. The re- 
sulting (Toxtra was 0.139 ^Q^oag MeV, which is equivalent 
to cTfrac =2.0 ± 1.0% at threshold. Floating the ^"Bi 
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electrons and neutrons independently also allowed a ver- 
ification of the Monte Carlo prediction for the photodis- 
integration rate. The results were in good agreement, 
giving 0.91 ± 0.13 times the Monte Carlo predicted rate. 

The external spike was fit with two PDFs in just the 
energy dimension, due to lower statistics. The elec- 
tron to neutron ratio in the ^^''Bi PDF was fixed to the 
Monte Carlo prediction and the overall normalization of 
this PDF was taken as a free parameter, along with the 
quiet data normalization. The Monte Carlo events were 
again convolved with a Gaussian, whose width was al- 
lowed to vary in the fit. The resulting value for (Textra 



was 0.273 



+0.030 
-0.035 



MeV, which gives crft-ac =7.6 ± 1.9% at 



threshold. The broader resolution for external events, 
which were generated in the H2O region but either trav- 
eled or were misreconstructed into the D2O, is not unex- 
pected since the detector's energy response was modeled 
less well in the outer detector regions. 

These four measures were combined to give the resolu- 
tion correction and associated uncertainty for electron- 
like events in Phase II. Since the two measure- 
ments are not independent, they were not used together. 
The weighted mean of the MPE fit and the two spike 
points was used to give the correction, with an associ- 
ated uncertainty. The difference of that value from the 
weighted mean of the "'^^N 7 point and the two spike 
points was then taken as an additional one-sided (neg- 
ative) uncertainty, to take into account the difference 
in the two ^^N measurements. This results in a final 
value of CToxtra = 0.168lg'Q8o MeV, which was applied as 
a constant smearing across the energy range. The four 
measurements and the resulting one sigma band on the 
final correction value for Phase II electron-like events are 
shown in Figure [TUl 

The MPE fit was also applied to the ^Li source but 
this was not included in the calculation due to the low 
statistics of the measurement. However, the energy de- 
pendence of both the *Li and the ^^N MPE fits were used 
to demonstrate that the use of a constant (Xoxtra across 
the energy spectrum was consistent with the data avail- 
able. 
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FIG. 10: (Color online) Measurements of energy resolution 
in Phase II. The solid area shows the one sigma band on the 
energy resolution correction applied to Phase II electron-like 
events. The ^^^Cf and muon follower points show the mea- 
surements of the energy resolution for neutron capture events, 
and were not used to evaluate the total shift for electron-like 
events. 



Gaussian and fitting the resulting form to the data. The 
mean and width of the convolving Gaussian were allowed 
to vary, in order to take into account possible correlations 
between energy scale and resolution effects. The result 
was (Textra = 0.153 ± 0.018 MeV. The observed energy 
scale from this measurement agreed very well with that 
evaluated in Sec. 18.21 

The statistics of the muon follower data set were low, 
and the resulting uncertainty on the measurement was 
therefore relatively large. Nevertheless, a similar analysis 
was performed, giving a Ccxtra of 0.237 ± 0.144 MeV. 

The weighted mean of the two points was used for the 
final correction to the energy resolution of neutron cap- 
ture events in Phase II, with its associated uncertainty, 
with the value dominated by the ^^^Cf measurement: 
Cextra = 0.154 ± 0.018 MeV. Both points are also shown 
on Fig. [ini 



8.3.2. Energy Resolution Uncertainties for Phase II 
Neutron Capture Events 

The energy resolution for neutron capture events in 
Phase II was measured using the ^^^Cf source, with a 
verification performed using a 'muon follower' data set, 
consisting of neutron capture events occurring within a 
defined time window after a muon passed through the 
detector. 

There are fewer uncertainties associated with the neu- 
tron measurement since the ^^^Cf source produced neu- 
trons whose captures on ^^Cl and deuterium resulted 
in the same 7 cascades as those from NC events. The 
measurement was performed by numerically convolving 
a spline-fit of the Monte Carlo energy distribution with a 



8.3.3. Energy Resolution Uncertainties for Phase I 
Electron-like Events 

No radon spikes were deployed in Phase I, and so only 
the two ^^N measurements were available. Both the MPE 
fit and the Gaussian convolution to the 7-ray energy dis- 
tribution were performed for Phase I ^^N runs, in the 
same manner as for Phase II (Sec. I8.3.ip . The central 
correction value was taken from the MPE fit directly, 
giving (Textra = 0.155 ± 0.036 MeV. The small number 
of energy resolution measurements in Phase I provides 
fewer handles on the uncertainty than the much-better 
calibrated Phase II. The uncertainties in Phase I were 
therefore chosen to match those of Phase II. The width 
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of the convolving Gaussian for Phase I events was there- 



fore taken as tTp 



0.155 



+0.041 
'0.080 



MeV. This was also 



applied to neutron capture events in Phase I, since the 
event topologies were similar. 



8.4. Energy Linearity 



The corrections derived in Sec. 18.21 were done primar- 
ily using the ^^N source, and therefore the uncertainty 
in the energy scale at the ^^N energy is very small. An 
additional uncertainty was included to account for pos- 
sible differential changes in the energy scale that were 
not correctly modeled in the Monte Carlo simulation. 
Such changes could be caused by residual crosstalk hits, 
or mis-modeling of the multi-photon PMT hit probabil- 
ities in the energy reconstruction algorithm. The differ- 
ential changes were determined relative to the point, 
and used calibration sources whose energies were sub- 
stantially higher. 

The pT source provided 7s roughly 14 MeV higher in 
energy than those from ^^N, resulting in a good lever-arm 
on any non-linear effects. This source was only deployed 
in Phase I since deployment in Phase II would have re- 
sulted in an overwhelming neutron signal. The difference 
between data and Monte Carlo-reconstructed event en- 
ergies was measured to be —1.36 ± 0.01% at the energy 
of the pT source. 

The MPE fit described in Sec. 18.21 was applied here 
to the *Li source, including an additional term in the 
parameterization to model first-order differential changes 
in the energy scale. The fit was done to both data and 
Monte Carlo events, and a difference of just —0.011 ± 
0.004% was found, evaluated at the same energy as the 
pT source 7 rays. 

Giving the pT and ®Li sources equal weight, the aver- 
age shift in energy scale at the energy of the pT source 
was found to be —0.69%. Using this as a measure of the 
degree by which the Monte Carlo energy scale could vary 
differentially from the data and assuming a linear inter- 
polation between the ^^N and pT energies, the linearity 
uncertainty was parameterized in terms of the difference 
of an event's energy from the ^^N source (~5.05 MeV). 
This results in a scaling factor that can be applied to the 
energy of each Monte Carlo event used to build the PDFs 
in the signal extraction procedure. Conservatively, this 
was applied as a two-sided uncertainty: 



T' 

off 



1.0 ±0.0069 X 



5.05 



19.0-5.05 



Tcff, 



(9) 



where 19 MeV is the effective energy of the pT source. 
Toff is the original effective kinetic energy of an individual 
event and T^jj is the modified energy. 

Tests using both the ^Li and ^^^Cf sources suggested 
no evidence for any linearity shift in Phase II. We expect 
any source of linearity shift to be common across the two 
phases, however, and therefore the results from Phase I 



were conservatively taken to apply to both phases in a 
correlated fashion. 



9. EVENT ISOTROPY 

As discussed in Sec. [21 we used a measure of event 
'isotropy' as one dimension of our PDFs to help distin- 
guish different types of events. By 'isotropy' we mean 
the degree of uniformity in solid angle of the hit PMTs 
relative to the fitted event location. 

Single electron events, like those created in neutrino 
CC and ES reactions, had a Cherenkov cone that, at solar 
neutrino energies, was somewhat diffuse due to electron 
multiple scattering in the water. Nevertheless, even with 
the multiple scattering, these events were characterized 
by a fairly tight cluster of PMT hits in a cone aligned 
with the forward direction of the electron. 

Neutron capture events on deuterium in Phase I led 
to a single 6.25 MeV 7 ray. Although these events could 
produce multiple Compton electrons and, hence, a num- 
ber of Cherenkov cones that distributed hits more widely 
than single electrons. Phase I neutron capture events in 
the data set were dominated by single Compton scatters 
and, thus, isotropy was not useful in distinguishing them 
from CC or ES events. 

In contrast, in Phase II neutrons captured primarily on 
•^^Cl, which typically led to a 7 cascade that looks very 
different from single electrons. Neutron capture on '^^Cl 
typically produced several 7 rays, with energies totaling 
8.6 MeV, which distributed PMT hits more uniformly in 
solid angle. The isotropy distribution for these events is 
thus a convolution of the isotropy distribution of single 
7-ray events with the directional distribution of the 7 
rays emitted in the possible 7-decay cascades. 

The isotropy of background events can also be signifi- 
cantly different from that of single electron and neutron 
events. Decays of ^'^^Tl, for example, produce both a /3 
and a 2.614 MeV 7 ray and, thus, resulted in a different 
distribution of hit PMTs than either single electrons or 
single 7s. 

The measure of isotropy was therefore critical to the 
analysis, helping us to separate CC and ES events from 
NC events, and both of these from low-energy back- 
ground events. 

We examined several measures of isotropy, including 
a full correlation function, the average angle between all 
possible pairwise combinations of hit PMTs, and con- 
structions of several variables using Fisher discriminants. 
We found that, for the most part, they all had compa- 
rable separation power between the single electron (CC 
and ES) and the neutron (NC) signals. As in our pre- 
vious Phase II publications @, we opted to use a linear 
combination of parameters, /3i4 = Pi -\- 4/34, where: 



N-l N 



N{N - 1) 



(10) 



= 1 J=i+1 
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FIG. 11: /3i4 isotropy distributions for ^^^Ci data and MC and 
data and MC. There is a very small shift of the Monte 
Carlo ^^^Cf /3i4 distribution toward higher (less isotropic) val- 
ues. 

In this expression, P; is the Legendre polynomial of order 
I, Oij is the angle between triggered PMTs i and j relative 
to the reconstructed event vertex, and N is the total 
number of triggered PMTs in the event. Very isotropic 
events have low (even negative) values of /3i4. 



Phase II, we treated the uncertainties on the mean and 
width of the /3i4 distribution for NC events and photo- 
disintegration neutrons separately from the other event 
types. Uncertainties on the /3i4 distributions of /3s and 
7S from radioactive background events were treated the 
same as for CC and ES events. The one exception to this 
was PMT /3-7 events, whose location at the PMT array 
led to effects on the /3i4 distribution that are not present 
in the other signals. The /3i4 distribution and associated 
uncertainties for PMT fi-^s are discussed in Sec. 113.31 

As usual in this analysis, we derived uncertainties on 
the mean, width, and energy dependence of the /3i4 dis- 
tribution by comparing calibration source data to Monte 
Carlo simulations of the calibration source runs. When 
we found a difference that was corroborated by more than 
one source, or was caused by known errors in the simula- 
tion, we adjusted the simulated distribution by shifting 
the mean of the distribution and/or convolving the dis- 
tribution with a smearing function to better match the 
calibration data. In such cases, additional uncertainties 
associated with the correction were included. 



9.1.1. /3i4 Uncertainties for Phase II Electron-like Events 



9.1. Uncertainties on the Isotropy Measure 

We parameterized the difference between the predicted 
(3i4 PDF and the true PDF by a fractional shift in the 
mean, /3i4, and a broadening of the width, ap-^^. We also 
allowed for an energy dependence in the shifts. 

Figure [TT] shows /3i4 distributions of Phase II data from 
^^^Cf and ^^N sources and from corresponding MC simu- 
lations. The ^^N source emitted a single 6.13 MeV 7 ray, 
which usually underwent Compton scattering and pro- 
duced one or more electron tracks, while neutrons from 
the ^^^Cf source were typically captured in Phase II by 
the chlorine additive, leading to a cascade of several 7 
rays. It is clear from the figure that the data and 
Monte Carlo agree very well, while the Monte Carlo simu- 
lation of the ^^^Cf source shows a very small shift toward 
higher /3i4 values (less isotropic events than in the data). 
This shift is discussed in Sec. 19.1.31 

Errors in the simulated distributions of /?i4 can have 
several sources: incorrect modeling of the detector optics 
or photomultiplier tubes, unmodeled event vertex recon- 
struction errors, errors in the model of the production 
of Cherenkov light (including the interactions of 7 rays 
and electrons in the detector) and, for neutrons captured 
on -^^Cl, uncertainties in our knowledge of the 7 cascade 
sequences and correlations between the directions of the 
multiple 7 rays. 

Except for the last item, these errors affect all event 
types. For Phase I, in which neutrons were captured 
on deuterons, we allowed for correlations among the un- 
certainties on all signals and most backgrounds. For 



The primary measure of isotropy uncertainties for 
Phase II electron-like events comes from comparisons of 
^^N calibration source data to Monte Carlo simulation. 
We fit Gaussians to both the data and simulated events 
for each run, and calculated the fractional difference be- 
tween the fitted parameters. Figure [12] shows the frac- 
tional difference in the means as a function of R^. Each 
point shown is the fractional difference for a single run, 
with the error bar evaluated as the combination of the 
uncertainty on the fit parameters for data and Monte 
Carlo events. The detector region in which the source 
was deployed has been identified for each run. 

Also shown in Fig. [12] are the averages of these differ- 
ences, in several radial bins. The uncertainty on each av- 
erage is the standard deviation of the points in that bin, 
weighted by the volume represented by the bin (smaller 
volumes have larger uncertainties). The overall weighted 
average within the entire 550 cm radius fiducial volume 
is consistent with zero, with an uncertainty of ±0.21%. 
The calibration data were collected at a high rate relative 
to normal neutrino data runs and, so, we added to this an 
uncertainty to account for the difference in /3i4 between 
high rate and low rate data (±0.1%) by comparing low- 
rate and high-rate ^^N source runs, as well as a small 
uncertainty of ±0.002% associated with a possible un- 
modeled time-dependence obtained by comparing data 
and Monte Carlo differences over time. The quadrature 
combination of these uncertainties on the mean of the 
/3i4 distribution totals ±0.24%. A similar analysis was 
performed for the width of the /3i4 distribution, yielding 
a total fractional uncertainty of ±0.54%. 
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FIG. 12: (Color online) Fractional differences in the mean 
of the /3i4 distributions for data and Monte Carlo, for the 
Phase II ^^N calibration source. Also shown in the figure are 
the averages in each radial bin, with the bands indicating the 
volume-weighted uncertainty in each bin. 



FIG. 13: (Color online) Fractional differences in the mean 
of the /3i4 distributions for data and Monte Carlo, for the 
Phase I calibration source. Also shown in the figure are 
the averages in each radial bin, with the bands indicating the 
volume-weighted uncertainty in each bin. 



9.1,2. /3i4 Uncertainties for Phase I Electron-like Events 

We applied an identical analysis to the Phase I ^^N 
data but, as shown in Figure [T31 we found a difference 
of —0.81 ± 0.20% between the means of the /3i4 distribu- 
tions for source data and source simulations. Compar- 
ison of ^^N data between Phase I and Phase II showed 
them to be consistent, and the (data-Monte Carlo) differ- 
ence seen in Fig. [T3]to be due to a shift in the simulated 
events. Further investigation showed that the difference 
was caused by the value of the Rayleigh scattering length 
used in the Phase I simulation. Explicit measurements of 
the Rayleigh scattering had been made and used in the 
simulation for Phase II but no such measurements ex- 
isted for Phase I. Use of the Phase II Rayleigh scattering 
length in Phase I simulations was found to produce the 
desired magnitude of shift, and we therefore corrected 
the /3i4 values of all simulated Phase I events by a factor 
of (1 - 0.0081) = 0.9919. 

We included three uncertainties associated with this 
correction. The first was 0.20% on the correction itself, 
evaluated from the volume-weighted average of the data 
and Monte Carlo differences for Phase I, as shown in 
Fig. [131 To take into account the fact that we used the 
consistency in the ^^N data between the two phases to 
support the correction of —0.81%, we added in quadra- 
ture the uncertainty on the difference between the means 
of the Phase I and Phase II ^^N /3i4 distributions, which 
was 0.34%. Finally, because we used the consistency of 
the Phase II data with the Monte Carlo simulation as evi- 
dence that the Phase I /3i4 distribution was correct, aside 
from the Rayleigh scattering correction, we included the 
volume- weighted Phase II uncertainty on the offset of the 
mean (0.21% from Fig. [H] in Sec. igXTj) . 

The evaluations of the uncertainties associated with 



rate dependence and time dependence in Phase I were 
0.08% and 0.03%, respectively, and the overall uncer- 
tainty on the mean of the j3n distribution in Phase I 
thus totaled 0.42%. 

We evaluated the uncertainty on the width of the fin 
distribution for Phase I in the same way as for Phase II, 
finding a fractional uncertainty which also totaled 0.42%. 



9.1.3. /3i4 Uncertainties for Phase n Neutron Capture 
Events 

Neutron capture events in Phase II were distinct from 
other neutrino-induced events and backgrounds in that 
the 7 cascade was more isotropic than a single electron or 
7 ray. The primary measurement of the uncertainty on 
the mean of the /3i4 distribution comes from deployments 
of the ^^^Cf source, which produced several neutrons per 
fission decay. The /3i4 distribution of the resulting neu- 
tron capture events was noticeably non-Gaussian, and we 
therefore derived uncertainties on the mean and width by 
fitting the distributions from simulated ^^^Cf runs 
directly to the distributions of data. The fit allowed for 
scaling as well as convolution with a Gaussian smearing 
function. Figure [Til shows the fit of a simulated ^^^Cf 
run to data, in which the fitted scaling was -1.2% and 
the smearing was an additional 1.8% of the width of the 
Monte Carlo distribution. 

We derived scaling factors from fits like that in Fig. [141 
for all ^^^Cf runs, and then volume- weighted them in 
the same way as for the ^^N data. The average of the 
volume-weighted differences showed an overall offset be- 
tween the means of the /3i4 distributions for data and 
Monte Carlo of ^ —1.4%. This result was not consis- 
tent with that from the ^^N data for Phase II (which. 
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FIG. 14: (Color online) Fit of Monte Carlo simulated ^14 
distribution for neutron capture events from ^^^Cf to data 
taken with the ^®^Cf source. The fitted shift for this sam- 
ple is —1.2%, and the additional smear is 1.8%, before any 
corrections for bias. 



FIG. 15: (Color online) Fractional difference in mean /3i4 
between data and Monte Carlo events for several neutron 
sources. The horizontal band indicates the error on the overall 
— 1.44% correction. 



as discussed above, had no significant offset), which in- 
dicated that the shift was not due to a detector effect. 
To check whether the shift was caused by mis-modeling 
of the ■^''^^Cf source in the simulation, we performed the 
same analysis on several types of neutron capture events: 
neutrons produced by passage of a muon through the de- 
tector ('muon followers'), neutrons from a tagged Am-Be 
source, and neutrons produced by deuteron photodisin- 
tegration during the deployment of a radon spike in the 
detector. Figure [T5] shows results from these sources. An 
energy-dependent fit to all sources except ^^^Cf showed 
an offset of —1.12 ± 0.31%, consistent with the data 
from the ^^^Cf source. This indicated that the offset 
was likely not a source effect but was instead associated 
with the simulation of the 7 cascade from neutron cap- 
tures on chlorine, possibly with some contribution from 
the energy- dependent correction of the Monte Carlo value 
for /3i4 presented in Sec. 19.1.51 All sources taken together 
gave an overall offset of —1.44%, and we therefore cor- 
rected the /3i4 PDF by multiplying each simulated event's 
/3i4 value by (1 + Sp.J = (1 - 0.0144) = 0.9856. 

The uncertainties on this correction came first from 
the uncertainty on the overall average, which was 0.17%. 
To this we added in quadrature the same rate- and 
time-dependent uncertainties as were calculated for the 
Phase II ^^N sources. We also added an uncertainty as- 
sociated with the multiplicity of neutrons from the ^^^Cf 
source of 0.09% (neutrons produced by either photo- 
disintegration of deuterons or the NC reaction are sin- 
gles, whereas the ^^^Cf source produces multiple neu- 
trons/decay), and 0.03% uncertainty to account for the 
relatively sparse sampling of the detector, giving a to- 
tal of 0.22%. Conservatively, we included a further un- 
certainty based on the difference between ^^^Cf and the 
other neutron-source data, a one-sided uncertainty of 
0.31%. The total uncertainty on the mean of the /3i4 



distribution for Phase II neutron captures was therefore 

+0.38 (w 
-0.22 '0. 

As well as a measure of any required shift, the fit de- 
scribed above also allowed for the widths of the data and 
Monte Carlo distributions to differ. A resolution param- 
eter was varied in the fit, as the standard deviation of the 
Gaussian by which the Monte Carlo distribution was an- 
alytically convolved. The results for each ^^^Cf run were 
volume- weighted using the procedure described above to 
result in an average overall smearing value. The same 
fit was performed on a sample of Monte Carlo-generated 
data, and the bias determined from these fits was sub- 
tracted from the overall average. The result was a frac- 
tional smearing correction to be applied to the PDFs 
of 0.43%, with an uncertainty (including all sources de- 
scribed above: time, rate, multiplicity, and sampling) of 
0.31%. 



9.1.4- /3i4 Uncertainties for Phase I Neutron Capture 
Events 

Neutrons created in Phase I captured on deuterons, 
releasing a single 6.25 MeV 7 ray. The uncertainties on 
the mean and width of the /3i4 distribution were therefore 
well-estimated by the measurements made with the "'^^N 
6.13 MeV 7-ray source, already discussed in Sec. 19.1.21 
We therefore used the same uncertainties for both event 
types, applied in a correlated fashion. 



9.1.5. Energy Dependence of I3\4, Uncertainties 

A final systematic uncertainty on the /3i4 distributions 
is their energy dependence. In Figure [12] we show the 
energy dependence of the fractional difference between 
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FIG. 16: (Color online) Fractional shift in mean /?i4 between 
data and Monte Carlo simulation in Phase II for several cali- 
bration sources as a function of kinetic energy, with the fit to 
Eq. shown. 



Monte Carlo predictions of the mean of the P14 distribu- 
tion and data from several different sources: the Phase II 
radon spike, low and high energy ^^N source events, the 
^^^Cf source (with the data corrected by the 1.44% shift 
discussed above) , and ^Li-source /? events in three energy 
bins. There clearly is an energy dependence in the data, 
which we fit with a function of the form: 

/ = Sfj,, + mp,,{T,s - 5.6 MeV), (11) 

where Tcs is kinetic energy, and 5.6 MeV is the kinetic 
energy at the high-energy -'^^N point (the point used to 
determine the offset in the mean of the Phase II electron 
/3i4 distribution). With this parameterization, the offset 
(Sfj^^) and the slope (TO/314) uncorrelated. Given that 
all the sources exhibited the same trend, we applied the 
same slope to all event types, but used the different off- 
sets and uncertainties for dp-^^ described in the previous 
sections. We performed a similar analysis for Phase I, 
although less calibration data were available, and found 
that the same slope fit the ^^N and ^Li data in this phase. 

We found no energy dependence in the broadening of 
the width of the distributions. These uncertainties 
were therefore treated as independent of energy. 

The corrections and uncertainties to the /3i4 distribu- 
tions are listed in Tables IIVI and |Vl 



Phase/Particles 


'5/314 


mpi4 (10- 


^ MeV"^) 


Il/electrons 


0.0 ±0.0024 


2.76 ± 


0.696 


Il/neutrons 


U.Ui** _0.0022 


2.76 ± 


0.696 


I/electrons 


-0.0081 ±0.0042 


2.76 ± 


0.696 


I / neutrons 


-0.0081 ±0.0042 


2.76 ± 


0.696 



TABLE IV: Summary of uncertainties on the (3-i4 scale. The 
/?i4 of each event was corrected by: /?i4 — >■ Pi4{l + (Spj^^ + 
m^i4(Teff -5.6 MeV))). 



Phase/Particles Correction (%) Uncertainty (%) 



Il/electrons 


0.0 


±0.42 


Il/neutrons 


0.43 


±0.31 


I/electrons 


0.0 


±0.42 


I/neutrons 


0.0 


±0.42 



TABLE V: Summary of uncertainties on the P14 width. 



10. CUTS AND EFFICIENCIES 

The data set contains two main types of background 
events: physics backgrounds, due to radioactive decays, 
and instrumental backgrounds, caused by the detector 
itself. Two sets of cuts were developed to remove these 
events, described in Sections 110.11 and 110.21 Each set 
of cuts had an associated level of signal loss, which was 
taken into account in the measurement of neutrino flux 
and spectra as described in Sec. 110.51 



10.1. Low-Level (Instrumental) Cuts 

There were many sources of instrumentally-generated 
events in the SNO detector, which produced hits origi- 
nating either in the PMTs or in the electronics channels. 
Static discharges in the nitrogen in the neck of the acrylic 
vessel and 'flasher' PMTs, in which discharges occurred 
within a photomultiplier tube itself, produced light in 
the detector. Electronic pickup generated by noise on 
the deck above the detector or by high-voltage break- 
down could produce hits in electronics channels. We re- 
moved these instrumental backgrounds with a suite of 
loose 'low-level' cuts that rejected events before event 
reconstruction. The cuts were based on event charac- 
teristics such as the distribution of PMT hit times, the 
presence of unusually low or high PMT charges, or un- 
usual time correlations between events (such as bursts of 
events with large numbers of hits) . More details on these 
low- level cuts can be found in 0, Si ■ We used the same 
cuts and cut criteria here, with the exception that the 
simple burst cut used in |7| was not used in this analysis 
because it was redundant with other burst cuts. 

The acceptance of these cuts was re-evaluated for this 
analysis, particularly in the low-threshold region (below 
Teff = 5.0 MeV) where the cuts had not previously been 
examined in detail. We discuss the results of these cut 
acceptance measurements in Sec. 110.5.11 



10.2. High-Level Cuts 

Background radioactivity events were produced pri- 
marily by the decays of ^^^Bi and 208'pj_ Lower energy 
(Teff < 3 MeV) decays of these nuclei in the heavy water 
could appear above our T^s = 3.5 MeV threshold because 
of the broad energy resolution intrinsic to a Cherenkov 
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detector. Decays within the waUs of the acryhc vessel, 
the light water surrounding the vessel, and the photomul- 
tiplier tube array could pass the energy cut and have mis- 
reconstructed vertex positions which falsely placed them 
within the fiducial volume. The PMT array was, by far, 
the radioactively hottest component of the SNO detec- 
tor and, consequently, the largest source of background 
events. We designed a suite of 13 loose cuts that used 
'high-level' information (reconstructed event position, di- 
rection, and energy) to remove events whose likely ori- 
gin was either outside the fiducial volume or whose true 
energy was below our threshold. All of the cuts were 
adjusted based exclusively on simulated events and cal- 
ibration data. Several of the cuts had a high degree of 
redundancy in order to maximize background rejection. 
The acceptance of the cuts was therefore evaluated col- 
lectively, as described in Sec. 110.51 

Five of the high-level cuts removed backgrounds using 
Kolmogorov-Smirnov (KS) tests of the hypothesis that 
the event had a single Cherenkov-electron track. Two 
of these tests compared azimuthal and two-dimensional 
(polar vs azimuthal) angular distributions to those ex- 
pected for Cherenkov light produced by an electron, and 
two others did the same for hits restricted to a narrow 
prompt time window. The fifth of these KS tests was a 
comparison of the distribution of fitted PMT time residu- 
als (see Eq. ([2])) with the expected distribution for direct 
Cherenkov light. 

Three more of the cuts applied event 'isotropy' to re- 
move misreconstructed events. Events whose true ori- 
gins were well outside the fiducial volume but which re- 
constructed inside tend to appear very anisotropic. For 
one of these cuts we used the mean angle between pairs 
of PMTs, (Oij), and for another the isotropy parameter 
/3i4, which is described in Sec. HI Both of these have been 
used in previous SNO analyses 0, [1] . The third of these 
cuts was based on the charge-weighted mean pair angle, 
Oij, in which each pair angle is weighted by the product 
of the detected charges of the two PMTs in the pair. 

Further cuts used information from the energy recon- 
struction algorithm discussed in Sec. 18.11 Two cuts re- 
moved events whose reported energy uncertainty was well 
outside of the range expected from the known energy res- 
olution. These are referred to in Sections ri0.5.3l - [T0.5.5l as 
the 'energy-uncertainty' cuts. The third was a compar- 
ison of the energy estimated with FTK (which used all 
hits) with that from a prompt-light-only energy estima- 
tor. Events whose origins were outside the acrylic vessel 
and which pointed outward often had a larger fraction 
of prompt hits because the direct light was not attenu- 
ated by the acrylic vessel. Such an event would have a 
higher energy as measured by a prompt-light energy es- 
timator than by the total-light energy reconstruction of 
FTK. We normalized the ratio of these two energy esti- 
mates by the ratio of prompt to total hits in the event. 
The cut itself was two-dimensional: events were removed 
if the normalized ratio of energy estimates was unusually 
large and the charge-weighted 9ij was unusually low (the 



latter indicating an outward-pointing event with a tight 
cluster of hits) . 

The last two high-level cuts were also used in deter- 
mining the PDFs for radioactive backgrounds from the 
PMTs. The first of these, the in-time ratio (ITR) cut, re- 
moved events based on the ratio of the prompt hits to the 
total hits. The prompt time window for the ITR cut ex- 
tended from 2.5 ns before the reconstructed event time 
to 5.0 ns after, and the full-event window was roughly 
250 ns long. The mean of the ITR distribution for SNO 
events is at 0.74. Events that were reconstructed at po- 
sitions far from their true origin tend to have small ITR 
values, because the PMT hits were spread across the en- 
tire time window. In previous analyses 0-11] we used 
the ITR cut with a fixed threshold, rejecting events with 
an in-time ratio smaller than 0.55. For the lower-energy 
events included in this analysis, the lower number of hits 
caused the distribution of ITR to broaden and introduced 
a large, energy- dependent bias in the acceptance of the 
cut. We therefore changed the cut threshold to scale 
with the number of hits (iVhit) in an event. The fixed 
value of 0.55 used in earlier publications corresponded to 
cutting events that fell more than 2.7a below the mean 
of the distribution, and we retained this criterion, so 
that the new version of the ITR cut rejected events that 
were more than 2.7a below the mean of 0.74, where now 
a = O.iS/VN^t- 

The last cut was aimed directly at removing events 
produced by radioactive decays in the PMTs themselves. 
Such events produced light either in the PMT glass or 
in the light water, just in front of the PMTs. Although 
only a tiny fraction of such events were misreconstructed 
inside the fiducial volume, the PMT array was relatively 
hot, with a total decay rate from uranium and thorium 
chain daughters of a few kHz. Because of their origin 
within or near the PMTs, these events were character- 
ized by a large charge in one PMT (or distributed over 
a few nearby PMTs) with hit times that preceded the 
reconstructed event time. The 'early charge' (EQ) cut 
therefore examined PMT hits in a window that ran from 
— 75 ns to —25 ns before the event time. If a PMT hit in 
this window had an unusually high charge, or there was 
an unusually large number of hits in this window, then 
the event was cut. To account for variations in PMT 
gain, 'unusually high charge' was defined by using the 
known charge spectrum of the PMT in question to cal- 
culate the probability of observing a charge as high as ob- 
served or higher. If more than one hit was in the window, 
a trials penalty was imposed on the tube with the lowest 
probability, and an event was cut if this trials-corrected 
probability was smaller than 0.01. We defined 'unusually 
large number of hits' in a similar way, by comparing the 
number of hits observed in the early time window to the 
expected number, given the total number of hits in the 
event. If the Poisson probability of having the observed 
number in the early time window was below 0.002, the 
event was cut. 
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10.3. Burst Removal 

Atmospheric neutrinos, spontaneous fission, and 
cosmic-ray muons could all produce bursts of events 
that were clearly not due to solar neutrinos. Most 
of these bursts had a detectable primary event (like a 
high-energy atmospheric-neutrino event) followed by sev- 
eral neutron events. In addition, many instrumentally- 
generated events came in bursts, such as those associated 
with high- voltage breakdown in a PMT. 

We therefore applied several cuts to the data set to re- 
move most of these time-correlated events. Four of these 
were part of the suite of instrumental cuts described in 
Sec. 110. ll The first removed events that were within 5 fis 
of a previous event and, therefore, eliminated events as- 
sociated with PMT afterpulsing or Michel electrons from 
decays of stopped muons. The second removed all events 
within 20 seconds of an event that had been tagged as 
a muon. Most of these 'muon followers' were neutrons 
created by passage of a cosmic-ray muon through the 
heavy water, which captured either on deuterons or, in 
Phase II, on '^^Cl, but the cut also removed longer- lived 
cosmogenic activity. The muon follower cut resulted in 
a very small additional overall detector deadtime be- 
cause of the very low rate of cosmic rays at SNO's depth. 
Atmospheric neutrinos could also produce neutrons, ei- 
ther directly or by creating muons which, in turn, dis- 
integrated deuterons. We therefore removed any event 
within 250 ms of a previous event that had iVhit > 60 
(Phase I) or A^hit > 150 (Phase II). The fourth cut was 
aimed primarily at residual instrumental bursts, and re- 
moved events that were part of a set of six or more with 
A'hit > 40 that occurred within an interval of six seconds. 

Because of the relatively loose criteria used, after these 
cuts were applied there were still time-correlated events 
in the SNO data set that were very unlikely to be so- 
lar neutrinos, but were primarily low-multiplicity neu- 
trons created by atmospheric neutrino interactions. We 
therefore applied a final 'coincidence cut' that removed 
events if two or more occurred within a few neutron cap- 
ture times of each other. For Phase I this window was 
100 ms; a shorter window of 15 ms was used for Phase II 
because of the shorter neutron capture time on chlorine 
compared to deuterium. The cut was 'retriggerable', in 
that the window was extended for its full length past the 
last event found. If a new event was thus 'caught', the 
window was again extended. We calculated that this cut 
removed less than one pair of events from each data set 
due to accidental coincidences. 



10.4. Cut Summary 

The numbers of events in the data sets after successive 
application of each set of cuts are shown in Table IVll The 
burst cuts described in Sec. 110.31 are included in instru- 
mental cuts, except for the final coincidence cut, which 
appears in the last line of the table. 



Events Phase I Phase II 

Full data set 128421119 115068751 

Instrumental 115328384 102079435 

Reconstruction 92159034 77661692 

Fiducial volume (<550 cm) 11491488 8897178 

Energy range (3.5-20 MeV) 25570 40070 

High-level cuts 9346 18285 

Coincidence cut 9337 18228 



TABLE VI: Number of events remaining in the data set after 
successive application of each set of cuts. 



10.5. Cut Acceptance 

As in previous analyses Q , the fraction of signal events 
expected to pass the full set of analysis cuts (the 'cut ac- 
ceptance') was determined by separating the cuts into 
three groups: instrumental, reconstruction, and high- 
level cuts. Correlations between these groups had been 
shown to be minimal [8] , and it was verified that this was 
still true after the addition of new high-level cuts for this 
analysis. 

The ^^N and ^Li calibration sources were used for the 
primary measurements of cut acceptance and the ^'"'^Cf 
source was used for neutron capture events in Phase II. 
Neutron events in Phase I are well-modeled by ^^N events 
since capture on deuterium resulted in a single 7 at 
6.25 MeV and "'^^N was a source of 6.13 MeV 7s. 



10.5.1. Instrumental Cut Acceptance 

The instrumental cuts were not simulated in the Monte 
Carlo code and, therefore, we could not make a relative 
estimate of their acceptance by comparing simulation to 
data. Instead, an absolute measure of their acceptance 
was made using calibration data and applied as a correc- 
tion (with uncertainties) to the PDFs. 

Being a near-perfect source of CC-like electron events, 
the ^Li source was used to evaluate the signal loss for 
electron-like events, and ^^^Cf was used for Phase II neu- 
tron capture events. The ^^N source was used as a check 
and any difference in the values obtained was conserva- 
tively taken as a two-sided systematic uncertainty. Fig- 
ure [T7] shows the ^^N and ^Li measurements in Phase I. 
The weighted mean of the ^Li signal loss shown in the 
figure was taken as the correction to the PDFs, and the 
median deviation of the points from this value was used 
to represent the energy-dependent uncertainty. 

The ^^N source, which was deployed more frequently 
and at more positions than ^Li, was used to determine 
time- and position-dependent uncertainties. Runs were 
binned by position and date, and the median deviation 
of the bin values from the best-fit value was taken as the 
measure of systematic uncertainty. 
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1.0 was taken as the value of the uncertainty. The ^Li 
source was used to investigate energy dependence. As ex- 
pected, the signal loss decreased at higher energies, where 
more information was available to reconstruct an event. 
The simulation was shown to reproduce this effect very 
accurately and the uncertainty was therefore treated in 
the same manner as the position-dependent uncertainty. 

Combining the systematic uncertainties in quadrature, 
we obtained the final uncertainties associated with recon- 
struction acceptance: 

• Phase I: ±0.034% (stat) ±0.060% (syst) 

• Phase II e^: ±0.037% (stat) ±0.090% (syst) 

• Phase II n: ±0.000% (stat) ±0.009% (syst) 



FIG. 17: (Color online) Signal loss due to the instrumental 
cuts for the ^''N and *Li calibration sources as a function of 
reconstructed kinetic energy, in Phase I. 



After combination of the systematic uncertainties in 
quadrature, the final estimates of signal loss due to the 
instrumental cuts were: 

• Phase I: 0.214% ±0.026 (stat) ±0.094 (syst) 

• Phase II e": 0.291% ±0.028 (stat) ±0.202 (syst) 

• Phase II n: 0.303% ±0.003 (stat) ±0.186 (syst) 

where "e~" refers to electron-like events and "n" to neu- 
tron captures. The acceptance is given by one minus the 
fractional signal loss and was applied as an adjustment 
to the normalization of the PDFs. 



10.5.3. High-Level Cut Acceptance 

To take into account the acceptance of the high-level 
cuts, the ratio of the cut acceptance for data and Monte 
Carlo events was calculated and applied to the PDFs as 
a normalization correction. This ratio was evaluated as 
a function of energy, position and time. 

The energy-uncertainty cuts described in Sec. 110.21 
were observed to have much stronger variations in signal 
loss as a function of position and energy than the other 
high-level cuts and were therefore treated separately. It 
was verified that the correlations between the two re- 
sulting subsets of high-level cuts were minimal, so that 
treating them independently was a valid approach. The 
following sections describe the analysis for each subset of 
cuts, where 'reduced high-level cuts' refers to the subset 
that does not include the energy-uncertainty cuts. 



10.5.2. Acceptance of Reconstruction 

Occasionally, the reconstruction algorithm failed to 
converge, and returned no vertex for an event. In past 
analyses, an upper bound was placed on the resulting sig- 
nal loss, by using calibration source data, but a different 
approach was used in this analysis. What is important 
is how well the effect is reproduced in the simulation. 
Therefore, a comparison was made of the acceptance of 
data and Monte Carlo events and the difference of the 
ratio from unity was taken as a systematic uncertainty 
on the PDF normalization. 

Results from the ^^N source, and the ■^''^^Cf source for 
Phase II neutrons, demonstrated that the signal loss in 
the data was reproduced by the simulation to within the 
statistical uncertainties. Analysis of runs taken during 
the two phases showed no significant deviation with time. 
A position-dependent uncertainty was evaluated by tak- 
ing the ratio of the acceptance of ^^N data and Monte 
Carlo events as a function of source deployment position. 
The difference of the weighted average of the points from 



10.5.4. Reduced High-Level Cut Acceptance 

The data/Monte Carlo acceptance ratio and its uncer- 
tainty were calculated for each calibration source run. 
The runs were divided into radial bins, and the error- 
weighted mean and standard deviation were calculated 
in each bin. Finally, the volume-weighted average of the 
bin values was calculated. 

The energy dependence of the acceptance ratio was in- 
vestigated using ^^N and ^Li data for electron-like events 
and ^^^Cf for Phase II neutron capture events. The ^^N 
data were restricted to the energies below 9 MeV to avoid 
complications associated with event pile-up caused by the 
high rate of the calibration source. 

The measurements from ^^N and ^Li were in very good 
agreement, and were both consistent with the acceptance 
ratio having no dependence on energy. The normaliza- 
tion correction for the PDFs was therefore evaluated us- 
ing the -'^^N source data by taking the weighted mean of 
the values in each energy bin. The median deviation of 
the ^Li points from the best-fit was taken as a systematic 
uncertainty on the energy-dependence. 



25 



The acceptance ratio for Phase II neutron capture 
events was evaluated using ^^^Cf data. To avoid pile- 
up of fission 7s, the events were required to have energies 
in the interval 4.5-9.5 MeV. An energy-dependent un- 
certainty was included to account for any variation of 
individual energy bins from the overall average. 

The stability of the acceptance as a function of time 
was studied using ^^N runs taken in the center of the 
detector. No trend was observed, but the time variability 
was incorporated as an additional systematic uncertainty. 

The ^^N source was also used to evaluate a system- 
atic uncertainty associated with a possible position de- 
pendence of the acceptance ratio. Runs were binned by 
position in the detector, the volume-weighted average of 
the bins was found and the mean deviation of the ratio 
in each bin from this value was calculated. A comparison 
of ^^N and ^'^^Cf source data showed that they exhibited 
statistically equivalent position dependences, so the more 
widely deployed ^^N source was used to quantify this ef- 
fect for both electron-like and neutron capture events. 

The acceptance corrections and associated uncertain- 
ties derived from the difference between the high-level cut 
acceptances for data and Monte Carlo events are summa- 
rized in Table EUl 





Phase I Phase II e 


Phase II n 


Correction 


0.9945 


0.9958 


0.9983 


Stat uncert (%) 


0.0273 


0.0159 


0.0196 


Energy dep (%) 


0.1897 


0.1226 


0.0005-2.3565 


Position dep (%) 


0.1630 


0.3144 


0.3144 


Time dep (%) 


0.0805 


0.0130 


0.0130 



TABLE VII: Correction and associated uncertainties for the 
high-level cut acceptance ratio. The Phase II neutron energy- 
dependent uncertainty was treated differentially with energy; 
the quoted range covers the value across the energy spectrum. 



position-dependent uncertainties using the more widely- 
deployed ^^N source. -^^N source data were divided into 
50 cm slices along the z-axis, and the acceptance ra- 
tios calculated in the slices were combined in a volume- 
weighted average. The uncertainty on this average was 
derived from the deviation of the points from unity. 

The energy-uncertainty cuts were even more sensitive 
to the effects of pile-up than were the other high-level 
cuts. Therefore, to evaluate an energy-dependent uncer- 
tainty on the acceptance ratio for electron-like events, 
events from the ^^N source were restricted to energies 
below 7 MeV, and the lower rate ^Li source was used for 
measurements at higher energies. ^^^Cf data were used 
for Phase II neutron capture events, with the deviations 
from unity measured in the 8.5-9 MeV bin also applied to 
higher energy events. This resulted in energy-dependent 
uncertainties for both electron-like and neutron capture 
events. 

The uncertainties in acceptance were applied as uncer- 
tainties in normalization of the PDFs. The values are 
summarized in Table I Villi 





Phase I 


Phase II e 


Phase II n 


Stat uncert (%) 


0.0377 


0.0668 


0.0322 


Position dep {+) 


(%) -H0.0750 


-H0.0838 


-H0.0838 


Position dep ) 


(%) -1.0760 


-0.9897 


-0.9897 


Time dep (%) 


0.0834 


0.0531 


0.0531 



TABLE VIII: Uncertainties on the energy- uncertainty cut ac- 
ceptance ratio. Energy-dependent uncertainties were treated 
differentially with energy and are not shown. The uncertainty 
in position is asymmetric. 



10.5.6. Overall Cut Acceptance 



10.5.5. Energy- Uncertainty Cut Acceptance 

We expect that the effect of placing cuts on the uncer- 
tainty on the estimate of an event's energy reported by 
the energy reconstruction algorithm should be the same 
for data and Monte Carlo events. Nevertheless, uncer- 
tainties on this assumption were evaluated using the ^®N 
and ^^^Cf source data, applying the same energy ranges 
as in the reduced high-level cut analysis fSec. 110. CT]) . 

Differential uncertainties were evaluated using the 
same method as for the reduced high-level cuts. The 
stability over time was measured using ^^N data. The 
acceptance ratio was observed to be stable, but an ad- 
ditional uncertainty was included based on the spread of 
the points. 

The ^^N and ^^^Cf data showed statistically equiv- 
alent position-dependent behavior in the acceptance of 
the energy-uncertainty cuts, and we therefore evaluated 



The final correction to the PDF normalization comes 
from combination of the high-level cut correction (Ta- 
ble lVIII) and the instrumental cut correction fSec. lI0.5TT|) . 
The various contributions to uncertainty on signal loss 
were treated as uncorrelated and combined in quadra- 
ture to give the final uncertainty on the cut acceptance 
correction. Table HXl lists the final corrections and uncer- 
tainties. 



Phase I Phase II e Phase II n 
Correction 0.9924 0.9930 0.9954 

Pos uncertainty (%) 0.34-0.45 0.41-0.80 0.38-2.70 

Neg uncertainty (%) 1.12-1.17 1.07-1.08 1.06-1.65 



TABLE IX: Corrections applied to the Monte Carlo-generated 
PDFs due to cut acceptance. The uncertainties were evalu- 
ated differentially with energy; the quoted range covers their 
values across the energy spectrum. 
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Figure [TH] shows a comparison of the cut acceptance 
for data and Monte Carlo events from a single ^^^Cf run 
in Phase II. The full set of analysis cuts was applied to 
both data and simulation, and the Monte Carlo-predicted 
acceptance was corrected by the value from Table HXl As 
the figure shows, the Monte Carlo simulation reproduces 
the shape of the data distribution very closely. 
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FIG. 18: (Color online) Acceptance of the full set of analysis 
cuts for both data and Monte Carlo events from a single ^^^Cf 
run in Phase II, as a function of kinetic energy. 



11. TRIGGER EFFICIENCY 

As discussed in Sec. IH the primary trigger for SNO 
was a coincidence of PMT hits within a 93 ns time win- 
dow, set to iVcoinc = 18 hits for the early part of Phase I 
and to A^coinc = 16 hits for the remainder of Phase I and 
all of Phase II. We define the 'efficiency' of the trigger 
as the probability that an event with A^coinc hits actually 
triggered the detector. Small shifts in the analog (DC- 
coupled) baseline, noise, and disabled trigger electron- 
ics channels could all lead to a non-unity efficiency. We 
measured the efficiency using the isotropic laser source, 
by triggering on the laser pulse and comparing an offline 
evaluation of the trigger (by counting hits in a sliding 
93 ns window) to the output of the hardware trigger. We 
found that for the A"coinc = 18 hit threshold, events with 
23 or more hits in coincidence triggered the detector with 
an efficiency greater than 99.9% and, for the A^coinc — 16 
hit threshold, the efficiency reached 99.9% at 21 hits. 
Figure [TO] shows the efficiency measured as a function of 
A'coinc, for Phase I at the higher A'coinc ~ 18 threshold, 
and for Phase II at the lower A'coinc = 16 hit threshold. 

For events at our T = 3.5 MeV analysis threshold, the 
mean number of hits in an event over the full 400 ns event 
window was ~ 30 for Phase I and ~ 27 for Phase II, 
with RMS's of 1.8 hits and 1.7 hits, respectively. The 
numbers of hits in the 400 ns event window and in the 



W0.8 



0.6 



0.4 



0.2 



• Phase I, Threshold at N . =18 

coinc 

o Phase II, Threshold at N =16 

come 



10 12 14 16 1^ 



20 22 24 26 28 30 
N .„ 



FIG. 19: Comparison of the trigger efficiencies in the two 
data-taking phases and for the two different thresholds used. 



93 ns trigger coincidence window differed primarily in the 
contribution from random PMT noise which, for both 
phases, contributed on average roughly 1 additional hit 
in the 400 ns event window. Thus, for both phases, the 
trigger efficiency was above 99.9% for all but a negligible 
fraction of events with a high enough A'coinc to pass the 
analysis cuts. 

Because our PDFs and overall normalization were de- 
rived from simulation, we compared the trigger-efficiency 
estimate from the data to the simulation's prediction. We 
also compared the idealized simulated trigger to a sim- 
ulation that included variations in the trigger baseline 
as measured by an online monitor. We found that the 
Monte Carlo simulation's prediction of trigger efficiency 
was in excellent agreement with our measurement for 
both SNO phases, and that the measured variations con- 
tributed a negligible additional uncertainty to our overall 
acceptance. 



12. UNCERTAINTIES ON THE NEUTRON 
CAPTURE EFFICIENCIES 



In Phase I, neutrons produced through the NC reaction 
and background processes were captured on deuterons 
within the heavy water, releasing a single 6.25 MeV 7 
ray. In Phase II, the neutrons were captured primarily 
on ■^^Cl, releasing a 7 cascade of total energy 8.6 MeV. 
The absolute cross sections for these capture reactions, 
along with detector acceptance, determined the rate of 
detected neutron events. The uncertainty on the neutron 
capture efficiency for Phase II overwhelmingly dominates 
that for Phase I in the final flux determinations because 
of the larger capture cross section. 

In this analysis, we used the Monte Carlo simulation 
to define the central values of the neutron capture effi- 
ciencies. Included in our simulation were the measured 
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isotopic purity of the heavy water, as well as its density 
and temperature and, for Phase II, the measured density 
of salt added to the D2O. 

To assess the systematic uncertainties on the neutron 
capture efficiencies, we used data taken with the ^^^Cf 
source deployed at many positions throughout the detec- 
tor, and compared the observed counting rates to simu- 
lations of the source runs. The differences between data 
and simulated events provide an estimate of the simu- 
lation's accuracy. The Phase I and Phase II data sets 
are noticeably different in their neutron detection effi- 
ciency because of the much larger capture cross section in 
Phase II, and the higher energy 7 cascade from neutron 
capture on chlorine. We therefore assessed the uncer- 
tainties in the two phases slightly differently, as discussed 
below. We also compared the results of this 'direct count- 
ing' approach with a 'time series analysis', in which the 
relative times of events were used to extract the capture 
efficiency. The two methods were in excellent agreement 
for both phases. Our capture efficiency uncertainty for 
Phase II is ±1.4%, and for Phase I it is ±2%. 



12.1. Phase II Neutron Capture Efficiency 
Uncertainties 

For the Phase II analysis, neutron events from the 
^^^Cf source were selected using the same burst algo- 
rithm that was used in previous SNO publications [8|. 
Neutrons were identified by looking for prompt fission 
7 events from the ^^^Cf decay, and tagging subsequent 
events that occurred within 40 ms. Figure [20] plots the 
neutron detection efficiency for each source run as a func- 
tion of radial position of the source in the detector, for 
both data and Monte Carlo simulated events. The source 
position for a run was determined by finding the mean 
reconstructed position of the prompt fission 7 events, to 
eliminate the large positioning uncertainties of the source 
deployment mechanism. The efficiencies shown in Fig. 1201 
were each fitted to a phenomenologically-motivated neu- 
tron detection efficiency function; 



;(s) = ^(tanh(B(s-C))-l), 



(12) 



where e(s) gives the neutron capture efficiency at source 
radius s. 

To determine the uncertainty on the simulation's pre- 
diction of capture efficiency, we first calculated the mean 
capture efficiency in the D2O volume, given the two func- 
tions shown in Fig. [20l as follows: 



Jo 



e{s)ds 



r 

Jo 



600.5 



(13) 



We took the difference of 0.8% between data and sim- 
ulation as a baseline uncertainty. (The mean detection 
efficiency measured this way was 35.6%). 

The normalization of the curves shown in Fig. 1201 de- 
pends on the strength of the ^^^Cf source, which we know 
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FIG. 20: (Color online) Data and Monte Carlo neutron detec- 
tion efRciencies in Phase II fitted to the phenomenologically- 
motivated neutron detection efficiency function. 

to 0.7% based on ex-situ measurements. An overall shift 
in reconstructed event positions, discussed in Sec.[Sl also 
changed the measured efficiency in data relative to the 
simulation results. By varying the value of this shift 
within its range of uncertainty we found it resulted in an 
additional 0.3% uncertainty in capture efficiency. The 
uncertainty in the fit parameters of the neutron detec- 
tion efficiency function was included conservatively by 
taking the entire statistical uncertainty on the data ef- 
ficiency measurements of Fig. 1201 which yields another 
0.9%. Lastly, we included a 0.1% uncertainty to account 
for the fraction of ^^°Cf in the ^^^Cf source (only ^'^^Cf 
is simulated by the Monte Carlo code). The overall un- 
certainty on the neutron capture efficiency, calculated by 
adding these in quadrature, was 1.4%. 

We checked these results by performing an indepen- 
dent time series analysis, in which we fit directly for the 
efficiency at each source deployment point based on the 
rates of neutron capture and 7 fission events (the source 
strength is not an input parameter). The fit included 
parameters associated with the overall fission rate, back- 
grounds from accidental coincidences, and the mean cap- 
ture time for neutrons. We obtained the efficiency as 
a function of source radial position, to which we fit the 
same efficiency function from Eq. [121 a-nd extracted the 
volume-weighted capture efficiency directly (rather than 
by comparison to Monte Carlo) . The mean efficiency cal- 
culated this way was 35.3 ± 0.6%, in excellent agreement 
with the value of 35.6% from the direct counting method, 
and well within the uncertainties on both measurements. 



12.2. Phase I Neutron Capture Efficiency 
Uncertainties 



The measurement of neutron capture efficiency uncer- 
tainty for Phase I is more difficult than for Phase II, pri- 
marily because the lower capture cross section in Phase I 
made identification of neutron events from the ^"'^^Cf 
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FIG. 21: (Color online) Comparison of the fit functions to the 
data and Monte Carlo in Phase I. 



source difficult. The number of detected neutrons per 
fission was small (less than one on average) , and the long 
capture time (roughly 50 ms) made coincidences more 
likely to be accidental pile-up of prompt fission 7s than 
neutrons following the 7s. 

Instead of using the burst algorithm, we separated neu- 
tron events from fission 7s based on their differing ener- 
gies and mean free paths in D2O. Events were required 
to be more than 150 cm from the source position and 
to have energies above the mean energy expected for a 
neutron capture event, for both data and Monte Carlo 
events. The detected rate of events after these cuts was 
used for the data and Monte Carlo simulation compari- 
son. 

An additional parameter was added to the neutron de- 
tection efficiency function for these data, as follows: 



e(s) = A(tanh(B(s - C)) - D), 



(14) 



and the resulting fits to data and Monte Carlo are shown 
in Figure [211 

The difference of the volume- weighted integrals of the 
two curves is just 0.9%, but the small value is clearly 
due to cancellation differences at different radii. The 
shape difference is driven by small differences between 
the data and Monte Carlo fits at large radii, which are 
likely due to unassessed systematic errors on the data 
points themselves. We included additional uncertainties 
to account for these. In particular, we included a 0.6% 
uncertainty associated with the statistical uncertainties 
of the data and Monte Carlo neutron detection efficiency 
function parameters, and an additional 0.6% uncertainty 
associated with knowledge of the source position. We 
also included a further uncertainty of 0.9% to account 
for data and Monte Carlo differences in the energy cut 
applied to select neutrons. 

We applied the same source-strength uncertainties as 
for the Phase II analysis, namely the 0.7% absolute 
source strength calibration, and 0.1% from the (unmod- 
eled) contamination of ^^''Cf in the ^^^Cf source. The 
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FIG. 22: (Color online) Mean neutron capture time from the 
time series analysis in Phase I as a function of source position. 
The Une shows the best fit to the simulation using a cubic 
polynomial. 

total uncertainty on the neutron capture efficiency for 
Phase I comes to 2%. 

To check our estimates, we also performed a time series 
analysis of the ^^^Cf data. Unlike Phase II, for Phase I 
we cannot extract the absolute efficiency to compare with 
that derived from the direct counting method because of 
the 150 cm reconstruction cut. Instead, we performed 
the time series analysis on both Monte Carlo and source 
data runs, and compared them. We found the fractional 
difference between the source-derived and Monte Carlo- 
derived efficiencies to be just 0.3%, well within the 2% 
uncertainty obtained from the direct counting method. 
One output of the time series analysis is the neutron 
capture time: the time between neutron emission from 
the ^^^Cf source and capture on a deuteron. Figure 
shows the neutron capture time as a function of source 
radial position for both data and Monte Carlo. As the 
^^^Cf source approaches the acrylic vessel and light wa- 
ter region, the capture time decreases significantly. The 
overall agreement between the measured capture times 
in data and Monte Carlo is very good throughout most 
of the volume. 



13. BACKGROUNDS 

Lowering the energy threshold opened the analysis 
window to additional background contamination, pre- 
dominantly from radioactive decays of ^-'^^Bi and ^os-pj 
in the 238^ g^j^jj 232rpj^ chains, respectively. In Phase II, 
neutron capture on ^^^Na produced a low level of ^"'Na 
in the detector which, in its decay to ^^Mg, produced a 
low energy (3 and two 7s. One of these 7s has an energy 
of 2.75 MeV, which could photodisintegrate a deuteron. 
The result was some additional electron-like and neutron 
capture background events. In addition, radon progeny 
that accumulated on the surface of the AV during con- 
struction could have created neutrons through (a,n) reac- 
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tions on isotopes of carbon and oxygen within the acryhc. 

In the past, most of these backgrounds were esti- 
mated using separate self-contained analyses and then 
subtracted from the measured neutrino fluxes. In this 
analysis, the Monte Carlo simulation was used to create 
PDFs for each of 17 sources of background events (except 
for PMT /3-7 events, for which an analytic PDF was used 
in each phase, as described in Sec. 113. 31) . and the numbers 
of events of each type were parameters in the signal ex- 
traction fits. Table 1x1 lists the sources of physics- related 
backgrounds that were included in the fits. 



Detector Region 


Phase I 


Phase II 


D2O volume 


Internal ^"Bi 


Internal ^^*Bi 




Internal 


Internal ^°^Tl 






^^Na 


Acrylic vessel 


Bulk 2"Bi 


Bulk ^i^Bi 




Bulk ^''^Tl 


Bulk ^''^Tl 




Surface (a,n) ns 


Surface {a,n) ns 


H2O volume 


External ^"Bi 


External ^^^Bi 




External ^"^TI 


External ^°^Tl 




PMT ^-7S 


PMT P-ys 



TABLE X: The sources of physics-related background events 
in the LETA analysis. 



All of the Monte Carlo-generated PDFs were verified 
using calibration sources. Ex-situ measurements [13, HI] 
of background levels in the D2O and II2O provided a 
priori information for several of them, which were used 
as constraints in the signal extraction fits. In addition, 
corrections were applied after the signal extraction fits 
to account for a number of background event types that 
contributed much smaller levels of contamination. The 
following sections describe these procedures. 

13.1. Background PDFs 

Most of the PDFs used in the signal extraction were 
created from Monte Carlo simulations of the specific 
event types. However, because of the limited number of 
simulated PMT /3-7 events available in the radial range 
of interest, an analytic parameterization of the PDF was 
used, as described in Sec. 113.31 This was verified by com- 
parison to the simulation and uncertainties associated 
with the value of each parameter were propagated in the 
signal extraction fits. 

The remainder of the background PDFs were verified 
by comparison of calibration data to simulated events. 
The D2O and H2O backgrounds were verified using the 
D2O- and H20-region radon spikes in Phase II and cali- 
bration sources deployed in these regions. Bulk AV back- 
grounds were verified using the and sources, 
and surface (a,n) neutrons using the ^^^Cf source de- 
ployed near the AV. 



In all cases, the data and Monte Carlo event distri- 
butions agreed to within the systematic uncertainties al- 
ready defined for the PDFs. Figure [23] shows the energy 
dimension of a fit to the internal radon spike. The fit 
was performed using the unbinned signal extraction code 
(see Sec. 114. Sp in a simplified configuration, as described 
in Sec. 18.3.11 The result is a good fit to the data, in par- 
ticular at low energy. Figure [M] shows a comparison of 
data to simulation for the ^"^^Th source deployed near the 
AV. A band is shown for the simulated events, represent- 
ing the quadrature sum of the statistical uncertainties 
with the effect of applying the dominant systematic un- 
certainties. The distributions in Tes, and /3i4 show 
good agreement within the la uncertainties. 
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FIG. 23: (Color online) One dimensional projection of the fit 
to the internal radon spike data. 

The cross section for photodisintegration affects the 
relative normalization of the neutron and electron parts 
of the background PDFs. The simulation used a theo- 
retical value for the cross section and the associated 2% 
uncertainty was propagated in the signal extraction fits. 

The simulation of ^''Na events used to generate a PDF 
was done under the assumption of a uniform distribution 
of events within the detector, since a primary source of 
^"^Na was the capture of neutrons produced by deployed 
calibration sources on ^■^Na. ^'*Na was also introduced 
via the neck, and via the water systems, which connected 
near the top and bottom of the AV. Therefore, the signal 
extraction fits were redone with different spatial distri- 
butions, in which the events originated either at the neck 
of the AV or at the bottom, with a conservatively cho- 
sen 10% linear gradient along the z-axis. The difference 
from the baseline (uniform distribution) fit was taken as 
a systematic uncertainty. 

13.2. Low Energy Background Constraints 

Several radioassays were performed during data taking 
to measure the concentrations of radon and radium in the 



30 



c 

I Ik' 



§10-1 



10- 





• Data 




MC 




(a) 


— •— 




" 


— •— 

, < 1 < < , < 1 < , < < 1 < , < , 



Teff(MeV) 



o 



- — ^ Data 












L (b) 


■ 

• 


/ 

— •— 

-•— 


, , , 1 , , , , 1 


. . 1 . . . . 1 . . . . 1 . . . 


1 , , . . 1 . , , 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 



2.5: 



— ^ Data 

(c) 



« 1 5- 

B - 

i-i 
O 

^ 1- 



0.5- 



background events in the analysis window. During Phase 
II, there was a leak in the assay system used to measure 
the 238y chain contamination that was not discovered 
until after data taking had ended, so there is no accu- 
rate constraint on the ^38^ level in the D2O during that 
phase. Other limits based on secondary assay techniques 
were found to be too loose to have any impact on the 
signal extraction results and so were disregarded. The 
results of the assays are given in tables IXII and IXIIl 



Phase Isotope Concentration (x 10 g/g of D2O) 


I 238U 


10 1 +^'' 
^"■^ -2.0 


232 rpj^ 


2.09±0.21(stat) io;9i(syst) 


II 238u 




232rp-|_^ 


1.76±0.44(stat) tagl^syst) 



TABLE XI: ^^*U and ^""^Th concentrations in the D2O vol- 
ume, determined from ex-situ radioassays in Phases I and II. 



Phase Isotope Concentration (g/g of II2O) 


I 238u 


29.5 ± 5.1x IQ-" 




s.itilx 10-^^ 


II 238u 


20.6 ± 5. Ox 10"^* 


232 rpj^ 


5.2± 1.6x 10"^" 



TABLE XII: ^^^'U and ^^^Th concentrations in the H2O vol- 
ume, determined from ex-situ radioassays in Phases I and II. 



These concentrations were converted into an expected 
number of events and were applied as constraints in the 
signal extraction fits, as described in Sec. 114.51 

In-situ analyses (36| were used to predict the number 
of background events from ^^Na decays in Phase II. The 
predicted value of 392 ± 117.6 events was applied as a 
constraint in the signal extraction fits. 
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FIG. 24: (Color online) Comparison of data to simulation 
for ^•^^Th source runs near the AV in Phase II, in (a) Tcfi, 
(b) 7?^, and (c) /3i4. The band represents the la uncertainty 
on the Monte Carlo-prediction, taking the quadrature sum 
of the statistical uncertainties with the effect of applying the 
dominant systematic uncertainties. 



D2O and II2O regions, as described in previous publica- 
tions .34 . .35j . Although equilibrium was broken in the 
decay chains, the results are expressed in terms of equiva- 
lent amounts of ^^^U and ^^^Th assuming equilibrium for 
ease of comparison with other measurements. The results 
were used to place constraints on the expected number of 



13.3. PMT /3-7 PDF 

We use the term "PMT events" to refer to all radioac- 
tive decays in the spherical shell region encompassing 
the PMTs and the PSUP. These events were primarily 
208rpj decays originating from ^^^^Th contamination in the 
PMT/PSUP components. 

PMT events occurred at a high rate, but only a tiny 
fraction of them reconstructed inside the signal box and 
within the fiducial volume: in Phase I, the acceptance 
was only 1.7 x 10^^ and in Phase II it was 5.9 x 10^^. 
Therefore, an enormous amount of computer time would 
be needed to generate enough events to create a PDF. 
Creation of a multi-dimensional PDF based entirely on 
simulation was therefore deemed to be impractical. 

A high rate thorium source was deployed near the 
PSUP in both phases to help model these events. How- 
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ever, interpretation of this data was complicated by the 
fact that a point source with a sufficiently high rate tends 
to produce significant 'pile-up' of multiple events that 
trigger in the same time window. This pile-up changes 
the topology of the events to the extent that they are 
not characteristic of PMT /3-7S, so they cannot be used 
directly as a model. 

Therefore, an analytic parameterization of the PDF, 
given in Eq. ([15]), was used. For this, the cos di- 
mension was assumed to be flat; the remaining three- 
dimensional PDF was of the form: 

PpMTiTeS, Pli ,R^)= e-^^"" X (e^«' + C) 

X M{pu\i3u = D + ER^,a = F)il5) 

where Af{x\x,a) is a Gaussian distribution in x with 
mean x and standard deviation a. The /3i4 dimension was 
determined from a Gaussian fit to Monte Carlo events, 
in which /3i4 was allowed a linear dependence on R^. 

The source location of the PMT events, their large 
number, and the fact that they must reconstruct nearly 
3 m from their origin to appear inside the fiducial volume 
means that they have features that distinguish them from 
other sources of backgrounds. Therefore, we were able to 
extract a prediction for the total number of PMT events, 
as well as for the shape of the energy and radial dimen- 
sions of the PDF, from the data itself, by performing a 
bifurcated analysis. 

In a bifurcated analysis, two independent cuts are se- 
lected that discriminate signal from background. The 
behavior of these cuts when applied both separately and 
in combination is used to assess the number of signal 
and background events in the analysis window. We as- 
sume that the data set consists of ly signal events and (3 
background events, so that the total number of events is 
S = P + V. The background contamination in the final 
signal sample is just the fraction of f3 that passes both 
cuts. If the acceptances for background and signal events 
by cut i are yi and Xi, respectively, the contamination is 
2/12/2/? and the number of signal events is xiX2V. 

Given the number, a, of events that pass both cuts, 
the number, b, that fail cut 1 but pass cut 2, and the 
number, c, that pass cut 1 but fail cut 2, we then relate 
these with a system of equations: 

a + c = xiv + yiP, (16) 

a + b = x2v + y2P, (17) 

a = X1X2V + yiy2P, (18) 

P + v = S, (19) 

which we solve analytically, using Monte Carlo- 
predictions for the cut acceptances, to determine the con- 
tamination, K = yij/2/5, in the signal sample. A feature 
of this method is that it produces a contamination esti- 
mate without including events from the signal box (those 
that pass both cuts) in the analysis. 

In this analysis, the 'background' comprised the PMT 
events and the 'signal' all other events, including both 



neutrino interactions and non-PMT radioactive decays. 
The cuts chosen were the in-time ratio (ITR) cut, because 
it selected events that were reconstructed far from their 
true origin, and the early charge (EQ) cut because it 
selected events in which a large amount of light produced 
hits early in time in a small number of tubes. These tend 
to be characteristics of PMT events (see Sec. I10.2"|l . 

For a bifurcated analysis to work, the probabilities of 
passing the cuts must be statistically independent. To 
demonstrate this, we loosened the cuts, and found that 
the increase in the number of background events agreed 
well with what would be expected if they were indepen- 
dent. 

One result of the bifurcated analysis is a prediction for 
the number of PMT events in the analysis window, which 
was used as a constraint in the binned likelihood signal 
extraction fits, as described in Sec. 114.51 

The acceptance of signal events {X1X2) ^ 1-0 and there- 
fore some non-PMT events were also removed by the 
cuts. Such events falsely increase the count of back- 
ground events in the three 'background boxes'. We lim- 
ited the impact of this effect by restricting the analysis to 
the 3.5-4.5 MeV region, which was overwhelmingly dom- 
inated by PMT events. We also included a correction 
for the number of non-PMT events in each of the back- 
ground boxes by using estimates from the Monte Carlo 
simulation for the acceptance of all other signals and 
backgrounds, and verifying these predictions with radon 
spike data. (^^*Bi, a radon daughter, is the dominant 
background other than the PMT events in this region) . 

To estimate the number of non-PMT events in each 
of the three background boxes, we multiplied the Monte 
Carlo-predicted acceptances of non-PMT events by the 
expected total number of these events in the data set. 
The procedure was therefore iterative: a PMT PDF was 
created using initial estimates for the total number of 
non-PMT events in the data set and their acceptances; 
the bifurcated analysis was used to predict the number 
of PMT events in the signal box; the data were re-fit 
with this new PMT constraint; the total number of non- 
PMT events in the data set, based upon the new fit, was 
then used to update the non-PMT event correction in 
the background boxes in the bifurcated analysis, and so 
on. In practice, the bifurcated analysis itself was simply 
included within the signal extraction fit, so the prediction 
for the number of PMT events could be recalculated as 
the fit progressed, and the penalty factor in the likelihood 
calculation from the resulting constraint could be varied 
accordingly. To determine systematic uncertainties on 
this overall procedure, we tested the analysis on sets of 
fake data and compared the prediction of the bifurcated 
analysis to the known true number of PMT events 
in the signal box. 

We verified the bifurcated analysis results by compar- 
ing the prediction of the total number of PMT P-j events 
in the signal box to an estimate made with an indepen- 
dent analysis performed outside the fiducial volume. This 
independent analysis looked for events that occurred at 
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high radius and were inward-pointing, which are charac- 
teristics of PMT /3-7 events, and extrapolated that count 
into the fiducial volume. The measurements agreed with 
the bifurcated analysis to well within the uncertainties 
on the two methods. 

To predict the shape of the PMT PDF, the bifurcated 
analysis was performed in discrete bins in Tcff and R^. 
Unlike the prediction for the total number of PMT events 
in the data set, this calculation was not included in the 
signal extraction, so a fixed estimate of the contamination 
of non-PMT events in the three background boxes was 
applied. This estimate was derived from a signal extrac- 
tion fit performed on a small subset of the data. To take 
uncertainties into account, bifurcated analyses were per- 
formed on Monte Carlo-generated 'fake' data sets with 
the dominant systematic and statistical uncertainties ap- 
plied in turn, to determine the effect of each on the ex- 
tracted shape for the PMT PDF. The differences of the 
results from the unshifted version were added in quadra- 
ture to obtain an additional uncertainty on the shape. 

A number of functional forms were fit to the T^g and 
distributions to determine the best parameterizations 
for the shapes. An exponential was found to be a good fit 
to the energy profile and an exponential plus a constant 
offset to the radial distribution (see Eq. ([15])). The fit 
results for Phase II are shown in Figure [511 

The parameters from the fits shown in Fig. [55] were 
varied in the signal extraction by applying a Gaussian 
penalty factor to the likelihood function, as described in 
Sec. 114.51 The mean of the Gaussian was the central fit 
value from Fig. [55] and the standard deviation was taken 
as the total uncertainty in this value, including both the 
fit uncertainty from Fig.[25]and the additional systematic 
uncertainties described above. Results for both phases 
are shown in Table IXllIl The fits to the bifurcated analy- 
sis prediction for the distribution showed a significant 
correlation between the exponent and the offset, with cor- 
relation coefficients of 0.846 and 0.883 in Phases I and II, 
respectively. This correlation was included in the Gaus- 
sian penalty factor in the signal extraction fits. 



Parameter 


Phase 


I 


Phase 


II 


Energy exponent, A (/MeV) 


-5.94 


± 


0.96 


-6.37 ± 


0.81 


exponent, B 


5.83 


± 


0.96 


5.28 ± 


0.79 


offset, C 


-0.40 


± 


1.43 


-0.32 ± 


1.16 



TABLE XIII: Parameters defining the PMT PDF shape, as 
defined in Eq. (fT5|l . 



13.4. Limits on Instrumental Backgrounds 

Because instrumental background events were not 
modeled by the simulation, their contamination in the 
analysis window was determined directly from the data. 
A bifurcated analysis was used, similar to that described 




FIG. 25: (Color online) Predicted shapes for the PMT PDF 
in (a) i?^ and (b) Tcs in Phase II. 



in Sec. 113.31 In this instance, two sets of cuts were used 
to define the analysis: the instrumental cuts and the 
high-level cuts, described in Sec. [TU] The numbers of 
events in the data set failing each and both sets of cuts 
were used to estimate the contamination by instrumental 
backgrounds. 

As was done in Sec. 113. 3[ a prediction of the number 
of good (physics) events that failed the instrumental cuts 
was used to correct the number of events in each of the 
background boxes. We obtained this prediction using the 
cut acceptances given in Sec. 110.5.11 and an estimate of 
the numbers of signal and radioactive background events 
in the data set. The analysis was performed at two energy 
thresholds in order to study the energy dependence of the 
contamination. Results are given in Table IXIVl 

Since these events were not modeled in the simulation, 
it is difficult to directly predict their effect on the signal 
extraction fit results. However, because virtually all of 
them fall into the lowest energy bin, they are unlikely 
to appear like neutron events. Since the Toff distribu- 
tions of CC and ES signals were unconstrained in the 
signal extraction fit, they could mimic these event types. 
Therefore, a conservative approach was taken, in which 
the estimated contamination from the 3.5 McV analysis 



33 





Threshold 


Phase 


3.5 MeV 4.0 MeV 


I 


2.64 ± 0.22 0.09 ± 0.42 


II 


4.48 ± 0.27 0.52 ± 0.23 



TABLE XIV: Estimated number of instrumental contamina- 
tion events in the full data set at different analysis thresholds. 



was applied as an additional uncertainty in the lowest 
energy bin for both the CC and ES signals. 

13.5. Atmospheric Backgrounds 

The NUANCE neutrino Monte Carlo simulation pack- 
age p??! was used to determine the contribution of at- 
mospheric neutrino events to the data set. The esti- 
mated number of atmospheric neutrino events was not 
large enough to merit introducing an additional event 
type into the already complex signal extraction proce- 
dure. Instead, 15 artificial data sets were created that 
closely represented the best estimate for the real data 
set, including all neutrino signals and radioactive back- 
grounds in their expected proportions. The NUANCE 
simulation was used to predict the distribution of atmo- 
spheric neutrino events in each of the four observable 
parameters used to distinguish events in the signal ex- 
traction fit (see Sec. [HI), and a number of such events 
were included in each artificial data set, drawn from the 
estimate for the number in the true data. Signal ex- 
traction was performed on these sets to determine which 
signal the events would mimic in the extraction. This re- 
sulted in a small correction to the NC flux of 4.66 ± 0.76 
and 17.27 ± 2.83 events to be subtracted in Phases I and 
II, respectively, and small additional uncertainties for the 
CC and ES rates, mostly at the sub-percent level. 

Atmospheric events were often characterized by a high- 
energy primary followed by several neutrons. Therefore, 
there was significant overlap with events identified by the 
'coincidence cut', which removed events that occurred 
within a fixed time period of each other. This overlap was 
exploited to verify the predicted number of atmospheric 
events. Without application of the coincidence cut, a to- 
tal of 28.2 ± 5.4 and 83.9 ± 15.9 atmospheric neutrino 
events were predicted in Phases I and II, respectively. 
The coincidence cut reduced these numbers to 21.3 ±4.0 
and 29.8±5.7 events, which were the numbers used in the 
creation of the initial artificial data sets. A second group 
of sets was created, using the pre-coincidence cut esti- 
mates for the number of events, to determine the change 
in the NC flux due to the additional events. The signal 
extraction was then performed on a subset of the real 
data, both with and without the application of the coin- 
cidence cut, and the observed difference in the NC flux 
was entirely consistent with the predictions, thus verify- 
ing the method used to derive the NC flux correction. 



13.6. Isotropic Acrylic Vessel Background (lAVB) 

Early in the SNO analyses, a type of instrumental 
background was discovered that reconstructed near the 
AV and was characterized by very isotropic events (/3i4 < 
0.15). At higher energies (A^hit > 60), these events form 
a distinct peak in a histogram of /3i4, and they are easily 
removed from the data by a combination of the fiducial 
volume and isotropy cuts. However, at lower energies, 
position reconstruction errors increase and the isotropy 
distributions of the lAVB and other events broaden and 
join, so that removal of the lAVB events by these cuts is 
no longer assured. 

Accurate simulation of these events is difficult because 
the physical mechanism that produces the lAVB events 
has not been identified and crucial lAVB event charac- 
teristics cannot be predicted. These include the light 
spectrum, photon timing distribution, location, and ef- 
fective event energy. To circumvent this problem, sim- 
ulated events were generated that covered a wide range 
of possibilities. Three event locations were modeled: on 
the exterior and interior AV surfaces, and uniformly dis- 
tributed within the AV acrylic. Events were generated at 
three different photon wavelengths that cover the range 
of SNO detector sensitivity: 335, 400, and 500 nm. The 
photons were generated isotropically, with the number 
of photons in an event chosen from a uniform distribu- 
tion with a maximum above the energy range used in 
the neutrino analysis. The photon time distribution was 
a negative exponential, with the time constant for an 
event chosen from a truncated Gaussian with mean and 
standard deviation of 5 ns. 

Using PDFs derived from the simulated event sam- 
ples, maximum likelihood signal extraction code was used 
to estimate the number of lAVB events in the data in 
the vicinity of the AV, between 570 and 630 cm from 
the detector center, in accompaniment with the CC, ES, 
and NC neutrino event types and ^oS'pj g^^j 2i4gj back- 
grounds in the D2O, AV, H2O, and PMTs. This was done 
separately for each of the nine simulated photon wave- 
length/event location combinations. Because the energy 
distribution of the lAVB events was unknown, the lAVB 
extractions were done as a function of Nhit in 11 bins. 
The ratio of the number of lAVB events that passed all 
the neutrino cuts to those that fit near the AV in each 
Nhit bin was calculated for each simulated lAVB case as a 
function of event energy. These ratios were used, together 
with the estimated numbers of such events near the AV, 
to estimate the lAVB contamination in the neutrino sam- 
ple as a function of energy for each of the simulated lAVB 
cases. 

The polar-angle distributions of hit PMTs in the sim- 
ulated lAVB events were studied in a coordinate system 
centered on the middle of the AV, with its z-axis along 
the radial vector through the fitted event location. There 
are marked differences in these distributions among the 
different simulated cases due to optical effects of the AV. 
Comparisons of these distributions were made between 
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simulated events and high Nhit, high isotropy events in 
the data that reconstruct near the AV (presumed to be 
lAVB events). A fit was made to find the weighted com- 
bination of the simulated cases that best fit the high Nhit 
data. The resulting weights were assumed to be valid at 
all energies, and were used together with the contamina- 
tion ratios discussed above: first, to estimate the total 
lAVB background expected in the neutrino analysis data 
set as a function of energy (totaling 27 and 32 events 
above 3.5 MeV in Phases I and II, respectively) and, sec- 
ond, to generate a set of simulated lAVB events repre- 
sentative of those expected to contaminate the neutrino 
data. 

A test similar to that described in Sec. 113.51 was per- 
formed. Fifteen artificial data sets were created that also 
contained simulated lAVB events based on estimates of 
the weighted contributions of the simulated cases and 
their energy distributions. It was found that the major- 
ity of the lAVB events fit out as other background event 
types, so that the result of adding the simulated lAVB 
background was only small additional uncertainties for 
each of the neutrino flux parameters, with no required 
corrections. The increase in uncertainty for the NC fiux 
was evaluated at 0.26%. The increases of the CC uncer- 
tainties were also mostly at the sub-percent level, and the 
increase in uncertainties on the ES rates were so small as 
to be negligible (< 0.01%). 

13.7. Additional Neutron Backgrounds 

A full study of other possible sources of neutron back- 
ground events, such as from events such as {a,n) re- 
actions and terrestrial and reactor antineutrino interac- 
tions, was presented in previous publications 0, The 
full set of simulated NC events was used to adjust these 
numbers for the lowered energy threshold and for the live 
times and detection efficiencies in the two phases to give 
a final correction to the NC fiux of 3.2±0.8 and 12.0±3.1 
neutron capture events in Phases I and II, respectively. 

14. SIGNAL EXTRACTION METHODS 

An extended maximum likelihood method was used to 
separate event types based on four observable parame- 
ters: the effective electron kinetic energy, Tes; the angle 
of the event direction with respect to the vector from 
the Sun, cos Oq ; the normalized cube of the radial posi- 
tion in the detector, R^; and the isotropy of the PMT 
hits, /3i4. Two independent techniques were used, as 
described in sections 114.21 and 114.31 One method used 
binned PDFs and the other an unbinned, "kernel estima- 
tion" approach. 

We performed two distinct types of fit. The first ex- 
tracted the detected electron energy spectra for CC and 
ES events in individual Toff bins, without any model con- 
straints on the shape of the underlying neutrino spec- 



trum. We refer to this as an 'unconstrained' fit. The 
second fit exploited the unique capabilities of the SNO 
detector to directly extract the energy-dependent i^e sur- 
vival probability (Sec. 114. 4p . The survival probability 
was parameterized as a polynomial function and applied 
as a distortion to the ^B neutrino energy spectrum (taken 
from HI). The shapes of the CC and ES T^s spec- 
tra were recomputed from the distorted ^B spectrum as 
the fit progressed, allowing the polynomial parameters to 
vary in the fit. The overall fluxes were also constrained 
in this fit through the requirement of unitarity. The fea- 
tures in common for the two signal extraction approaches 
are described below. 

The types of events included in the fit were the three 
neutrino interaction types (CC, ES and NC) and 17 back- 
ground event types across the two phases of data, as de- 
fined in Table |Xl The likelihood was maximized with 
respect to the number of events of each signal type, and 
several systematic parameters affecting the shapes of the 
PDFs, as described in Sections [Ti?^ and [Ti31 

To extract energy spectra for the CC and ES neutrino 
signals in the unconstrained fits, CC and ES PDFs were 
created in discrete T^g intervals and the fitted numbers 
of events in these intervals were allowed to vary inde- 
pendently. The energy spectra for events from the NC 
interaction and from radioactive backgrounds have no de- 
pendence on the neutrino oscillation model, and so the 
shapes of these spectra were fixed within their systematic 
uncertainties. 

The flux of solar neutrinos was assumed to be constant, 
so a single set of neutrino-related fit parameters was ap- 
plied to both phases. Therefore, the neutrino signal pa- 
rameters varied in the fit were an NC rate and a number 
of CC and ES rates in discrete energy intervals, as defined 
in Sections 114.21 and 114.31 Although SNO was primarily 
sensitive to the ^B chain of solar neutrinos, we included a 
fixed contribution of solar hep neutrinos, which was not 
varied in the fit. Based on results from a previous SNO 
analysis we used 0.35, 0.47, and 1.0 times the Solar 
Standard Model (SSM) prediction for CC, ES, and NC 
hep neutrinos, respectively. Taken together, these corre- 
spond to 16.4 events in Phase I and 33.3 events in Phase 
II. 

To take into account correlations between parame- 
ters, multi-dimensional PDFs were used for all signals. 
In the unconstrained fits, CC and ES were already di- 
vided into discrete energy bins, and three-dimensional 
PDFs were created in each bin for the other observ- 
ables: P(/3i4,i?^ , cosOq). In the survival probability fits, 
fully four-dimensional PDFs were used for CC and ES 
events. For the NC and background PDFs the cos 6*0 
distribution is expected to be fiat, since there should 
be no dependence of event direction on the Sun's po- 
sition, but correlations exist between the other observ- 
ables. For these event types, the PDFs were factorized 

as P(reff,/3i4,i?3) X P(COS00). 

Uncertainties in the distributions of the observables 
were treated as parameterized distortions of the Monte 
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Carlo PDF shapes. The dominant systematic uncertain- 
ties were allowed to vary in the fit in both signal extrac- 
tion methods. Less significant systematics were treated 
as in previous SNO analyses 0], using a 'shift-and-refit' 
approach: the data were refit twice for each systematic 
uncertainty, with the model PDFs perturbed by the es- 
timated positive and negative 1 u values for the uncer- 
tainty in a given parameter. The differences between the 
nominal flux values and those obtained with the shifted 
PDFs were taken to represent the 68% C.L. uncertainties, 
and the individual systematic uncertainties were then 
combined in quadrature to obtain total uncertainties for 
the fluxes. 



14.1. Systematic Uncertainties: Phase Correlations 

Uncertainties related to theoretical quantities that are 
unaffected by detector conditions (such as the photo- 
disintegration cross section uncertainty) were applied to 
both phases equally. Uncertainties in quantities depen- 
dent on detector conditions (such as energy resolution) 
were treated independently in each phase. Uncertain- 
ties in quantities that partly depend on the operational 
phase (such as neutron capture efficiency, which depends 
both on a common knowledge of the ■^^^Cf source strength 
and on the current detector conditions) were treated as 
partially correlated. For the latter, the overall uncer- 
tainty associated with each phase thus involved a com- 
mon contribution in addition to a phase-specific uncer- 
tainty. Since neutron capture events were more similar 
to electron-like events in Phase I than in Phase II, several 
of the neutron-related uncertainties applied to Phase II 
only. The correlations are summarized in Table IX VI 



14.2. Binned-Histogram Unconstrained Fit 

In this approach, the PDFs were created as three- 
dimensional histograms binned in each observable dimen- 
sion, as summarized in Table IXVII For CC and ES, 
three-dimensional PDFs were created in each Tog inter- 
val, to fully account for correlations between all four ob- 
servable dimensions. Fifty rate parameters were fitted: 
the CC and ES rates in each of 16 spectral bins, the 
NC normalization and 17 background PDF normaliza- 
tions. Dominant systematic uncertainties were allowed 
to vary within their uncertainties, or 'fioated', by per- 
forming one-dimensional scans of the likelihood in the 
value of each systematic parameter. This involved per- 
forming the fit multiple times at defined intervals in each 
systematic parameter and extracting the value of the like- 
lihood, which included a Gaussian factor whose width 
was defined by the independently estimated uncertainty 
on that parameter, as described in Sec. 114.51 This com- 
bined a 'priori knowledge from the calibration data and 
Monte Carlo studies used to parameterize systematic un- 
certainties with information inherent in the data itself. If 



Systematic uncertainty 


Correlation 


Energy scale 


Both 


Electron energy resolution 


Uncorrelated 


Neutron energy resolution 


Phase II only 


Energy linearity 


Correlated 


/3i4 electron scale 


Correlated 


/3i4 neutron scale 


Phase II only 


/3i4 electron width 


Correlated 


/3i4 neutron width 


Phase II only 


/?i4 energy dependence 


Correlated 


Axial scaling 


Uncorrelated 


z scaling 


Uncorrelated 


X, y, z offsets 


Uncorrelated 


X, y, z resolutions 


Uncorrelated 


Energy dependent fiducial volume Uncorrelated 


cos 00 resolution 


Uncorrelated 


PMT Tcff exponent 


Uncorrelated 


PMT exponent 


Uncorrelated 


PMT R^ offset 


Uncorrelated 


PMT /3i4 intercept 


Uncorrelated 


PMT /?i4 radial slope 


Uncorrelated 


PMT /?i4 width 


Uncorrelated 


Neutron capture 


Both 


Photodisintegration 


Correlated 


^*Na distribution 


Phase II only 


Sacrifice 


Uncorrelated 


lAVB 


Uncorrelated 


Atmospherics backgrounds 


Uncorrelated 


Instrumental contamination 


Uncorrelated 


Other neutrons 


Uncorrelated 



TABLE XV: Phase correlations of the systematic uncertain- 
ties. "Correlated" refers to a correlation coefficient of 1.0 
between the phases and "uncorrelated" refers to a coefficient 
of 0.0. "Both" means an uncertainty was treated as partially 
correlated between the phases. 



a new likelihood maximum was found at an offset from 
the existing best estimate of a particular systematic pa- 
rameter, then the offset point was defined as the new 
best estimate. An iterative procedure was used to take 
into account possible correlations between parameters. 
The final uncertainties on each parameter were defined 
by where the log likelihood was 0.5 less than at the best- 
fit point, and the differences in each fitted flux parameter 
between these points and the best-fit point were taken as 
the associated systematic uncertainties for that parame- 
ter. For more details of this approach, see f39| . 

The parameters floated using this approach, along with 
their relevant correlations, as described in Sec. 114.11 were: 

• Energy scale (both correlated and uncorrelated in 
each phase) 
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Observable 


Min Max 


Bins 


Bin width 


CC, ES Tcff 


3.5 MeV 11.5 MeV 


16 


0.5 MeV 


Other Tcff 


3.5 MeV 5.0 MeV 


6 


0.25 MeV 








5.0 MeV 11.5 MeV 


13 


0.5 MeV 


cos 6q 


-1.0 1.0 


8 


0.25 




0.0 0.77025 


5 


0.15405 


014 


-0.12 0.95 


15 


0.0713 



TABLE XVI: PDF configurations used for the binned- 
histogram signal extraction approach. 



• Energy resolution (uncorrelated in each phase) 

• /3i4 scale for electron-like events (correlated be- 
tween phases) 

• PMT /3-7 exponent (uncorrelated in each phase, 
see Sec. 113. 3|) 

• PMT f3-j offset (uncorrelated in each phase, see 
Sec. MM 

• PMT (3-j Tcs exponent (uncorrelated in each 
phase, see Sec. 113. 3|) 

The remaining systematic uncertainties were applied us- 
ing the 'shift-and-refit' approach. 

14.3. Unbinned Unconstrained Fit Using Kernel 
Estimation 

In this approach, the PDFs were created by kernel esti- 
mation. Like standard histogramming techniques, kernel 
estimation starts with a sample of event values, drawn 
from an unknown distribution, P{x). Based on this finite 
sample, the parent distribution is approximated by P{x), 
which is a sum of kernel functions, Ki{x), each centered 
at an event value from the sample: 

1 " 

Pix)^-y^Mx-U). (20) 

i=l 

The most common choice of form of kernel functions is 
the normalized Gaussian distribution, 

K{x/h) = ^e-(-/'')V2, (21) 

where h is called the bandwidth of the kernel. One can 
pick a different bandwidth, hi, for the kernel centered 
over each event. 

Kernel-estimated density functions have many useful 
properties. If the kernel functions are continuous, then 
the density function will also be continuous. In one di- 
mension, kernel estimation can also be shown to con- 
verge to the true distribution slightly more quickly than 



a histogram with bin size the same as the kernel band- 
width. Generalizing the kernel estimation method to 
multiple dimensions is done by selecting a kernel with 
the same dimensionality as the PDF. We used a multi- 
dimensional Gaussian kernel that was simply the product 
of one-dimensional Gaussians. We followed the prescrip- 
tion given in [i^] for the selection of bandwidths for each 
event in each dimension. 

By varying the values associated with the events in the 
PDF sample individually, kernel estimation can very nat- 
urally be extended to incorporate systematic variation 
of PDF shapes. For example, energy scale is incorpo- 
rated by a transformation of the simulated event values, 
ti ^ {1 + a) X ti, where a is a continuously variable pa- 
rameter. Such transformations preserve the continuity 
and analyticity of the PDF. We can then add these sys- 
tematic distortion parameters to the likelihood function, 
and also optimize with respect to them using a gradi- 
ent descent method. This allows correlations between 
systematics and neutrino signal parameters, as well as 
between systematics themselves, to be naturally handled 
by the optimization algorithm. In addition, the infor- 
mation in the neutrino data set itself helps to improve 
knowledge of detector systematics. 

Three kinds of systematic distortions can be repre- 
sented within this formalism. Transformations like en- 
ergy scale and position offset have already been men- 
tioned. A Gaussian resolution systematic can be fioated 
by transforming the bandwidth, h, through analytic con- 
volution. Finally, re-weighting systematics, such as the 
neutron capture efficiency, are represented by varying the 
weight of events in the sum. 

The main challenge in using kernel estimation with 
large data sets is the computational overhead associated 
with repeatedly re-evaluating the PDFs as the parame- 
ters associated with detector response vary. We made 
several algorithmic improvements to make kernel estima- 
tion more efficient and did much of the calculation on 
off-the-shelf 3D graphics processors. For more detail on 
the implementation of the fit on the graphics processors, 
see [41| . 

The kernel-estimated PDFs had the same dimensional- 
ity over the same ranges of the observables as the binned 
fit, except with an upper energy limit of 20 MeV instead 
of 11.5 MeV. CC rates were extracted in 0.5 MeV in- 
tervals up to 12 MeV, with a large 12-20 MeV interval 
at the end of the spectrum. To reduce the number of 
free parameters in the fit, ES rates were extracted in a 
3.5-4.0 MeV interval, in 1 MeV intervals from 4 MeV 
to 12 MeV, and in a final 12-20 MeV interval. The CC 
and ES PDFs were fixed to be flat in the Toff dimension 
within each T^s interval. During fitting, the following 
parameters, corresponding to the dominant systematic 
uncertainties, were allowed to vary continuously: 

• Energy scale (both correlated and uncorrelated in 
each phase) 

• Energy resolution (uncorrelated in each phase) 
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• /3i4 electron and neutron scales 

• PMT /3-7 exponent (uncorrelated in each phase) 

• PMT j3-j offset (uncorrelated in each phase) 

• PMT /3-7 Teff exponent (uncorrelated in each 
phase) 

Altogether there were 18 CC parameters, 10 ES param- 
eters, 1 NC parameter, 17 background normalization pa- 
rameters, and 16 detector systematic parameters. The 
remaining systematic uncertainties were applied using 
the 'shift-and-refit' approach. 



14.4. Energy-Dependent i^e Survival Probability Fit 
Using Kernel Estimation 

The unique combination of CC, ES, and NC reactions 
detected by SNO allowed us to fit directly for the energy- 
dependent Ve survival probability without any reference 
to flux models or other experiments. Such a fit has sev- 
eral advantages over fitting for the neutrino mixing pa- 
rameters using the NC rate and the 'unconstrained' CC 
and ES spectra described in the previous sections. 

The unconstrained fits described in Sees. 114.21 and [T473l 
produce neutrino signal rates for CC and ES in intervals 
of reconstructed energy, Toff, with the free parameters in 
the fit directly related to event counts in each T^ff in- 
terval. Although this simplifies implementation of the 
signal extraction fit, physically-relevant quantities, such 
as total ®B neutrino fiux and neutrino energy spectra, 
are entangled with the energy response of the SNO de- 
tector. Comparing the unconstrained fit to a particu- 
lar model therefore requires convolving a distorted ®B 
neutrino spectrum with the differential cross sections for 
the CC and ES interactions, and then further convolving 
the resulting electron energy spectra with the energy re- 
sponse of the SNO detector to obtain predictions for the 
Toff spectra. 

Moreover, the unconstrained fits of Sees. 114. ^ and flX^ 
have more degrees of freedom than are necessary to de- 
scribe the class of MSW distortions that are observable 
in the SNO detector. For example, the RMS width of 
Teff for a 10 MeV neutrino interacting via the CC pro- 
cess is nearly 1.5 MeV. Therefore, adjacent Toff bins in 
the unconstrained fit are correlated, but this information 
is not available to the minimization routine to constrain 
the space of possible spectra. By fitting for an energy- 
dependent survival probability, we enforce continuity of 
the energy spectrum and thereby reduce covariances with 
backgrounds, most notably ^^^Bi events. Events from the 
CC reaction can no longer easily mimic the steep expo- 
nential shape of the background energy distribution. In 
addition, systematic uncertainties that are correlated be- 
tween the CC and NC events will naturally cancel in this 
approach within the fit itself. 



We therefore performed a signal extraction fit in which 
the free parameters directly described the total *B neu- 
trino flux and the energy-dependent Ve survival proba- 
bilities. We made the following assumptions: 

• The observed CC and ES Toff spectra come from 
a fixed distribution of neutrino energies, E^,, with 
the standard differential cross sections; 

• The Ve survival probability can be described by 
a smooth, slowly varying function of E^, over the 
range of neutrino energies to which the SNO detec- 
tor is sensitive; 

• The CC, ES and NC rates are directly related 
through unitarity of the neutrino mixing matrix; 

• Ve regeneration in the Earth at night can be mod- 
eled as a linear perturbation to the daytime i^e sur- 
vival probability. 

Given these assumptions, we performed a fit in which 
the neutrino signal was described by six parameters: 

• - the total ^B neutrino flux; 

• Co, ci, C2 - coefficients in a quadratic expansion of 
the daytime Ve survival probability around E^ ~ 
10 MeV; 

• oq, ai - coefficients in a linear expansion of the 
day/night asymmetry around E^, = 10 MeV. 

The day/night asymmetry. A, daytime survival proba- 
bility, P^P' , and nighttime i^e survival probability, P^l^^^ , 
that correspond to these parameters are: 



A{E^) = ao + ai{E^ - 10 MeV) 

Pi^'^iE,) = co + ci(S, -lOMeV) 

^2 



C2{E^ - 10 MeV)' 
l + A{E,)/2 
1 - A{E,)/2 



(22) 

(23) 
(24) 



The survival probabilities were parameterized in this way 
to reduce correlations between cq and the higher order 
terms by expanding all functions around the detected 
^B spectrum peak near 10 MeV. The simulated neutrino 
energy spectrum after application of the analysis cuts, 
shown in Figure [26l rapidly drops in intensity away from 
10 MeV. The broad Toff resolution of the detector in 
combination with the limited range of detectable neu- 
trino energies limits our sensitivity to sharp distortions. 
For this reason, we chose to fit for a smooth, polyno- 
mial expansion of the survival probability. By using a 
generic form, we allow arbitrary models of neutrino prop- 
agation and interaction to be tested, including standard 
MSW effects, as long as they meet the assumptions de- 
scribed above. Monte Carlo studies demonstrated that 
this analytical form was sufficient to model the class of 
MSW distortions to which the SNO detector was sensi- 
tive. We propagated the uncertainty in the shape of the 
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undistortcd '^B energy speetrum as an additional 'shift- 
and-refit' systematic uncertainty to ensure the extracted 
survival probability incorporated this model dependence. 




Neutrino Energy (MeV) 



FIG. 26: (Color online) Simulation of the undistorted energy 
spectrum of *B neutrinos that trigger the detector, before the 
application of the T^b threshold, and after a T^s > 3.5 MeV 
cut is applied, normalized to the SSM prediction. The sharp 
cut in Toff results in a smooth roll-off in detection efficiency 
for energies less than the peak energy. Also shown is the 
spectrum of incident neutrinos predicted by [s^l, arbitrarily 
normalized, to illustrate the effect of detector sensitivity. 

To implement this fit, we performed a 'four phase' 
signal extraction, dividing the data and the PDFs into 
Phase I-day, Phase I-night, Phase Il-day, and Phase II- 
night groups. Background decay rates from solid me- 
dia, such as the acrylic vessel and the PMTs, were con- 
strained to be identical day and night. Decay rates in 
the D2O and II2O regions were free to vary between day 
and night to allow for day/night variations in the water 
circulation and filtration schedules. We floated the same 
detector systematics as in the unconstrained fit described 
in Sec. 114.31 The fit has 6 neutrino parameters, 26 back- 
ground normalization parameters, and 16 detector sys- 
tematic parameters, for a total of 48 free parameters. 

We constructed the PDFs in the same way as described 
in Sec. 114. 3[ with the exception of the CC and ES sig- 
nals. Instead of creating a 3D PDF (/3i4, R^, cosOq) for 
intervals in Tgff in the undistorted spectrum, we created 
4D PDFs (Teff, /3i4, R^, cos9q) for separate E^, inter- 
vals in the undistorted spectrum. There were 9 CC and 
9 ES PDFs in each of the 4 day/night phases, with Ei, 
boundaries at 4, 6, 7, 8, 9, 10, 11, 12, 13, and 15 MeV. 

During optimization, the signal rates associated with 
the 76 CC, ES and NC PDFs were not allowed to vary 
freely, but were determined by the 6 neutrino parameters. 
We defined an 'ES survival probability': 

P^^^E^) = Pi^y + e(l-P^,^y{E^)) (25) 
P^"f *(£;,) = P-^^^ + e{l PZ^''\E,)) (26) 

where e — 0.156 is the approximate ratio between the 



j/^ T- and Ve ES cross sections. The ES cross-section ratio 
is not constant as a function of neutrino energy, so we 
took the variation with energy as an additional system- 
atic uncertainty. The signal rates were defined in terms 
of ^SB, Pee and Pes to be: 

i?NC = (27) 

Rch = J, r Pte'iE.) (28) 

i^ii = r ^^-^ p^'^'^iE,) (29) 
= p '^^i r ^^^^^^^ (^0) 
- -F^T— r (^-) (31) 

where Eq is 4 MeV and Ei is the upper energy boundary 
of the i-th E^ interval. 

The survival probability fit included the same 'shift- 
and-refit' systematics as the unconstrained fit, along with 
all of the day /night systematics used in previous analy- 
ses 0, @ . These systematics accounted for diurnal varia- 
tions in reconstructed quantities, such as energy scale and 
vertex resolution, as well as long-term variation in detec- 
tor response which could alias into a day /night asym- 
metry. In addition, the non-uniformity of the cos 9q 
distributions of CC and ES events can also alias into 
a day/night asymmetry, so we incorporated additional 
day/night systematic uncertainties on all observables in 
the CC and ES PDFs. 



14.5. Application of Constraints 

A priori information from calibrations and background 
measurements was included in the fits to constrain some 
of the fit parameters, in particular several of the radioac- 
tive backgrounds (discussed in Sec. 113. 2p and any sys- 
tematic parameters floated in the flt. 

The extended likelihood function had the form: 

£(a, /3) = Ceiata{a\fi)CeaUb[P) (32) 

where a represents the set of signal parameters being fit 
for, j3 represents the nuisance parameters for the system- 
atic uncertainties that were floated in the fits, Cdata{a-\P) 
is the extended likelihood function for the neutrino data 
given the values of those parameters, and CcaiibiP) is a 
constraint term representing prior information on the sys- 
tematic parameters, obtained from calibration data and 
ex-situ measurements. The contribution to CcaUb{P) for 
each systematic parameter had the form: 

Ceamm = e "^^r' (33) 
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where Xi is the value of parameter i, and /i^ and ai are 
the estimated value and uncertainty determined from ex- 
ternal measurements (with asymmetric upper and lower 
values for ai where required). This results in a reduction 
of the likelihood as the parameter value moves away from 
the a priori estimate. 

14.6. Bias Testing 

To verify that the signal extraction methods were unbi- 
ased, we used half the Monte Carlo events to create 'fake 
data' sets, and the remaining events to create PDFs used 
in fits to the fake data sets. A fit was performed for each 
set and the results were averaged to evaluate bias and 
pull in the fit results. 

We created 100 sets containing only neutrino events, 
45 sets also containing internal background events, and 
15 sets containing the full complement of neutrino events 
and internal and external backgrounds. The numbers of 
fake data sets were limited by the available computing 
resources. 

The two signal extraction methods gave results that 
were in excellent agreement for every set. The biases for 
the neutrino fiuxes were consistent with zero, and the 
Gaussian pull distributions were consistent with a mean 
of zero and standard deviation of 1. 

Additional tests were performed in which one or more 
systematic shifts were applied to the event observables 
in the fake data sets, and the corresponding systematic 
parameters were floated in the fit, using a priori inputs 
as in the final signal extraction fits, to verify that the two 
independent methods for propagating systematic uncer- 
tainties were also unbiased. In all cases, the true values 
for the neutrino fluxes were recovered with biases consis- 
tent with zero. 

14.7. Corrections to PDFs 

A number of corrections were required to account for 
residual differences between data and PDFs derived by 
simulation. An offset of the laserball position along the z- 
axis during calibration of PMT timing introduced an off- 
set to reconstructed positions along this axis in the data. 
A correction was therefore applied to all data events, as 
described in Sec. HI In addition, a number of corrections 
were applied to the reconstructed energy and isotropy of 
events (see Sections 18.21 and M respectively). 

The Monte Carlo simulation was used to link the neu- 
trino rates between the two phases, thus taking into ac- 
count variations in detector efficiency and livetime. Sev- 
eral corrections were applied to the Monte Carlo flux pre- 
dictions, as described below. 

The predicted number of events for signal type i per 
unit of incident flux, including all correction factors, is: 

= TVf*^ <5"™ (5^^= TVf ° TVf TVf i?, r, (34) 



where: 

• iVj^^ is the number of events predicted by the 
Monte Carlo simulation for signal i per unit inci- 
dent flux. This is recalculated as needed to account 
for any systematic shifts applied to the PDFs. 

• (5*'™ corrects for events aborted in the simulation 
due to photon tracking errors. This correction in- 
creases with the number of photons in an event. 

• Sf'^'^ corrects for differences in the acceptances of 
the instrumental and high level cuts for data and 
Monte Carlo events fSec. fTOS)) . 

• Nl^° is a correction to account for CC interactions 
on chlorine and sodium nuclei in the D2O volume 
that are not modeled in the simulation. This cor- 
rection is relevant only to the CC signal in Phase II. 

• Nf is a correction to the number of target 
deuterons and hence is relevant to CC and NC only. 

• iVf is a correction to the number of target electrons 
and hence is relevant to ES only. 

• Ri accounts for radiative corrections to the 
neutrino-deuteron interaction cross section for NC. 
Radiative corrections relevant to the CC and ES 
interactions were included in the simulation. 

• r corrects for deadtime introduced into the data set 
by the instrumental cuts. 

These corrections are summarized in Table IXVIIl 



Correction Phase CC ES NC 



^sim 


I, II 


(1.0 - 0.0006238 xTeft )- 


race 


I 


0.9924 


0.9924 


0.9924 


race 
Oi 


II 


0.9930 


0.9930 


0.9954 


Nr 


1— 1 


1.0002 








I, II 


1.0129 




1.0129 


m 


I, II 




1.0131 




Ri 


I, II 






0.977 


T 


I 


0.979 


0.979 


0.979 


T 


1— 1 


0.982 


0.982 


0.982 



TABLE XVII: Corrections applied to the expected number of 
CC, ES and NC events used in the signal extraction fits. 



15. RESULTS 

The detailed improvements made to this analysis, as 
described in previous sections, allow a more precise ex- 
traction of the neutrino flux parameters and, as a result. 
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of the MSW oscillation parameters. Results from the m> 
constrained fit are given in Sec. 115. 1] and from the energy- 
dependent fit to the Vf, survival probability in Sec. 115.21 
This new method for directly extracting the form of the 

survival probability from the signal extraction fit pro- 
duces results that are straightforward to interpret. A 
direct comparison can be made of the shape of the ex- 
tracted survival probability to model predictions, such as 
the LMA-predicted low-energy rise. 

Sec. 115.31 describes the measurements of the neutrino 
oscillation parameters. As has been observed in a number 
of recent publications (42| - |4^ , the different dependence of 
the survival probability on the mixing parameters Qyi 
and 013 between solar and reactor neutrino experiments 
means that a comparison of solar data to reactor antineu- 
trino data from the KamLAND experiment allows a limit 
to be placed on the value of sin^ ^13 . The new precision 
achieved with the LETA analysis in the measurement of 
tan^ 012 results in a better handle on the value of sin^ 6*13 
in such a three-flavor oscillation analysis. Results of this 
analysis are presented in Sec. 115.31 including a constraint 
on the value of sin^ 6*13. 

15.1. Unconstrained Fit 

Our measurement of the total flux of active ®B solar 
neutrinos, using the NC reaction ($nc) is found to be: 

• Binned-histogram method 

^binned ^ s.MQ t[J;i^g(stat) l^ j57(syst) X lO^cm-^s-i 

• Kernel estimation method 

<i>*^'c™-' =5.171 ±g:i^9(stat)10:i?2(syst) X 10*^cm-2s-i 

This represents l3'g% total uncertainty on the flux, 
which is more than a factor of two smaller than the best 
of previous SNO results. The statistical uncertainty has 
been reduced by nearly -\/2, to 3.1%. However, the largest 
improvement is in the magnitude of the systematic un- 
certainty, which has been reduced from 7.3% and 6.3% in 
previous analyses of Phase II and Phase III [oj data, 
respectively, to 2.4% (taking the average of the upper 
and lower values). 

Figure [17] shows a comparison of these results to those 
from previous analyses of SNO data. Note that the 
spectral shape used in the previous Phase I and 
Phase II analyses [1^ differs from that used here [38| . 
The bands represent the size of the systematic uncer- 
tainties on each measurement, thus illustrating the im- 
provements achieved with this analysis. 

Throughout this analysis, the quoted 'statistical' un- 
certainties represent the uncertainty due to statistics of 
all signals and backgrounds in the fit, with correlations 
between event types taken into account. Therefore, they 
include uncertainties in the separation of signal events 
from backgrounds in the fits. For example, the statisti- 
cal uncertainties on the quoted results for <I>nc include 
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FIG. 27: (Color online) Total *B neutrino flux results using 
the NC reaction from both unconstrained signal extraction fits 
in comparison to unconstrained fit results from previous SNO 
analyses. 'LETA 1' refers to the binned-histogram method 
and 'LETA 11' to the kernel estimation method. 



both the Poisson uncertainty in the number of NC events, 
and covariances with other event types. This is different 
from previous SNO analyses, in which the background 
events were not included in the signal extraction fits and 
any uncertainty in the level of background events was 
propagated as an additional systematic uncertainty. 

The two independent signal extraction fit techniques 
are in excellent agreement, both in the central NC flux 
value and in the magnitude of the uncertainties. The 
result from the binned-histogram method is quoted as 
the final unconstrained fit result for ease of comparison 
to previous analyses, which used a similar method for 
PDF creation. 

This result is in good agreement with the prediction 
from the BS05(OP) SSM of 5.69x 10^ cm^^ g-i [46], to 
within the theoretical uncertainty of ±16%. It is also 
in good agreement with the BS05(AGS,OP) model pre- 
diction of 4.51xl0^cm"^ s^^ ± 16% which was con- 
structed assuming a lower heavy-element abundance in 
the Sun's surface. 

The extracted CC and ES electron spectra from both 
signal extraction fits, in terms of the fraction of one 
unoscihated SSM, using the BS05(OP) model flux of 
5.69xl0^cm^^ s~^ [i^, are shown in Figure B51 An un- 
suppressed, undistorted spectrum would correspond to 
a flat line at 1.0. A greater suppression is observed for 
CC events than ES, since the ES spectrum includes some 
contribution from and Vr whereas CC is sensitive only 
to . Both spectra are consistent with the hypothesis of 
no distortion. The results from the two independent sig- 
nal extraction fits are again in excellent agreement for 
both the central fit values and the uncertainties. 

Figure dS] shows the CC electron spectrum extracted 
from the binned-histogram signal extraction fit with the 
errors separated into the contributions from statistical 
and systematic uncertainties. As for the NC flux re- 
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FIG. 28: (Color online) Extracted a) CC and b) ES electron 
spectra as a fraction of one unoscillated SSM (BS05(OP)), 
from both signal extraction fits, with total uncertainties. The 
final 12-20 MeV bin in the kernel estimation fit is plotted at 
the mean of the spectrum in that range. Both spectra are 
consistent with the hypothesis of no distortion (a flat line). 



suit, the uncertainties are dominated by those due to 
statistics (which includes the ability to distinguish sig- 
nal from background). This demonstrates the effect of 
the significant improvements made both in the determi- 
nation of the individual systematic uncertainties, as pre- 
sented in previous sections, and in the improved treat- 
ment of the dominant systematic uncertainties, whereby 
the self-consistency of the data itself was used to further 
constrain the allowed ranges of these parameters. It is 
worth noting that correlations between bins, which are 
not shown, tend to reduce the significance of any ob- 
served shape. Fitting to an undistorted spectrum (the 
flat line on Fig. gives a value of 21.52 for 15 de- 
grees of freedom, which is consistent with the hypothesis 
of no distortion. The prediction for the Teff spectrum 
for CC events taken from the best fit LMA point from a 
previous global analysis of solar data [9] is also overlaid 
on Fig. [29l The value of the fit of the extracted spec- 
trum to this prediction is 22.56 for 15 degrees of freedom, 
demonstrating that the data are also consistent with the 
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FIG. 29: (Color online) Extracted CC electron spectrum 
as a fraction of one unoscillated SSM (BS05(OP)) from the 
binned-histogram signal extraction fit, with the uncertain- 
ties separated into statistical (blue bars) and systematic (red 
band) contributions. The predictions for an undistorted spec- 
trum, and for the LMA point Am^i = 7.59 x lO"'^ eV^ 
and tan^ 612 = 0.468 (taken from a previous global so- 
lar -fKamLAND fit and floating the **B flux scale) are over- 
laid for comparison. 



The one-dimensional projections of the fits in each ob- 
servable parameter from the binned-histogram signal ex- 
traction are shown for each phase in Figures [301 a-nd 1311 
Of particular note is the clear ES peak observed in the 
cos 9q fits for both phases (Figs.lSCTc) andl5TTc)). demon- 
strating the extraction of ES events over the integrated 
energy spectrum, even with the low 3.5 MeV threshold. 
The error bars represent statistical uncertainties; system- 
atic uncertainties are not shown. Figure [32] shows the 
one-dimensional projection in Teff from Phase II (as in 
Fig. l3lTa)) but with the fitted contributions from indi- 
vidual signal types separated into six categories: CC, ES, 
and NC neutrino events, internal backgrounds (within 
the D2O volume), external backgrounds (in the AV, II2O, 
and PMTs) and hep neutrino events. 
The for the one-dimensional projections of the fit are 
given in Table IXVIIII These were evaluated using statis- 
tical uncertainties only and are, therefore, a conservative 
test of goodness-of-fit in the one-dimensional projections. 
In all dimensions, the final result is a good fit to the data. 

Table I XXIII in Appendix [X] shows the extracted num- 
ber of events for the neutrino fit parameters from the 
binned-histogram signal extraction fit, with total statis- 
tical plus systematic uncertainties. 

Table IXIXI shows the total number of background 
events extracted by each signal extraction in each phase, 
and a breakdown of the number of background neutron 
events occurring within each region of the detector. The 
two methods are in good agreement based on expecta- 
tions from studies of Monte Carlo-generated 'fake' data 
sets. For comparison, the total number of events in each 
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FIG. 30: (Color online) One dimensional projections of the 
fit in each observable parameter in Phase I, from the binned- 
histogram signal extraction. The panels show the fit projected 
onto (a) energy [Teg), (b) radius cubed (R^), (c) direction 
(cos^o), and (d) isotropy (/3i4). 



FIG. 31: (Color online) One dimensional projections of the 
fit in each observable parameter in Phase II, from the binned- 
histogram signal extraction. The panels show the fit projected 
onto (a) energy (Teff), (b) radius cubed (J?^), (c) direction 
(cos^o), and (d) isotropy (/3i4)- 



data set is also given (taken from Table |VT|. Due to the 
exponential shape of the energy spectra of most sources 
of background in this fit, the majority of the background 
events fit out in the lowest two bins in T^s, illustrating 
one of the major challenges of the low energy analysis. 

Tables IXXIIltiXXTVl in Appendix S show the effects 
of the individual systematic uncertainties on the ex- 
tracted NC rate, the CC rate in two energy intervals (4.0- 



4.5 MeV and 9.5-10.0 MeV) and the ES rate in the 3.5- 
4.0 MeV interval, all taken from the binned-histogram fit. 
The dominant source of uncertainty on the total neutrino 
fiux measured with the NC reaction is the neutron cap- 
ture uncertainty. Further significant contributions come 
from the Phase II energy resolution, the /3i4 scale for neu- 
tron capture events, the energy-dependent fiducial vol- 
ume, and the cut-acceptance uncertainties. 
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— ^ Data 

Total fit result 

— CC 

ES 

NC 

Internal bkg 
External bkg 

hep 






Phase I 


Phase II 


Background 


LETA I LETA II LETA I LETA II 


Total background events 


6148.9 


6129.8 


11735.0 11724.6 


D2O neutrons 


29.7 


34.0 


122.4 133.5 


AV neutrons 


214.9 


191.4 


295.7 303.4 


H2O neutrons 


9.9 


8.4 


27.7 26.3 


Total data events 


9337 


18228 



TABLE XIX: Number of background events extracted from 
the signal extraction fits for each method. 'LETA I' refers to 
the binned-histogram signal extraction, and 'LETA IF refers 
to the kernel estimation method. The total number of events 
in each data set is also given, taken from Table IVII 



FIG. 32: (Color) One dimensional projection of the fit in 
Teff in Phase II from the binned-histogram signal extraction, 
with the individual signals separated into the three neutrino 
interactions, internal backgrounds (within the D2O volume), 
external backgrounds (in the AV, H2O, and PMTs) and hep 
neutrino events. 



Phase Observable 


(data points) 


I Teff 


8.17 (16) 


COS ^0 


3.69 (8) 


p 


2.61 (5) 


/3l4 


20.99 (15) 


II Tefi 


13.64 (16) 


COS 9q 


3.07 (8) 


P 


2.98 (5) 


Pl4 


26.25 (15) 



TABLE XVIII: x values for the fit of the extracted signals 
from the binned-histogram signal extraction to the data set 
for one-dimensional projections in each of the four observ- 
ables, in each phase. These were evaluated using statistical 
uncertainties only. The number of data points used for the 
calculations are given afterwards in parentheses. Because 
these are one-dimensional projections of a fit in four observ- 
ables, the probability of obtaining these x^ values cannot be 
simply evaluated; these are simply quoted as a qualitative 
demonstration of goodness-of-fit. 



tainties on the shape of the CC electron spectrum are en- 
ergy resolution and the shape of the PMT /3-j PDF, par- 
ticularly as a function of Teff. The /3i4 scale for electron- 
like events is also a significant contributor. It is worth 
noting that the contribution from the fiducial volume un- 
certainty, which was significant in previous analyses [Sj], 
is now relatively small. 
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FIG. 33: (Color) Effect of systematic uncertainties on the 
extracted CC electron spectrum. The inset shows the same 
plot on a larger scale. 



Figure [33] shows the effects of several groups of sys- 
tematic uncertainties on the extracted CC electron spec- 
trum, taken from the binned-histogram fit. Four groups 
cover systematic effects that apply to the observables 
(Teff, cos 9q, and /3i4), in which the individual con- 
tributions are summed in quadrature (for example the 
Teff group includes the effect of energy scale, resolution 
and linearity); 'normalization' uncertainties include neu- 
tron capture, cut-acceptance, energy-dependent fiducial 
volume and photodisintegration uncertainties; the final 
group consists of uncertainties in the shape of the PMT 
/3-7 PDF. The dominant sources of the systematic uncer- 



The two signal extraction methods are in excellent 
agreement for all the neutrino flux parameters, as well 
as the sources of background events. This is a stringent 
test of the result, since the two methods differed in sev- 
eral fundamental ways: 

• Formation of the PDFs 

The methods used to create the PDFs were entirely 
independent: one using binned histograms, and 
the other using smooth, analytic, kernel-estimated 
PDFs. 

• Treatment of systematic uncertainties 
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The dominant systematics in the fits were 'floated' 
using different approaches: in the kernel method 
they were floated directly, whereas an iterative like- 
lihood scan was used in the binned-histogram ap- 
proach. 

• PMT /3-7 constraint 

In the binned-histogram method, a constraint on 
the total number of PMT events was implemented 
using a bifurcated analysis of the data (Sec. 113.51) . 
whereas no constraint was applied in the kernel 
method. 

That these independent approaches give such similar re- 
sults demonstrates the robust nature of the analysis and 
the final results. 



15.2. Survival Probability Fit 

Under the assumption of unitarity (for example, no os- 
cillations between active and sterile neutrinos), the NC, 
CC, and ES rates can be directly related. Based on this 
premise, a signal extraction fit was performed in which 
the free parameters directly described the total neu- 
trino flux and the survival probability. This fit there- 
fore produces a measure of the total flux of ^B neutrinos 
that naturally includes information from all three inter- 
action types. Applying this approach, the uncertainty 
on the flux was reduced in comparison to that from the 
unconstrained fit fSec llS.!]) . The total flux measured in 
this way ($sb) is found to be: 



SB 



5.046 



(stat) " 



-0.123 



(syst) X 10^ 



which represents l3g% total uncertainty. This is the 
most precise measurement of the total flux of ^B neutri- 
nos from the Sun ever reported. 

The survival probability was parameterized as a 
quadratic function in E^^ representing P^^^ , and a linear 
day/night asymmetry, as deflned in Eqs. (|22|) and (|23|) of 
Sec. 114.41 The best-flt polynomial parameter values and 
uncertainties are shown in Table IXXVIl and the correla- 
tion matrix is shown in Table IXXVIIl both presented in 
Appendix \K\ For all the extracted parameters, the total 
uncertainty is dominated by that due to statistics. 

Figure [M] shows the RMS spread in the best flt survival 
probabilities, P^^^{E^) and P^^^'^\E^), and day /night 
asymmetry, A{E^). The bands were computed by sam- 
pling the parameter space 1000 times, taking into account 
the parameter uncertainties and correlations. Overlaid 
on Fig. [331 are the predicted shapes of the day and night 
survival probabilities and the day /night asymmetry for 
the best-fit point from a previous global analysis of solar 
data [9|. 

The advantage of this direct parameterization for 
the survival probability is that model testing becomes 
straightforward. We can test the goodness-of-fit to an 
undistorted spectrum by setting ci = C2 = 0.0 in Eq. B51 



and we can test the goodness-of-fit to a model with no 
day/night asymmetry by setting = ai — 0.0 in Eg. 1221 
Requiring both simultaneously, we find a Ax^ = 1.94 for 
4 degrees of freedom, demonstrating that the extracted 
survival probabilities and day/night asymmetry are con- 
sistent with the hypothesis of no spectral distortion and 
no day /night asymmetry. For comparison, the A^^ value 
of the fit to the LMA point shown in Fig. [Ml is 3.9 for 4 
degrees of freedom, showing that the data are also con- 
sistent with LMA. 

^0.6 
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FIG. 34: (Color online) Best fit and RMS spread in the (a) 
Pee^{E^), (b) P^^'^^'-iE^), and (c) A{E„) functions. The sur- 
vival probabilities and day /night asymmetry for the LMA 
point Amai = 7.59 x 10"^ eV^ and tan^ 012 = 0.468, taken 
from a previous global solar+KamLAND fit 9] , are shown for 
comparison. 



This method for parameterizing the day/night asym- 
metry differs from previous SNO analyses, which quoted 
an asymmetry for each interaction type: 



A = 2 



i't'N + 4>d) ' 



(35) 



where (pD and 4>iq are the interaction rates measured for 
the day and night data sets, respectively. A combined 
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FIG. 35: (Color online) One dimensional projections of the 
fit in Phase I-day, from the polynomial survival probability 
fit. The panels show the fit projected onto (a) energy (Teff), 
(b) radius cubed (R^), (c) direction (cos 9q), and (d) isotropy 
(/3i4). The binning of data is purely for display purposes; the 
fits were performed unbinned. 



FIG. 36: (Color online) One dimensional projections of the 
fit in Phase Il-night, from the polynomial survival probability 
fit. The panels show the fit projected onto (a) energy (Teff), 
(b) radius cubed (R^), (c) direction (cosS©), and (d) isotropy 
(/3i4). The binning of data is purely for display purposes; the 
fits were performed unbinned. 



analysis of the results from Phase I and Phase II, assum- 
ing an undistorted neutrino spectrum, gave a result of 
A — 0.037 ±0.040 ^8j. For comparison, the current anal- 
ysis made no assumption about the shape of the underly- 
ing neutrino spectrum, except that it is a smooth, slowly 
varying function of over the range of neutrino energies 
to which the SNO detector is sensitive. The value of qq 



extracted under this assumption was — 0.032 ± 0.040. 
Uncertainty on the day/night asymmetry measurement 
has always been dominated by statistics, so the improve- 
ments made to systematic uncertainties in this analysis 
have a small effect. The effect of the additional statistics 
gained by going lower in energy appears to be balanced 
by the additional degrees of freedom allowed in the shape 
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of the neutrino energy spectrum. 

The one-dimensional projections of the fits in the ob- 
servable parameters for Phase I-day and Phase Il-night 
are shown in Figures ES] and EH 

15.3. Mixing Parameters 

A three-flavor, active solar neutrino oscillation model 
has four parameters: 9i2 and ^13, which quantify the 
strength of the mixing between flavor and mass eigen- 
states, and Ani^i and Ato|]^, the differences between the 
squares of the masses of the neutrino propagation eigen- 
states. The approximation of Am|]^ ~ Am|2 can be 
made because |Am32| ^ |Ato2]^|, while the remaining 
mixing angle, 623, and the CP-violating phase, 5, are ir- 
relevant for the oscillation analysis of solar neutrino data. 

For the sake of comparison with other oscillation anal- 
yses, this work employed tan^ 612 to quantify the leading 
effects of the mixing angles for solar neutrino oscillations. 
Smaller effects due to ^13 are quantified with sin^6'i3. 
The value of Am§i was fixed to -^2.3 x lO'^ eV^ (471 . 
an assumption that was necessary for the numerical de- 
termination of the three-flavor survival probabilities, but 
whose precise value had very little impact on our calcu- 
lation. 

The parameters describing the Pee{E^) function for so- 
lar neutrinos are, in order of importance, ^12, Am^i, 6*13, 
and Am^i. For experiments sensitive to neutrinos from 
terrestrial sources, near the detector, the survival proba- 
bilities were accurately calculated using a formula with- 
out the effect of matter. The inclusion of matter effects in 
the survival probability calculation for solar neutrino ex- 
periments involves the numerical integration of a system 
of coupled differential equations: 

i-^i>a{x) = HfTpa{x) , (36) 

ax 

where Hf is the Hamiltonian in flavor space, including 
matter effects in both the Sun and the Earth, x is the 
position along the propagation direction, and ipa{x) is 
a vector containing the real and imaginary coefficients 
of the wave function, where a = (e, /x, r). The system 
was solved for each new value of x as the wave func- 
tion was propagated from the Sun to a given detector on 
the Earth. The probabilities were then calculated from 
the magnitudes of the wave function coefficients. The 
integration was performed with the adaptative Runge- 
Kutta algorithm. Radial profiles of the electron density 
and neutrino production in the Sun were taken from the 
BS05(OP) model (46)]. The matter density inside the 
Earth was taken from the Preliminary Reference Earth 
Model (PREM) 'is'], which is the most widely accepted 
data since the density profile is inferred from seismologi- 
cal considerations. For more details on the survival prob- 
ability calculation, see [4Qj. 

Constraints on neutrino mixing parameters can be de- 
rived by comparing neutrino oscillation model predic- 



tions with experimental data, as has been done in pre- 
vious SNO analyses 0-i|- The approach for the inter- 
pretation of the solar and reactor neutrino data used the 
covariance method. From a series of observables with 
an associated set of measured parameters from a num- 
ber of experiments, the corresponding theoretical expec- 
tations were calculated for a given neutrino oscillation 
parameter hypothesis. In order to calculate the model 
prediction for the neutrino yield at a given detector, each 
of the neutrino fluxes that the detector was sensitive to 
was weighted with the neutrino survival probabilities, 
convolved with the cross-sections for the neutrino-target 
interactions as well as with the detector response func- 
tion, and then considered above the experiment's energy 
threshold. The function quantifies the difference be- 
tween the experimental data and theoretical model ex- 
pectation for the observable under study. 

In the results presented here, the free parameters were 
the neutrino mixing parameters and the total fiux of the 
*B and hep neutrinos. The survival probabilities and, 
hence, the fluxes and spectra of solar neutrinos and re- 
actor antineutrinos were fully constrained by the mixing 
parameters. The function in each case was minimized 
over a fine grid of points with respect to tan^ 9i2, sin^ 6'i3, 
and Am2i- The Ax^ = ~ Xmin differences were the 
indicators of the confidence levels (C.L.) in the one- and 
two-dimensional projections. The 68%, 95%, and 99.78% 
C.L. regions in two-dimensional parameter projections 
were drawn following the standard definitions: Ax^ = 
2.279, 5.99, and 11.83, respectively. For one-dimensional 
projections the errors on the parameter were the stan- 
dard Icr C.L. at Ax^ = 1. For all projections shown 
in this section, the ^^-s minimized with respect to the 
undisplayed parameters at each point in the MSW space. 

The information from the LETA survival probability 
measurement was included by evaluating the polynomial 
survival probability and day /night asymmetry (as de- 
fined in Eqs. ^ and ^ of Sec. HAA^ that best repre- 
sented the model prediction at each point in the MSW 
plane. To do this, it was necessary to take into account 
the sensitivity of the SNO detector (including effects such 
as the energy dependence of the cross sections, reaction 
thresholds, and analysis cuts) so that the parameteriza- 
tion of the model prediction at each point in the MSW 
plane sampled the neutrino energy spectrum in the same 
manner and over the same range as the data. We calcu- 
lated the number of detected events that passed all the 
cuts as a function of neutrino energy using the Monte 
Carlo simulation, and what was thus equivalent to a 'de- 
tected neutrino energy spectrum' (given in Table lXX Villi 
in Appendix \^ was distorted by the model-predicted 
survival probability at each point in the MSW plane. 
This was fit to a similarly obtained spectrum, now dis- 
torted by the polynomial parameterization, allowing the 
five polynomial parameters to vary in the fit. At each 
point in the plane, we then calculated the x^ value of the 
fit of the model-predicted polynomial parameters (cq, ci, 
C2, ao, and ai) to the result from the signal extraction. 
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taking into account all uncertainties and correlations as 
output by the signal extraction fit. The SNO rates from 
Phase III '9] were treated as a separate data set. 

Figure [37] shows the allowed regions of the 
(tan^ 012, Am^i) parameter space when the LETA data 
were analyzed in combination with the rates from 
Phase III Q- The 21/ contours were projected from the 
parameter space at a constant value of sin^ — 0.0, 
making them equivalent to an effective two-flavor anal- 
ysis. While the best fit point falls in the so-called 
'LOW region, with Am^^ = 1.15 x 10 "'^(eV^) and 
tan^ 612 = 0.437 l^osg, the significance levels of the LOW 
and the higher mass Large Mixing Angle (LMA) regions 
are very similar. The predicted shape for the survival 
probability is very flat in both regions, and the day /night 
asymmetry is expected to be small, so the SNO-only anal- 
ysis has little handle on distinguishing the two regions. 
A notable difference between LOW and LMA is in the 
predicted sign of the slope of the energy dependence of 
the day/night asymmetry, with LOW predicting a neg- 
ative slope, as was extracted in the polynomial survival 
probability signal extraction fit reported in Sec. 115.21 
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FIG. 37: (Color) SNO (LETA + Phase III) two-flavor oscil- 
lation parameter analysis. 



As described above, the observables from the SNO 
LETA fit used in the hypothesis testing were the poly- 
nomial parameters of the survival probability. In a full 
global analysis, event yields were used for the other so- 
lar neutrino experiments, including the SNO Phase III 
results. For each set of parameters, the oscillation model 
was used to predict the rates in the Chlorine [l|. Gal- 
lium [2I, and Borexino [lol | experiments, the Super- 
Kamiokande Phase I zenith spectra [13] and Phase II 



day /night spectra [6|, and the KamLAND rates and spec- 
trum [iJI, as well as the SNO rates and spectra. The 
expected rates and spectra were divided by the respec- 
tive predictions, calculated without oscillations, to re- 
move the effects of the model scaling factors. The unitless 
rates were then used in the global calculation. 

Although the <I>8b scale was determined in the LETA 
signal extraction, we re-introduced it as a free parameter 
in the minimization at each point in the parameter 
space to constrain it with all solar data. The uncertainty 
of the scale was retrieved from its marginal distribution, 
as was done for the oscillation parameters. 

The SNO LETA covariance matrix was taken from the 
signal extraction output given in Table IXX VIII as before. 
For other experiments, the total covariance matrix was 
assembled from the individual statistical and systematic 
components, as described in [8(. Correlations between 
SNO's LETA and other solar experimental results were 
allowed via the floated $83 scale parameter. 

The KamLAND rates and spectrum were predicted us- 
ing three-flavor vacuum oscillations. Publicly available 
information about the KamLAND detector and nearby 
reactors were included in our calculation, which repro- 
duced the unoscillated spectrum of Fig. 1 of Ref. [l3| 
with good accuracy. To include the effects of three-flavor 
oscillations, we then compared the obtained with non- 
zero values of 61^ with those obtained with ^13 = 0, for 
each set of (tan^ 6112, Atoji) values. In this way, we built 
a Ax^ function to parameterize the change of the x^ 
map in Fig. 2 of Ref. [l^ due to a non-zero value of 
6*13. This allowed us to include the KamLAND experi- 
ment in our three-flavor neutrino oscillation analysis and 
to precisely reproduce KamLAND 's two- flavor neutrino 
contours. When including the KamLAND antineutrino 
spectrum we assumed CPT invariance, and we used the 
KamLAND data only to constrain the oscillation param- 
eters (as opposed to the flux scale), whereas all other 
solar neutrino rates were used to collectively determine 
the absolute scale of the ^B neutrino flux as well as the 
oscillation parameters. 

Figure [35] shows the allowed regions of the 
(tan^ 012, Amji) parameter space when the global solar 
data and the KamLAND data were analyzed, both sep- 
arately and together, in a two-flavor analysis. It is in- 
teresting to note that the global solar analysis does not 
significantly alter the constraints in the LMA region rel- 
ative to the SNO-only analysis. 

Figure [39] shows the results of a three-flavor oscil- 
lation analysis. Fig. [5Wa) shows an overlay of the 
global solar and the KamLAND allowed regions in 
(tan^ 012, ATO21) parameter space, under a two-flavor hy- 
pothesis. Fig. lSW b) shows the same overlay for the three- 
flavor hypothesis. Allowing the value of sin^ 0i3 to be 
non-zero clearly brings the two regions into much better 
agreement. The three-flavor contours show the effect of 
allowing both <I>sb and sin^ 6*13 to float at each point in 
space. Allowing these extra degrees of freedom worsens 
the uncertainties on the two dominant oscillation param- 
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FIG. 38: (Color) Two-flavor oscillation parameter analysis for 
a) global solar data and b) global solar + KamLAND data. 
The solar data includes: SNO's LETA survival probability 
day/night curves; SNO Phase III integral rates; CI; SAGE; 
Gallex/GNO; Borexino; SK-I zenith and SK-II day/night 
spectra. 



FIG. 39: (Color) Solar and KamLAND oscillation parameter 
analysis for a) a two-flavor oscillation hypothesis and b) a 
three-flavor hypothesis. The solar data includes SNO's LETA 
survival probability day/night curves, SNO Phase III integral 
rates, CI, SAGE, Gallex/GNO, Borexino, SK-I zenith and 
SK-II day/night spectra. The is minimized with respect 
to all undisplayed parameters, including sin^ &13 and "I>8b- 



eters, tan^ 6*12 and Am^i. The regions obtained with all 
solar data are consistent with the SNO-only data and 
show an extension of the space towards larger values of 
tan^ di2 when sin^ ^13 is allowed to vary. In contrast, 
the three-flavor KamLAND contours show an extension 
towards smaller values of tan^ 0i2. 

Figure 1301 shows the confidence regions in the 
(tan^ 012, sin^ ^^l3) space. The directionality of the 
contours explains the excellent agreement of tan^ 612 
between the solar and KamLAND experiments when 
sin^ 613 is allowed to vary in the fit. 

Tables IXXI and IXXII summarize the oscillation param- 
eter results from the various two- and three-flavor os- 
cillation analyses, respectively. When all solar experi- 
ments are combined with data from the KamLAND re- 



actor antineutrino experiment in a two- flavor fit, the best 
fit point is found to be at 9i2 = 34.06 ^o g4 degrees and 
Am|i = 7.59 loJ? x 10"^ eV^. The uncertainty on the 
mixing angle has been noticeably reduced in compari- 
son to SNO's previous analyses, resulting in the world's 
best measurement of 612 to date. The global value of 
<i>SB from this fit is extracted to a precision of ^2 95^- 
The combination with KamLAND in a three-flavor flt 
has allowed us to constrain sin^^is, giving a value of 
sin^ 6*13 = 2.00tlsl ^ 10^^. This implies an upper bound 
of sin2 6li3 < 0.057 (95% C.L.). 
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Oscillation analysis 
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6.31 X 10 


Solar+KamLAND 
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+0.042 
-0.033 
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Xmin 


/ndf 
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TABLE XXI: Best-fit neutrino oscillation parameters and ex- 
tracted *B flux from a three-flavor oscillation analysis. Un- 
certainties listed are ±la after the '^^^ minimized with 
respect to all other parameters. 



FIG. 40: (Color) Solar oscillation parameter analysis, identi- 
cal to Fig. l39f b). but projected in the mixing angle space. The 

including Amii and $ 



^ is minimized with respect to all undisplayed parameters, 



Oscillation analysis 


tan^ 6i2 


Amii(eV^) 


SNO (LOW) 


n. .07+0.058 
U.U/ -0.058 


-r -|r+0-38 
1- -0.18 


X 10"^ 


SNO (LMA) 


n. .(-7 +0.038 
-0.042 


5-50 _i 62 


X 10"^ 


Solar 


-0.041 


5.891^;}^ 


X 10"* 


Solar+KamLAND 


+0.040 
U.'iOl _o.029 


1 .oy _o.2i 


X 10"* 




Xmin/ndf 


$SB (xlO'^cm-^s-^) 


SNO (LOW) 


6.80/9 


t- Q-|o +0.176 
O.UlO _o.l99 


SNO (LMA) 


8.20/9 


4.984 


Solar 


67.5/89 


5.104 


Solar+KamLAND 


82.8/106 


= Q-|o +0.119 
O.UlO _o.l48 



TABLE XX: Best-fit neutrino oscillation parameters and ex- 
tracted *B flux from a two-flavor oscillation analysis. The 
'SNO' results are from the combined LETA + Phase III oscil- 
lation analysis. Uncertainties listed are ilo" after the was 
minimized with respect to all other parameters. 



16. SUMMARY AND CONCLUSIONS 

We have described here a joint low energy threshold 
analysis of SNO's Phase I and Phase II data sets down 
to an effective kinetic energy threshold of TeS = 3.5 MeV. 
The low threshold increased the statistics of the CC and 
ES events by roughly 30%, and of NC events by --70%. 
A new energy estimator improved the energy resolution 
by 6%, thus reducing the number of background events 
reconstructing above threshold by ~6Q%. Separation of 
electron-like and neutron capture events was improved 
by the joint fit of data from Phases I and II, due to 
the difference in neutron detection sensitivity in the two 
phases. In addition, use of calibration data to correct 



the Monte Carlo-generated PDF shapes, and reduction 
of systematic uncertainties, have all contributed to in- 
creased precision on both the total solar neutrino 
flux and the derived neutrino mixing parameters. Fit- 
ting our data without constraints on the shape of the 
underlying neutrino energy spectrum or the unitarity 
of the mixing matrix gives a total neutrino flux of 
0NC = 5.14^o;2o (stat©syst) x lO^cm"^ s"\ measured 
by the NC reaction only, where © refers to the quadra- 
ture sum. This is in good agreement with the predictions 
of recent Standard Solar Models. The uncertainties on 
this result are more than a factor of two better than in 
our previous publications. The CC and ES reconstructed 
electron spectra for this fit are consistent with the hy- 
pothesis of no spectral distortion, and with the best fit 
LMA point. 

We have also used the unique capabilities of the SNO 
detector to perform the first direct fit to data for the 
energy-dependent survival probability, without any 
reference to flux models or other experiments. The fit for 
the survival probability assumes unitarity of the neutrino 
mixing matrix, and that the underlying neutrino spec- 
trum follows a smoothly-distorted *B shape. We have 
parameterized the survival probability as a second-order 
polynomial, allowing for a linear energy-dependent asym- 
metry between day and night spectra. The fit gives us a 
total ^B neutrino fiux of = S.OStoio (stat © syst) x 
lO^cm^^ s^^. No evidence for either a significant spectral 
distortion or a day/night asymmetry was found. 

With the results of the survival probability fit, we have 
created contours that show the allowed regions of the 
mixing parameters, finding that for SNO data alone the 
best fit point is in the LOW region of parameter space, 
but consistent with the LMA region at the 68.3% con- 
fidence level. Combining all solar experiments and the 
KamLAND reactor antineutrino experiment in a two- 
flavor fit, we find the best fit point is at 9i2 = 34.06 ^J g4 
degrees and Amii = 7.59 t^JJ? x 10"^ eV^. The uncer- 
tainty on the mixing angle has been noticeably reduced 
from SNO's previous analyses, resulting in the world's 
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best measurement of Oi2- The global value of $8b from 
this fit was extracted to a precision 

oft^-38% In a three- 
flavor fit, we find sin^ 6*13 = 2.00;^;°^x IQ-^. This implies 
an upper bound of sin^ ^13 < 0.057 at the 95% confidence 
level. 
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Appendix A: SNO Results: Data Tables 



Table IXXIII shows the extracted numbers of neutrino 
events from the binned-histogram signal extraction fit in 
each phase. The total statistical plus systematic uncer- 
tainties are given. 

Tables IXXIIItiXXVl show the effects of the individual 
systematic uncertainties on the extracted NC rate, the 
CC rate in two energy intervals (4.0-4.5 MeV and 9.5- 
10.0 MeV) and the ES rate in the 3.5-4.0 MeV inter- 
val, taken from the binned-histogram unconstrained sig- 
nal extraction fit. 

The direct signal extraction fit to the Ve survival prob- 
ability parameterized the neutrino fluxes as: 

• <1>SB - the total neutrino flux; 

• Co, Ci, C2 - coefficients in a quadratic expansion of 
the daytime v^. survival probability around Ei, — 
10 MeV; 

• oo, ai - coefHcients in a linear expansion of the 
day/night asymmetry around E,y = 10 MeV. 



Parameter 


Number 


of Events 




Phase I 


Phase II 



CC 3.5-4.0 MeV 
CC 4.0-4.5 MeV 
CC 4.5-5.0 MeV 
CC 5.0-5.5 MeV 
CC 5.5-6.0 MeV 
CC 6.0-6.5 MeV 
CC 6.5-7.0 MeV 
CC 7.0-7.5 MeV 
CC 7.5-8.0 MeV 
CC 8.0-8.5 MeV 
CC 8.5-9.0 MeV 
CC 9.0-9.5 MeV 
CC 9.5-10.0 MeV 
CC 10.0-10.5 MeV 
CC 10.5-11.0 MeV 
CC 11.0-11.5 MeV 
ES 3.5-4.0 MeV 
ES 4.0-4.5 MeV 
ES 4.5-5.0 MeV 
ES 5.0-5.5 MeV 
ES 5.5-6.0 MeV 
ES 6.0-6.5 MeV 
ES 6.5-7.0 MeV 
ES 7.0-7.5 MeV 
ES 7.5-8.0 MeV 
ES 8.0-8.5 MeV 
ES 8.5-9.0 MeV 
ES 9.0-9.5 MeV 
ES 9.5-10.0 MeV 
ES 10.0-10.5 MeV 
ES 10.5-11.0 MeV 
ES 11.0-11.5 MeV 
NC 



147.00 : 
154.37; 
180.52; 
174.63 ; 
175.52; 



47+76.96 

69.98 tltll 

+ 17.26 
-16.42 
+ 17.35 
-16.53 
-16.93 
-16.59 
-14.99 
-14.94 
-13.53 
-13.63 

163.28 tllZ 
161.09 
142.23 till 
119.69 tin 
101.34 t^.^^ 
84.03 

82.49^,11 
58.75 tlfs 
25.90 till 

I'i.l-'J -20.76 

55.00 ti%r 

42.92 
35.90 tlil 

20.25 tt fr 
15.25 I5.73 
19.731^:^1 
23.971^:^1 
19.72 ttll 
18.75 Its? 
16 16+'*°i 

ID. ID _4 02 
11 47+3.49 

10.23 1^3^^ 
4.38 

q qy +2.26 
•J-J' -2.26 



-21.78 



[-108.36 
-107.09 

98.39 tlttl 

+24.16 
-22.98 
+24.27 
-23.11 



205.70 
215.89 
252.94 tllll 
244.55 ■ 
245.02 
227.47 tltil 
224.83 ti^:^5g 
198.41 ^ 
167.51 
142.44 



+21.00 
-20.92 
+ 18.88 
-19.03 



+ 13.57 
-13.92 
+ 12.05 
-12.40 
+ 10.89 
-11.29 
9.71 
-10.09 



+8.07 
-8.48 
+5.28 
-5.46 



116.39 tut 
83.36 
36.88 
104.30 tllil 
77.34 
60.32 
50.37 
28.331^:^* 



-14.54 
-14.04 
-10.72 
-10.72 
-10.22 
-10.21 



78 



21.33 t^.g? 
27.58 1^:™ 



33.69 



+7.46 



27.79 ttil 
26.54tt:lt 

+5.61 
-5.63 
+4.98 
-4.99 



22.65 

16.38 



14 64 

11. -4 50 

+3.72 
-3.72 



6.27 



fl.94 
-1.95 



4.831^3:^4 



5.18 
870.171^3:°^ 



7.44 
3257.04 



+2.79 
-2.80 
131.26 
124.61 



TABLE XXII: Extracted number of events for each neutrino 
parameter from the binned-histogram signal extraction fit, in 
each phase, with total uncertainties. 



Where the day/night asymmetry, A, daytime i>e survival 
probability, P^^^ , and nighttime i^e survival probability, 
PggSht^ that correspond to these parameters are: 



A{E,,) ^ ao + ai{E^ - 10 MeV) 
Pf^^[E,) = co-Hci(.B, -lOMeV) 
+C2{E^ - 10 MeVf 
1 + A{E,)/2 



1 - A{E,)/2 



(Al) 

(A2) 
(A3) 



The best-fit polynomial parameter values and uncer- 
tainties are shown in Table IXXVI| and the correlation 
matrix is shown in Table IXXVIII 
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Systematic Phase Effect on rate /% 









NC 


CCl 


CC12 


ESO 




scale (+) 


I, II 


-0.293 


-2.037 


-2.144 


-0.156 


'Teff 


scale (— ) 


I, II 

A, AA 


0.137 


0.475 


0.913 


0.035 


Teff 


scale (+) 


I 


0.030 


-0.956 


-0.337 


-0.148 


Tcff 


scale ( — ) 


I 


-0.084 


1.659 


0.652 


0.236 


Teff 


scale (+) 


II 


-0.307 


0.317 


-1.094 


0.105 


Teff 


scale ( — ) 


II 


0.177 


-0.493 


0.584 


-0.133 


Teff 


resn (elec) (+) 


I 


0.008 


-3.999 


-0.013 


-0.439 


Teff 


resn (elec) ( — ) 


I 


-0.030 


7.656 


0.017 


1.399 


Teff 


resn (elec) (+) 


II 


0.653 


-5.005 


-0.006 


-0.531 


Teff 


resn (elec) ( — ) 


II 


-0.716 


6.597 


0.027 


0.480 


Teff 


resn (neut) (+) 


I, II 


0.065 


-0.054 


-0.023 


-0.006 


Teff 


resn (neut) ( — ) 


I, II 


-0.041 


-0.058 


0.046 


0.013 


Teff 


linearity (+) 


I, II 


0.130 


-0.160 


0.379 


-0.125 


Teff 


linearity (~) 


I, II 


-0.132 


0.287 


-0.372 


0.301 


/3l4 


elec scale (+) 


I, II 


0.634 


-5.064 


-0.082 


-0.648 


I3l4, 


elec scale (— ) 


I, II 


-0.622 


5.559 


0.086 


0.607 


Pli 


neut scale (+) 


I, II 


0.719 


-1.962 


-0.040 


-0.068 


/3l4 


neut scale (— ) 


I, II 


-0.411 


1.204 


0.029 


0.048 


^14 


elec width (+) 


I, II 


0.306 


-1.263 


-0.079 


-0.027 


/3l4 


elec width (— ) 


I, II 


-0.286 


2.342 


0.058 


0.099 


^4 


neut width (+) 


I, II 


0.067 


-0.240 


-0.002 


-0.014 


Pl4 


neut width ( — ) 


I, II 


-0.054 


0.217 


0.012 


0.017 


Pli 


E-dep (+) 


I, II 


0.227 


1.661 


-0.054 


0.299 


Pl4 


E-dep (-) 


I, II 


-0.246 


-0.999 


0.068 


-0.228 



TABLE XXIII: Effect of systematic uncertainties in TcS and 
^14 on the NC rate, the CC rate in the intervals 4.0-4.5 MeV 
("CCl") and 9.5-10.0 MeV ("CC12"), and the ES rate in the 
interval 3.5-4.0 MeV ("ESO"). Systematics shown as apply- 
ing to both phases were treated as 100% correlated between 
the phases. The '(+)' and '( — )' labels refer to the result of 
applying the positive and negative side of each double-sided 
uncertainty. 



Table IXXVIIII lists the Monte Carlo-generated neu- 
trino energy spectrum for events that passed all the stan- 
dard analysis cuts (the "detected neutrino energy spec- 
trum" ) . Events are separated into those occurring during 
the daytime and during the nighttime. 



Systematic 


Phase 


Effect on 


rate /% 








NC 


CCl 


CC12 


ESO 


Angular resn (-I-) 


I 


-0.032 


-0.688 


-0.075 


1.176 


Angular resn ( — ) 


I 


0.039 


0.648 


0.128 - 


-1.477 


Angular resn (+) 




-0.058 


-0.458 


-0.172 


3.219 


Angular resn (— ) 


II 


0.065 


0.298 


0.194 - 


-3.488 


Axial scale (+) 


I 


-0.030 


0.261 


0.128 


0.047 


Axial scale ( — ) 


I 


0.188 


-2.377 


-0.746 - 


-1.344 


Axial scale (-f ) 


II 


0.030 


-0.366 


0.079 - 


-0.037 


Axial scale ( — ) 


II 


-0.320 


-1.981 


-0.493 - 


-0.892 


Z scale (-)-) 


I 


-0.052 


0.377 


0.151 


0.018 


Z scale (— ) 


I 


0.000 


0.000 


0.000 


0.000 


Z scale (-|-) 


II 


0.004 


0.007 


0.018 


0.044 


Z scale (— ) 


II 


-0.070 


-0.906 


-0.130 - 


-0.112 


X offset (-I-) 


I 


-0.002 


-0.075 


-0.010 - 


-0.444 


X offset (-) 


I 


0.004 


-0.103 


-0.000 


0.032 


X offset {+) 


II 


0.009 


-0.538 


0.009 - 


-0.075 


X offset (-) 


II 


-0.007 


0.002 


0.003 


0.022 


Y offset (+) 


I 


-0.035 


-0.034 


0.000 


0.009 


Y offset (-) 


I 


0.005 


-0.084 


0.002 - 


-0.101 


Y offset (-I-) 


II 


-0.029 


-0.695 


0.035 - 


-0.279 


Y offset (-) 


II 


0.003 


-0.146 


0.007 


0.046 


Z offset (-I-) 


I 


0.011 


-0.275 


-0.032 - 


-0.642 


Z offset (-) 


I 


-0.003 


-0.060 


0.002 


0.112 


Z offset (-I-) 


II 


—0.168 


— 1.009 


0.006 - 


-0.317 


7 nffipt ( — \ 


II 


-0.013 


0.027 






X resn 


I 


-0.002 


-0.206 


-0.004 - 


-0.216 


X resn 


II 


0.052 


-0.732 


0.003 - 


-0.020 


Y resn 


I 


-0.007 


0.079 


-0.002 - 


-0.109 


Y resn 


II 


0.038 


-0.417 


0.019 - 


-0.201 


Z resn 


I 


-0.003 


0.173 


-0.002 - 


-0.224 


Z resn 


II 


0.115 


-1.354 


0.023 - 


-0.418 



TABLE XXIV: Effect of systematic uncertainties in cos^q 
and on the NC rate, the CC rate in the intervals 4.0- 
4.5 MeV ("CCl") and 9.5-10.0 MeV ("CC12"), and the ES 
rate in the interval 3.5-4.0 MeV ("ESO"). The '(-I-)' and '(-)' 
labels refer to the result of applying the positive and negative 
side of each double-sided uncertainty. 
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Systematic 


Phase 


Effect on 
NC CCl 


rate /% 
CC12 


ESO 


E-dep fid vol (+) 


I 


0.397 


-0.277 


-1.735 


0.378 


E-dep fid vol (-) 


I 


-0.230 


0.119 


1.027 


-0.233 


E-dep fid vol (+) 


II 


-0.698 


0.794 


-1.144 


0.322 


E-dcp fid vol (-) 


II 


0.825 


-0.994 


1.376 


-0.389 


Cut acceptance (+) 


I, II 


-0.357 


-0.519 


-0.434 


-0.451 


Cut acceptance (— ) 


I, II 


1.039 


1.299 


1.136 


1.171 


Photodisint.n (+) 


I, II 


-0.180 


0.134 


-0.002 


0.026 


Photodisint.n (— ) 


I, II 


0.183 


-0.100 


0.004 


-0.023 


neut cap (+) 


I 


-0.049 


-0.797 


0.003 


-0.074 


neut cap (— ) 


I 


0.044 


0.829 


-0.001 


0.084 


ncut cap (+) 


II 


-1.306 


0.616 


-0.001 


0.062 


neut cap (— ) 


II 


1.338 


-0.612 


0.003 


-0.060 


neut cap (+) 


I, II 


-0.759 


0.040 


-0.000 


-0.001 


neut cap (— ) 


I, II 


0.770 


-0.053 


0.001 


-0.011 


^^Na model (+) 


II 


0.028 


-0.751 


0.008 


-0.056 


2*Na model (-) 


II 


0.067 


-0.463 


0.003 


-0.182 


PMT Teff exponent (+) 


I 


0.009 


-6.482 


-0.003 


-1.469 


PMT Teff exponent (-) 


I 


0.002 


3.217 


0.004 


0.821 


PMT Teff exponent (+) 


II 


0.046 


-0.814 


0.001 


-0.196 


PMT Teff exponent (-) 


II 


0.011 


-0.328 


0.003 


0.010 


PMT exponent (+) 


I 


-0.048 


-2.875 


0.003 


-0.402 


PMT 7?^ exponent (-) 


I 


0.035 


1.746 


0.000 


0.238 


PMT R-^ exponent (+) 


II 


0.023 


-2.371 


0.002 


-0.185 


PMT R^ exponent (-) 


II 


0.004 


0.870 


-0.000 


0.440 


PMT R^ offset (+) 


I 


0.053 


5.674 


-0.004 


0.774 


PMT R:^ offset (-) 


I 


-0.016 


-2.113 


0.003 


-0.203 


PMT R^ offset (+) 


II 


-0.005 


0.735 


-0.000 


0.370 


PMT R^ offset (-) 


II 


0.001 


-1.014 


0.003 


-0.111 


PMT /3i4 mean (+) 


I 


-0.042 


-2.271 


0.002 


-0.714 


PMT j3i4 mean (-) 


I 


0.062 


0.559 


0.000 


0.509 


PMT /3i4 mean (+) 


II 


-0.516 


4.456 


0.029 


0.396 


PMT Pi4 mean (-) 


II 


0.524 


-4.102 


-0.027 


-0.802 


PMT /3i4 width (+) 


I 


0.075 


-1.388 


-0.001 


-0.008 


PMT ^14 width (-) 


I 


-0.070 


0.192 


0.005 


0.060 


PMT /3i4 width (+) 


II 


0.357 


-1.054 


-0.006 


0.257 


PMT /3i4 width (-) 


II 


-0.365 


1.394 


0.009 


-0.459 



Parameter Value Stat Syst D/N Syst 



TABLE XXV: Effect of relative normalization uncertainties 
and systematic uncertainties in background PDFs on the NC 
rate, the CC rate in the intervals 4.0-4.5 MeV ("CCl") and 
9.5-10.0 MeV ("CC12"), and the ES rate in the interval 3.5- 
4.0 MeV ("ESO"). Systematics shown as applying to both 
phases were treated as 100% correlated between the phases. 
The '(+)' and '(— )' labels refer to the result of applying the 
positive and negative side of each double-sided uncertainty. 



ao 

ai 

CO 
Cl 
C2 



+0.0366 
-0.0360 
+0.0279 
-0.0292 
+0.0205 
-0.0197 



0.0325 

-0.0311 

0.3435 

n nn7Qc; +0.00780 +0.00308 

U.UUOO -0.00745 -0.00335 

(-0.00302 +0.00148 

-0.00311 -0.00128 



+0.0059 
-0.0092 
+0.0104 
-0.0056 
+0.0111 
-0.0066 



+0.0145 
-0.0148 
+0.0140 
-0.0129 
+0.0050 
-0.0059 
+0.00236 
-0.00240 
+0.00057 
-0.00074 



TABLE XXVI: Extracted polynomial parameter values, star- 
tistical uncertainties, average systematic uncertainties, and 
day /night systematic uncertainties from the survival probar 
bility fit. 





$8b 


ao 


ai 


Co 


Cl 


C2 




1.000 


-0.166 


0.051 


-0.408 


0.103 


-0.246 


ao 


-0.166 


1.000 


-0.109 


-0.263 


0.019 


-0.123 


ai 


0.051 


-0.109 


1.000 


-0.005 


-0.499 


-0.031 


Co 


-0.408 


-0.263 


-0.005 


1.000 


-0.101 


-0.321 


Cl 


0.103 


0.019 


-0.499 


-0.101 


1.000 


-0.067 


C2 


-0.2 Ki 


-0.12:! 


-().0:u 


-0.:!2i 


-O.OGT 


i.OOO 



TABLE XXVII: Correlation matrix for the polynomial sur- 
vival probability fit. 
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Energy (MeV) 


D 


ay 




Ni 


ght 




Energy (MeV) 


Day 


Night 


2.2 


7.717 


xlO" 


6 


7.726 


XlO" 


6 


9.0 


7.675 


XlO" 


1 


9.699 xlO" 


1 


2.6 


7.211 


xlO" 


5 


8.505 


xlO" 


5 


9.4 


7.858 


xlO" 


1 


9.970 XlO" 


1 


3.0 


5.074 


xlO" 


4 


6.592 


XlO" 


4 


9.8 


7.882 


xlO" 


1 


1.000 




3.4 


2.168 


xlO" 


3 


2.992 


XlO" 


3 


10.2 


7.666 


XlO" 


1 


9.723 XlO" 


1 


3.8 


7.339 


xlO" 


3 


8.796 


xlO" 


3 


10.6 


7.298 


xlO" 


1 


9.251 XlO" 


1 


4.2 


1.599 


xlO" 


2 


1.971 


XlO" 


2 


11.0 


6.725 


XlO" 


1 


8.524 XlO" 


1 


4.6 


3.165 


xlO" 


2 


3.948 


XlO" 


2 


11.4 


5.974 


XlO" 


1 


7.573 XlO" 


1 


5.0 


6.130 


xlO" 


2 


7.632 


xlO" 


2 


11.8 


5.117 


XlO" 


1 


6.485 xlO" 


1 


5.4 


1.099 


XlO" 


1 


1.375 


xlO" 


1 


12.2 


4.137 


XlO" 


1 


5.256 xlO" 


1 


5.8 


1.768 


XlO" 


1 


2.221 


xlO" 


1 


12.6 


3.167 


XlO" 


1 


4.000 XlO" 


1 


6.2 


2.595 


XlO" 


1 


3.266 


XlO" 


1 


13.0 


2.211 


XlO" 


1 


2.807 XlO" 


1 


6.6 


3.491 


XlO" 


1 


4.403 


xlO" 


1 


13.4 


1.368 


XlO" 


1 


1.748 XlO" 


1 


7.0 


4.398 


XlO" 


1 


5.560 


XlO" 


1 


13.8 


7.208 


XlO" 


2 


9.023 XlO" 


2 


7.4 


5.260 


XlO" 


1 


6.667 


XlO" 


1 


14.2 


2.965 


XlO" 


2 


3.786 XlO" 


2 


7.8 


6.061 


XlO" 


1 


7.713 


xlO" 


1 


14.6 


9.843 


XlO" 


3 


1.248 XlO" 


2 


8.2 


6.761 


XlO" 


1 


8.508 


xlO" 


1 


15.0 


2.799 


XlO" 


3 


3.578 XlO" 


3 


8.6 


7.275 


XlO" 


1 


9.243 


XlO" 


1 


15.4 


2.008 


XlO" 


4 


2.086 XlO" 


4 



TABLE XXVIII: Monte Carlo-generated undistortcd neutrino energy spectrum for events that passed all the applied 
analysis cuts, divided into those occurring during the daytime and during the nighttime. The spectra have been normalized to 

the peak nighttime response, and the relative scales of the day and night spectra reflect the livetime and detector acceptance 
differences between day and night. The quoted energies are the central values of 0.4 MeV intervals. The spectrum is zero 
outside the displayed range. 



